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XIX. BAND DRITTES HEFT 1951 


Die unendlich ausgedehnte Scheibe mit gleichformig am Rand belastetem 
quadratischen Loch *. 


Von H. Montag. 


1, Einleitung. Die vorliegende Arbeit behandelt das quadratische Loch in einer unendlich 
ausgedehnten Scheibe mit doppelsymmetrischer gleichférmiger Randbelastung. Die Grundlagen 
zu dieser Theorie sind von Th .Schade! geschaffen worden. Es kann daher darauf verzichtet 
werden, im einzelnen auf die Entwicklung einzugehen. Es sollen nur die wesentlichen hier be- 
nutzten Ergebnisse angegeben werden. 

Das Scheibenproblem fiihrt bekannterweise auf die Differentialgleichung 


AA He 0 
mit vorgegebenen Randbedingungen. Aus der Funktion F lassen sich dann die Spannungen 
berechnen zu 
meee _ OF eel ea 
Ou Oat ” Ty Ox dy” 


ta aye? 
Zur Anpassung des Koordinatensystems an die Randbedingung werden die von R. Grammel? 
-angegebenen erweiterten Kreisfunktionen 

x = Wo} (n) v, 

y= Sin (n) v, 
die der algebraischen Gleichung x" + y"=— 1  geniigen, 
benutzt. 

Mit der Transformation 

x= 960} (n) v, 

¥ =o Sin (n) v 
(Abb. 1) erhalt man mit dem Separationsansatz zur Lésung 
der Potentialgleichung A ® = 0 

D, = 0 & (0) 


die allgemeine Liésung 


Abb. 1. 


, = 0? (Gin (n)? v +-Coj (n)? v) Al eat are tg (Tg (n) v) x + yl — oP fir 9 = 0,5; 0,7; 1. 
und mit n—> oo nach einiger Umformung die bekannten Funktionen 

A= — <0: 04, = Fel, 

A— 0) : Dy = arc tg, aie 

A>0 ae 


| 
| Die entsprechenden Bipotentialfunktionen WY, ergeben sich durch Multiplikation der Potential- 
jfunktionen mit zz. Fiir das hier vorliegende Problem interessieren nur die Funktionen, die 
_keinen Pol im Unendlichen haben. 


2. Spezielle Bedingungen und Lésung der Differentialgleichung. Fiir das in Abb. 2 schematisch 


dargestellte Scheibenproblem miissen neben der Singularitatenfreiheit der Funktionen im Aufen- 


* Dissertation Aachen 1949. 
1 Th. Schade, Ing.-Arch. 19 (1951), S. 118. 
2 R. Grammel, Ing.-Arch. 16 (1948), 5S. 188 sowie 18 (1950), S. 250. 
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bereich des Loches die Spannungen und damit die zweiten Ableitungen der Funktionen im Un- | 
endlichen verschwinden. Weiterhin mu die Lésung zu v= 0; 2;4;6 gerade sein, was man 
aber durch Zusammensetzen von Teilfunktionen leicht erreichen kann (wie Th. Schade in der er- | 
wahnten Arbeit gezeigt hat). Man braucht sich dann nur noch um die Lésung im Bereich I | 
zu kiimmern und hat dafiir zu sorgen, daB die Stetigkeit der zweiten Ableitungen an den Dia- | 
gonalen erfiillt wird, was auch die Stetigkeit der ersten Ableitungen nach sich zieht. 
Diese Bedingung ist auch physikalisch sofort einzu- 
sehen. Die Diagonalen miissen aus Symmetriegriinden 
Hauptspannungstrajektorien sein, was letzten Endes 
eine Deformationsbedingung darstellt, die aussagt, daB | 
die Diagonallinien gerade bleiben. Das 1aBt sich aber | 
nur mit der Bedingung (0,),—1 = (Ov)»—1 erreichen, | 
was mit oben Gesagtem identisch ist. Wie man spater 
sehen wird, ist die Erfiillung dieser Bedingung streng | 
nicht méglich, so daB man zu einer Naherung ge-_ 
zwungen ist. 
Aus Symmetrieeigenschaften kann man weiterhin | 
fordern, da8 t,, = 0 fiir v= 0 ist. 
Alle erwahnten Bedingungen mit Ausnahme der | 
Stetigkeitsbedingung werden von denim vorigen Ab- | 
schnitt abgeleiteten Potential- und Bipotentialfunk- | 
Abb. 2. tionen erfillt, wenn man nur Funktionen zulaBt, | 
die gerade in v sind. 
Die Stetigkeitsbedingung soll hier durch eine Gleichgewichtsbedingung fiir einen Schnitt 
durch die positive x- und y-Achse verbessert werden. Sie lautet 


(fo, dx), <9 +f. ay )ean= in oa | 


Es hat sich ergeben, daB bei einem Ansatz fiir F, der sich aus einer Reihe von Polynomlésungen 
allein aufbaut, der erste Ausdruck von (2,1) zu Null wurde. Es mute daher noch die Bipotential- 
funktion 


yare tg = = Q Varc tg Vv 


zu einem Lésungsansatz hinzugenommen werden, welche a priori den angefiihrten Bedingungen 
einschlieBlich der Stetigkeitsbedingung geniigt. 
Der Lésungsansatz fiir die Bipotentialgleichung 


AA F0 
wiirde also lauten 
N N+1 
F = by Qo @ v are tee a Oo Oey Oe (2,2) 
= f= 
wofiir die Randbedingungen zu befriedigen sind 
er OF : : 
Tn Byte NOt) 30 Fay etre, = ONE ae (2,300 


und die Gleichgewichtsbedingung 


(fo, d x), =o 1 (fs. ay), Se Oy 


oder auch 


(fo. deena + (Joy 042) .-o= 02 (2,4) 


sowie die zusatzlichen Bedingungen der Stetigkeitsforderung an den Diagonalen. Die Rand- 
bedingung (2,3) wiirde bei einem Koeffizientenvergleich in den v-Potenzen ein Gleichungssystem 
ergeben, das auBerordentlich kompliziert ist und bei einer groBen Anzahl von Funktionen eine 
groBe Rechenarbeit erfordern wiirde. 
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Um diese Schwierigkeit zu umgehen, fithrt man eine neue Variable V durch die Substitution 


1 


ein. Die sich hiermit ergebenden Potential- und Bipotentialfunktionen sowie ihre Ableitungen, 
die sich als Linearkombinationen dieser schreiben lassen, sind in der Originalhandschrift tabel- 
larisch zusammengestellt. 
Die Randbelastung soll in der Form angenommen werden 
M 
(Flex eg—=— 2 Am V"*?, (2,6) 
m+l1 
Durch Koeffizientenvergleich in den V-Potenzen erhalt man das Koeffizientenschema nach 


Tabelle 1. Die Auflésung dieses Schemas ist in Tabelle 2 angegeben. 


Tabelle 1. 
ge 
es se er es pee EES EEE SE DS DS SE 
Rechte 
by | a, b, ComleGs iiss || Os ay, b, a; b, a b, a, by ag by Seite 
| 
—l i 4 3 0 
leant 6 6) Hepes 18 5 0 
—8\—16| 24| 64|/—20/—80| 5] 40 = 0 
_—24|—48] 80] 200/ —90—300| 42) 210| —7| —70 9) 0 
| —8s|—16| 30| 72/—42/—126| 28| 112| —9| —sal 0 
—5| —10| 2]! 491235) os 30) 105) | 0 
Dgh egiya) Vageer kero th 
4 so Si 0 
eg ee ee 
irs 
fey sy 642 zay 29 3 | —A,/2 
Bese 11D 2961-2151 915 i295 ang wy 
= Gils -16|. 40) 0881-—-60,/--160], «.35)..130 14-7 | —49 7 —A,/6 
—24|—48| 120| 264|—210|\—540| 168] 546|-63|—204| 9] 81—44/8 
—8|—16] 42) 92|84|—210| 84) 252 |45 |174—4,/80 
—5| 410) 28) “61\=63'-4154| 75) 213'—.4,/384 
—3| —6| 18| 39|—45|—108'—4,/1792 
—1| —14 - 97 —A,/2048 
= —2| O 
Tabelle 2. 
A,/4 | 4,/16 | 4,/96 | 4,/768 | 4,/2560 | .4,/30720 | A,/860160 | 4,/458752 
= sas Berea oe pee ees Cr —630 6930 —180180 —90090 
= 2 9 48 345 1050 11655 304920 153153 
= | —1 = S85 | 203 —315 —3465 —90090 | —45045 
= = as: = 156 —140 —1470 — 36960 —18018 
sb. — 1 is 45 140 1575 41580 21021 
os — nar —30 —90 —980 95200 | —12474 
_ 2 20 70 840 23100 12012 
ues —12 AG 588 —16128 —8316 
Pe 6 20 420 12600 6930 
16 —280 —8960 —5040 
7 8 168 6048 3696 
a —80 —3840 —2592 
SS 40 2240 1680 
he —960 —1008 
ee 480 516 
ee: = 294 
2, 112 
os 


Zur Erfiillung der Gleichgewichtsbedingung (2,4) hat man noch die Integrale 
0° 20 
d Pr ax) 
(i do) un (J Cem) 
zu bilden. 


| 11* 
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Setzt man in (2,2) bzw. (2,6) N= M, so zeigt die Rechnung fiir jedes beliebig N > 1, daB- 


(Fes de), are eo 


20 


ist. Zur Ausrechnung des zweiten Integrals bedient man sich der Rekursionsformel 


1 


dv 1 2m—3 | 
In= { ‘3 == ile a ] 
Pay” 21) 2 (1 + v2) 
0 


Tabelle 3. Die Ergebnisse sind in Tabelle 3 
See 7 7 7 zusammengestellt. | 
ha | Ss | oe | i“ | ime ie Damit ergibt sich aus der Gleich- 
i 0,7853982 4 | 0,4746036 7 | 0,3532164 gewichtsbedingung (2,4) fiir - 
2 0.6426991 5 0,4230906 8 0.3285447 Neel 28 
3 | 0,5445243 6 | 0,3839066 9 | 0,3082548 
10 0.2912381 die Bedingungsgleichung zwischen | 


den Koeffizienten A; | 


1,426991 A, + 1,695243 A, + 1,621040 A, +-1,496531 A, +1,378129 4, +1,276305 4, | 
- 1,190725 A, +1,118742 A,= 0. | 


| 


Weitere Bedingungsgleichungen kann man nun wahlweise zur Annéherung der Belastungs- | 
funktion und als weitere Stetigkeitsbedingungen fiir die Spannungen an den Diagonalen fordern. 


3. Numerische Rechnung. Um die giinstigste Naherung zu finden sind mehrere Versuchs- | 
rechnungen durchgefiihrt worden, deren Extrakt hier in Tabelle 4 zusammengefaBt ist. 


Der Vollstandigkeit halber seien zunachst aber noch die Formeln fiir die Spannungen an den 
Diagonallinien, an den Mittellinien und fiir die Randspannung o, (0 = @)) fiir N = 8 angegeben 
werden. Fiir die Spannungen an den Diagonallinien ergibt sich fiir py = 1 ( 


b 2 3 5 4/3 33 5/3 3 
(x)= 12 = (9) oe (7) by (*) E Le by ! (7) (pats by) | (*) i a bs) — 
615 7/15 eal 8/21 a 00\9 a 7 
Carita ealig liane st iastea) leet hn 
b 4 3 > 5 3 7 
(o5)v—10 = (4) 3° — (0) "ba i lies) (sath Glesa 
615 15 15 8/21 iii 9/7 11 
+(G) ab + (3) (sete) +(8) (ea) +(8) Ge Zh). 


Als wertvolles Ergebnis ware zu erwaihnen, daf fiir jeden Lésungsansatz dieser Art 


o 


(0x)y=12 = (6y)y=1p (3,2), 

0 = 00 2= 06 

wird (exakte Erfiillung der Spannungsstetigkeitsbedingung in der Ecke). Fiir die Spannungen | 
an den Mittellinien folgt fiir N — 6 


Ge ()2 ies (f) 4 = (fe a, +10 b,) =e a, +18 b,) — a (12 a, +288,) — 
= al (20 a, +40 b,) — Gl (30 a, +54 b,) — a A2 dg, 


S 


| 


ip = eh (2 a, +2 b,) eel (6 a, +6 By) el (12 a, 1126.) Gh (20 a, +20 b,) + 
a0 (7) G0 a; +30 b,) Flea 42 ag 
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Die Formel fiir die Randspannung lautet fiir N — 6 


(Fy)o—o5= V (2 by + 4 by) + V?( 


+ V3 (16 b, — 136 b, + 264 b, — 180 b, +40 b, - 


ibs 320 bar 42 b— 18 by 


6 a, + 6 a3) - 


8 a, 


48 ay + 60 ay 


Tabelle 4.. Tabellarische Ubersicht der gerechneten Faille. 


potesnseen Poti 
Fall Randbedingungen n den : kee 
| wSaT Diagonalinien) Ren 
ai 
ee == 1 = 0 F 
i ‘ ie li s keine 9 
oe 0 = Qo 
ts) 
o= 1/4 
Mo ee == bh it = : 
x ; oe, ae. keine 1133 
viol | 
Txy = 0 C= 
2) 
»= 1/5 | 
Of — tL V2 O—= 0 
= Bye 
7 | 
‘ cee eS keine 17 
al 
Of or dv = —1 
0 
00~% Aye is 
Ox pela 0 = @o 
toa) 
v= 1/4 
0, = — 1 p= iN pe 0 = 0 
v = 3/4 
3* Dill a, = 0 17 
aes 2 =} Co 
1 
0, f O24 dv = — 1 
(e=9, | 
On = — li v= 1/2 )9 = eo ay 0 
eae | ar 13 
i ee pase eo 5 Go 
00x ; 
—- — 
(as v= 1/2 Bre 
eee G2— 0y = 
a ee 7 i 13 
2% Oox o0a1f° payee 
Ox | ve 
Uo == OR 0 = 0 


20 a4) + 


+ V4 (96 bs— 624-by+ 1280 b;—1140 bg + 462, + 48 a,— 240 a,+ 360.0,— 210 a5-+ 42 a4) + 
+ V* (384.b,— 23686; +5280 b—5712b, +192 a,— 960 a, +1680 a,—1344a,) +- 
+ V® (1280 b, — 8000 b, + 19488 b, + 640 a, — 3360 a; +6720 ag) + 
+ V7 (3840 bg — 24960 b, +1920 a; — 10752 a) + 
L V8 (10752 b, + 5376 ag) . 
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In Abb. 3 und 4 sind die Randspannungen und die Spannungen in den Mittellinien fiir den Fall 2? 
und 2* graphisch wiedergegeben. Zum SchluB sei noch in Abb. 5 das sich nach Fall 2 ergebende:» 
Trajektorienbild aufgezeichnet. 

Die Ergebnisse der angefiihrten Versuchsrechnung sind in der folgenden Kurventafel (Abb. 6)) 
dargestellt. Neben der Belastungsfunktion ist die Abweichung der Spannungen an den Dia-) 
gonalen (¢,, 0y),—1/2 vom arithmetischen Mittel der Spannung an jedem Punkt aufgetragen., 


Py 


0, fire=0o (Rand)@) P 
Oy tir g=@o(Rand) @) 


Oy fir V=1 (Mitte/ schmitt) @y 
\_& tir V-1 (Mitte/schnitt) (2) 


Or tir V=1 (2) 


Abb. 3. Oy am Rand. Abb. 4. Spannungen im Mittelschnitt. 


Man beachte die groBen ZeichenmaBstabe und — 
den eingetragenen Spannungsmaximalwert von der 
Einheit. 

4, Diskussion der Ergebnisse. Vergleicht man | 
Abb. la, 2a, 3a, so erkennt man, da eine Hin- 
zunahme von Potentialfunktionen eine erhebliche | 
Verbesserung in der Annaherung der Belastungs- | 
funktion liefert (in Abb. 3a ist der Unterschied | 
graphisch nicht mehr zu erkennen). Die Stetigkeit | 
der Spannungen an der Diagonalen wird jedoch nur 
fiir gréBere o-Werte verbessert, wahrend nahe an 
0 = 0) sogar eine geringfiigige Verschlechterung 
auftritt. Das Gebiet der nennenswerten Abwei- 
chungen der Normalspannungen wird jedenfalls verkleinert. Doch scheint dieser Weg fiir die 
weitere Rechnung nicht der giinstigste zu sein. 

Eine wesentliche Verbesserung der Stetigkeit der Normalspannungen an der Diagonalen durch > 
direkte Anniherung ihrer Potenzreihen, d.h. Gleichsetzen gewisser Potenzglieder (in 0/09) der 
Polynome der Spannungen zeigt Abb. 3*. Da die Anzahl der Potentialfunktionen beibehalten 
wurde, tritt eine Verschlechterung der Randspannungsnaherung ein, die aber nicht wesentlich — 
ist. In dieser Methode liegt auch der Weg fiir weitere Rechnungen, um zu noch genaueren Resul- 
taten zu gelangen. 

Daf eine punktweise Erfiillung der Stetigkeitsbedingung nicht zu empfehlen ist, zeigt Abb. 2* 
und 2**, Randbelastungs- und Spannungsdifferenzen sind gegeniiber denen der bisher besproche- 
nen Methoden sehr grof. Als Ergingung ist in Abb. 2** eine Randbelastungsfunktion angenommen 
worden, die in der Ecke verschwindet, wobei jedoch zwei Punktvorschriften fiir die Stetigkeit 
an der Diagonalen vorgegeben worden sind. Es zeigt sich etwa das gleiche Bild, abgesehen da- 
von, daB durch Punktvorschriften an der Diagonalen der Rechenaufwand merklich erhéht wird. 


Abb. 5. Trajektorienfeld. 


5. SchluBbemerkung. Die Anwendung der Methode der Polynomdarstellung fiir Lésungen 
des Scheibenproblems der vorliegenden Art hat gezeigt, daB man mit verhdltnismaBig wenig 
Rechenarbeit zu brauchbaren Ergebnissen gelangt. Die Gleichungssysteme zur Erfiillung der 
Randbedingungen lassen sich sukzessive leicht lésen, so daB die einmal durchgefiihrte Auflésung 
auch fiir andere doppelsymmetrische Randbelastungen bestehen bleibt. Da die Stetigkeit der 
Spannungen in den Ecken a priori erfiillt wird, besteht nur noch die Aufgabe darin, ein ver- 
niinftiges Verhaltnis zwischen der Annaherung der Belastungsfunktion und der Stetigkeit der 
Spannungen auferhalb der Ecken zu finden, was hier numerisch durchgefiihrt wurde. 
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Randbelastung 


if 
| 
| 
| 


Abweichung 
graphisch nicht 
mehr darstel/oar 


3 5 


Omax =r] 10 
5 
e~ se) 
v=17 b 


5 I 


Abb. 6. Graphische Darstellung. 


1 
a) der Randbelastung und b) der Spannungsabweichung a ee dy an der Diagonalen y =1 bei den durchgerechneten Beispielen. 


Institut fiir Mechanik der Technischen Hochschule Aachen. 
(Eingegangen am 26. September 1950.) 
Anschrift des Verfassers: Dr. H. Montag, Aachen, Technische Hochschule, Lehrstuhl fiir Mechanik. 
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Der eingespannte, achsial pulsierend belastete Stab als Stabilitatsproblem*. | 
Von F, Weidenhammer, 


1, Aufgabenstellung, Die vorliegende Arbeit untersucht das Stabilitatsverhalten eines achsial| 
pulsierend belasteten Stabes mit eingespannten Enden. Der freiaufgestiitzt gelagerte und in) 
gleicher Weise achsial belastete Stab mit seinem merkwiirdigen Resonanzverhalten ist mehrfach || 
ausfiihrlich untersucht worden»*? und hat Anla8 zu einer sehr allgemeinen Theorie der Stabilitat | 
der elastischen Bewegung*® gegeben. Dabei zeigt sich, da der achsial pulsierend belastete Stab | 
mit eingespannten Enden etwa der einfachste (weil eindimensionale) Fall ist, welcher den vollen | 
Umfang der allgemeinen Theorie erfordert, wahrend der ebenso belastete, jedoch freiaufgestiitzt | 
gelagerte Stab als leichter zuganglicher Sonderfall in der allgemeinen Theorie enthalten ist. 

Die vorliegende Untersuchung will das Stabilitatsverhalten der durch eine pulsierende Achsial- | 
last angeregten Langsschwingungen untersuchen. Zu dem Zweck sollen zusammengehérige Werte 
von Erregerfrequenz und Amplitude der pulsierenden Achsiallast gefunden werden, welche bei 
beliebigen Anfangsbedingungen zu unbegrenzt anwachsenden Querschwingungen des Stabes_ 
fiihren kénnen. Hierzu wird die ,,Allgemeine Theorie der Stabilitat erzwungener Schwingungen | 
elastischer Kérper‘‘ von E. Mettler? verwandt. Durch eine etwas andere mathematische Behand- | 
lung werden wesentliche Ergebnisse dieser Theorie auf anderem Wege bestatigt und durch eine 
unmittelbare Berechnung des fiir die Instabilitat maBgebenden charakteristischen Exponenten 
erganzt. Auf diesem Wege wurde auch die Einarbeitung einer geschwindigkeitsproportionalen 
Dampfung mdglich. 


2, Ausgangsgleichungen, Bei der systematischen Behandlung des Problems geht man vom 
Hamiltonschen Prinzip aus. Dieses Prinzip verlangt 


by 


[(A+V—T) dt= Extr., (1) 
ty 
worin mit den in der technischen Biegelehre iiblichen Vernachlassigungen 
EF y2 o\2 E 1/2 1/2 
ay re | (u.+"3) dx +72 f widx+P(0 u(3)— ff (ug tw?) dx (2) 
Ip. 1/2 —1/2 


zu setzen ist®, Hierin bedeuten (vgl. Abb. 1): 


z, x = die Ortskoordinaten, 


u (x,t), w (x,t) = die gesuchten Verschiebungskomponenten in Stablangsrichtung und quer zur 
Stabachse, 
t= die Zeit, 
J = Const. = das fiir die Biegeschwingung maBgebende 4 
Flachentragheitsmoment, Zz 
F = Const. = die Querschnittflaiche des Stabes, Pit) 
(4 = Const. = die Stabmasse je Langeneinheit, walt) = 
l= die Stablange, 
= den Elastizitatsmodul, Abb. 1. Der Stab im Koordinatensystem. 
P(t)= die pulsierende Achsiallast. 


Zur Beantwortung der Stabilitaétsfrage wird die Methode der kleinen Schwingungen heran- 
gezogen. Man iiberlagert den Komponenten der Grundbewegung (gekennzeichnet durch den 


e 

* Dissertation Clausthal 1950. Referent: Prof. Dr. E. Mettler, Korreferent: Prof. Dr. K. Jung. Der Ver- 
fasser ist Herrn Prof. Dr. FE. Mettler fiir die Anregung zu dieser Arbeit und fiir verstan 
zu herzlichem Dank verpflichtet. 

' E. Mettler, Mitt. Forsch.-Anst. Gutehoffnungshiitte-Konzerns 8 (1940) 5.1. 

2 K. Klotter, Forschg. Ing.-Wes. 12 (1941), S. 209. 

° S. Lubkin-J. J. Stoker, Quart. Appl. Math. 1 (1943), S. 215. 
4 FE. Mettler, Ing.-Arch. 16 (1947), S. 135. 
5 
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dnisvolle Férderung 


E. Mettler, Ing.-Arch. 17 (1949), S. 418. 
Vgl. E. Mettler, Ing.-Arch. 13 (1942), S. 97. 
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oberen Index Null) infinitesimale Zusatzverschiebungen (gekennzeichnet durch einen Quer- 
strich)?: 
u= ul +1 y, w= w0) + w, (3) 
Hier speziell ist der Grundbewegungszustand die reine Langsbewegung: w = 0. Geht man 
mit dem Ansatz (3) in (2) ein, so findet man fiir die Glieder nullten Grades in den u und w 


te ty 1/2 U2 
EF f[ 5 I 
| (A + V0 — T%) dt — | fe ul) dx + P(t) uw (=) -- io | u()* dx dt , (2)° 
ty ty Sp ip: 
und nach Variation ergibt sich fiir die Grundlingsbewegung als Eulersche Differentialgleichung 
— EF wu) + ul?) == () (4) 
mit den Randbedingungen 


~ 


(0) (»= = 


1)=0, (4’) 


bo} 


EFu(2=3,1)+ P(t) =0. (4””) 


Diese Differentialgleichung wird spiter zu lésen sein, um die Grundbewegung formelmafig zu 
beschreiben. 

Die Glieder ersten Grades in wu, w heben sich heraus, weil sie die ersten Variationen von (2)°sind. 
Die Glieder zweiten Grades hingegen liefern 


ty ty 1/2 1/2 1/2 


aye = LOTT ES us 
A i, — 1/2 =U/2 =I 
und nach Ausfiihrung der Variation hat man fiir die gesuchte Querbewegung w(x, t) als Eulersche 
Differentialgleichung 


— EF (u!) w,), + EJ We een pW, =0 (5) 
mit den Randbedingungen 


= L = l - 1 de l 
please =w.(2=—3:1) == 0"; w(x=z,1)=w,(2=3,1)=0 
\ 


und beliebigen Anfangsbedingungen gefunden. Da u° (x, t) zufolge (4) als bekannte Funktion 
anzusehen ist, hat man fiir die Querbewegung w (x, t) eine gemischte Rand- und Anfangswert- 
aufgabe fiir eine partielle Differentialgleichung mit zeit- und ortsabhangigen Koeffizienten. 

Die aus (2)? ebenfalls durch Variation herleitbare Eulersche Differentialgleichung fiir die 
Funktion wu beschreibt die freien Langsschwingungen des Stabes. Diese interessieren im weiteren 
nicht, da sie unabhangig von w sind und daher auf die Stabilitaétsuntersuchung keinen EinfluB 
haben. 

Da Gleichung (5) zwangslaufig aus dem Hamiltonschen Prinzip folgt, kann sie vollstandig 
mit Hilfe der allgemeinen Theorie? diskutiert werden. Man kann (5) jedoch auch nachtraglich 
als Gleichgewichtsbedingung am Stabelement fiir die Krafte in Querrichtung deuten. Der erste 
Summand ist bis auf den Faktor dx der Anteil der Querkomponente am Stabelement, welcher 
von der Langskraft herriihrt. Ebenso stellt der zweite Summand den Querkraftanteil und der 
dritte die Tragheitskraft am Stabelement dar. Man kann also sinnvoll in diese Gleichgewichts- 
bedingung eine geschwindigkeitsproportionale Dampfungskraft einfiigen, auch wenn diese nicht 
aus dem urspriinglichen Ansatz herzuleiten ist. 

Man darf zwar nicht erwarten, mit dem Dampfungsglied die in der Wirklichkeit vorliegenden 
Verhaltnisse vollstandig zu beschreiben, aber der Einflu®B einer Dampfung soll wenigstens in 
dieser Naherung in die Betrachtungen einbezogen werden. 

Die so erganzte Gleichung (5) lautet daher (unter Weglassung des Querstriches) 


EJ We2nx— EF (uO w,), + wor + Bo =0. (6) 


Um auch diese Differentialgleichung als eine Eulersche Gleichung auffassen zu kénnen, hat man 


1 Vgl. FufBnote 4 von S. 162. 
2 Vel. Fufinote 5 von 5S. 162. 
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i 
die Ableitung ungerader Ordnung fw, durch eine geeignete Transformation zu beseitigen. beta 
man zu dem Zweck | 


pea 
w(x, t) =e 24 -W(x,t), (7) 
so wird 
B? | 
ED Woes GB EUS We)o sc Wee oe (8))) 


mit VW = W,,—0 an den Stabenden und beliebigen Anfangsbedingungen. 
Die Funktion W(x, t) hat keine unmittelbare mechanische Bedeutung mehr. Man kann W) 
etwa auffassen als die Querverschiebung eines Ersatzstabes, bei welchem die geschwindigkeits- 
proportionale Dampfung durch eine ausschlagproportionale ersetzt ist. . 
Die Differentialgleichung (8) wird fiir die weitere rechnerische Behandlung als Eulersche}) 
Gleichung des Variations-Ausdruckes 
t, U2 : 
2L=—f f [Es We el AS We at Wet PE Wada de Extn (9)) 
t, —1/2 4 Mh 
aufgefaBt. Die hierin auftretende Funktion u() kennzeichnet die Grundbewegung und ist aus (4) 
zu bestimmen. Die in den Randbedingungen von (4) vorkommende pulsierende Achsiallast P(t)| 
soll im weiteren aus einem ruhenden Anteil P, und einem harmonisch schwingenden Bestandteil | 
mit der Amplitude P, und der Frequenz w bestehen: 


P(t) = P,+ P, cos wt. 
Dann ist die Grundbewegung nach (4) mit (4’) und (4) leicht anzugeben: 


l 2 l 
Py (x + = sin soa es) 
2 1% EF 2 
u\) = a cos wt. (10) 
EF EF (po wo | 
aor cos |/ FE 1 | 


| 


Um in der weiteren Rechnung von den speziellen Stababmessungen frei zu werden, sollen 
geeignete Veranderliche und Parameter eingefiihrt werden. AufSer einer dimensionslosen Orts- 
koordinate = x/l wird auch eine dimensionslose Zeitzahlung 


Sot (11) 


eingefithrt. Die Langsanregungsfrequenz w wird ebenfalls durch eine reine Zahl 6 in Einheiten 


f=VEJ|ul der Biegeschwingungsfrequenzen gemessen: 
he 8 gungsireq & 


ee are 88 ep | 

w= bf = 5/8. (129 

Der Kiirze wegen wird die Zahl 6 trotzdem als Frequenz bezeichnet. Ferner wird auch die Ampli- — 

tude des pulsierenden Lastanteiles P, in einem relativen Ma8 angegeben. Durch Vergleich von | 

P, mit der Eulerschen Knicklast Pr eines Stabes mit gleicher Lagerung (ein Ende frei, das 
andere eingespannt) werde die relative Lastamplitude 
= P, — By 

ans Ed (13) | 


eingefiihrt. Da bei einer normalen technischen Aufgabenstellung die pulsierende Achsiallast | 
weit unterhalb der Knicklast bleibt, ist ¢ im allgemeinen klein gegen 1. | 

Im weiteren ist es noch zweckmiabig, die Lingsanregungsfrequenz w mit der tiefsten Eigen- 
frequenz W 1, der Stablangsschwingungen zu vergleichen. Das so eingefiihrte Resonanzverhaltnis | 


/ 


OoL 


wird in der weiteren Theorie eine wesentliche Rolle spielen. Es sei schon hier darauf hingewiesen, 
dafs in den bisherigen Untersuchungen! stets 7 = 0 vorausgesetzt und damit der Resonanzfall 


1 Vgl. die FuBnoten 1, 2 und 3 von S, 162 sowie K. Kejwal, Diss. Hannover 1947. 
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74 = 1 vollstandig unterdriickt wurde. Legt man wieder die gleiche Lagerung fiir den Stab zu- 


grunde, so ist: 
au 1/EEF 
UR iy | ne (15) 


SchlieBlich soll die ruhende Vorlast P, durch den dimensionslosen Parameter 202 — Hot 
gelegt werden. Fiir den Stab ohne Vorlast (P,—= 0) ist also «= 0. Ey 

Geht man auf die neuen Koordinaten und Parameter iiber, so folgt zunachst aus (10) das 
im Variationsausdruck (9) bendtigte 


Tee el 
85° (a aye 
(0) FIG mas a 
foe EE EY on Seed me} 
No YI 
und damit 
if 
s, 1/2 cos 5 n £+7) 
pk 2 9 WV? 
EJ ae (20 +e:4n * coss a é 
s, —1/2 AUS GS} 


| (17) 


+ Wi.— & w:— 20? dé ds = Extr. 


Hieran kniipft die weitere Rechnung an. 

3, Ausgangsgleichungen, a) Ubergang zu einem System gewéhnlicher Differen- 
tialgleichungen. Zur Beurteilung des Stabilitatsverhaltens des ,,langsangeregten‘‘ Stabes 
mu man wissen, ob die Querverschiebung w(x, t) fiir unbegrenzte Zeit endlich bleibt oder nicht. 
Um zu einer formelmaBigen Naherungslésung fiir w(x, t) za kommen, aus welcher man dieses 
ersehen kann, geht man so vor!, daB man das Variationsproblem (17) durch einen Ansatz 


W (&, s) = es X; (€) T; (s) (18) 


genahert lést. In diesem Produktansatz sind die X; geeignete, fest gewahlte Ortsfunktionen und 
die Zeitfunktionen T;(s) sollen so bestimmt werden, da der Variationsausdruck einen statio- 
naren Wert annimmt. Man hofft, dabei mit einer méglichst kleinen Zah] N von Ansatzgliedern 
eine brauchbare Naherung zu erhalten. Auf diesem Mittelweg zwischen den direkten und in- 
direkten Methoden der Variationsrechnung gewinnt man eine Naherungslésung zunachst fiir 
W(é,s) und damit wegen (7) auch fiir w(x, t). 

Als Ortsfunktion im Ansatz (18) werden die Eigenfunktionen (Eigenschwingungsformen) des 
konstant gedriickten Stabes gewahlt, da diese ein bewahrtes vollstandiges Funktionensystem 
fiir Entwicklungen bilden. Aus der Differentialgleichung fiir die Eigenschwingungen 


pe 
Wadeye 20 Weare 0 Wi ays VO (19) 
mit den Randbedingungen W = W, = 0 ergeben sich durch den Produktansatz W = X(&) T (s) 
das Eigenwertproblem 


RP ogee = 36; X80 1 


mit ee x; ( 1) ; (20) 
cae ae 
und die Schwingungsgleichung 
627; 4 (* Tap) T= Zi, (21) 


Hierin bezeichnen wie immer im folgenden Striche Ableitungen nach der Ortskoordinate &, 
Punkte Ableitungen nach der ,,Zeit‘‘ s. Die Eigenwerte von (20) bestimmen nach (21) die Eigen- 
frequenzen der Biegeschwingungen 


2 2 B? 
wz = 07 f? — 4 * 


1 Vgl. FuBnote 5 von S. 162. 


| 
| 
| 
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Auf Grund der Selbstadjungiertheit von (20) sind die Eigenfunktionen X; orthogonal: 
1/2 
[ X:X,d—E=0 (i#m). 
= 4/2 


Sie werden im folgenden als durch die Bedingung 


1/2 
exe ded 
ye 
normiert angenommen. 
Multipliziert man (20) mit X,, und integriert partiell unter Beachtung der Randbedingungen 
iiber die Stablange, so folgt auf Grund der Orthogonalitat und Normierung 


1/2 1/2 al (mi), | 
f Xp Xi dé — 207 [ Xn Xia = 
— 1/2 


SYP) 


(22)) 
0; (==), | 
wovon spater Gebrauch gemacht werden wird. 

Der Ansatz (18) erfiillt bei der getroffenen Wahl der Ortsfunktionen die Randbedingungen) 
fiir alle Zeit und wird nun in das Variationsproblem (17) eingefiihrt. Nach kurzer Rechnung; 
wird dann 


;| > (2 ia r &2 Si + SF T na E (23) 
?—— —~___ | [? — 7— 6 cos S La dt ME Gia == [Bpaiir, 23) /1 
‘ A rf? — » 


fie deel i=1 m=1 
worin zur Abkiirzung 

; 1/2 

aU 1 , / 

a cos 7 n(é + 3) aE ORT (24) 


SOS tal 


S12 
gesetzt wurde. In den Koeffizienten von Tj kann man in bester Naherung den EinfluB einer | 
kleinen Dampfung auf die Eigenfrequenzen vernachlassigen. Von der Zulassigkeit dieser Nahe- 
rung iiberzeugt man sich leicht durch Einfiihrung des logarithmischen Dekrementes 

W,|__ Bs 

W, 2; ‘ 


eo; = In 


Mit w; = 0; f?— B?/4v? wird ein Vergleich der Eigenfrequenzquadrate mit und ohne Dampfung 
moglich: 


pi eae 


unged,. 


Wenn auch fiir #; keine genauen Zablenwerte vorliegen, so diirfte doch ein Dekrement von | 
9 = 0,1 bei technischen Aufgabenstellungen kaum erreicht werden. Dann kann man aber #?/z? 
gegen 1 mit einem Fehler von etwa 1°/ 9) vernachlassigen, so dah 


pe 2 

6} — — = 6? 

4 nef? | 

ist. Damit (17) einen stationéren Wert annimmt, sind die T; notwendig aus den Eulerschen 
Gleichungen 


N | 

6? 7; -- 0} T;— ecoss >| F,;T,, = 0 Ce ga (25) | 

m=1 | 

zu bestimmen und die Lésungstheorie dieser Gleichungen enthalt die Antwort auf die Stabili- | 

tatsfrage. Sie werden daher als die Stabilitatsgleichungen bezeichnet. 
Das System gewohnlicher Differentialgleichungen (25) mit periodischen Koeffizienten ist in 

gewisser Weise eine Verallgemeinerung der Mathieu-Gleichung. | 


N 
Die Funktion W(x, t) = 5) X;T; wird durch die Dampfung nicht beeinfluBt, da diese weder 
i=l 
in die Differentialgleichungen fiir X; noch in die fiir die T; eingeht. Zufolge der Transformation 
(7) ist w(x, t) = W(x, t), wenn keine Dimpfung vorhanden ist (8 = 0). Fiir den Fall nicht- 
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verschwindender Dampfung (8 4 0) tritt nach (7) vor W(x, t) ein mit der Zeit abklingender 
Exponentialfaktor, ohne da die weitlaufig zu berechnende Funktion W (<x, t) selbst sich dndert. 
b) Berechnung der Gleichungskonstanten fiir den Stab ohne Vorlast (P,=-0). 
Um die Stabilitatsgleichungen (25) im erforderlichen Umfange lésen zu kénnen, muB man zu- 
nachst die darin auftretenden Gleichungskoeffizienten 6; und F,,,; berechnen. Bei dieser Rech- 
nung wird es notwendig, den Stab ohne und mit ruhender Vorlast getrennt zu behandeln. 
Die Theorie des Stabes ohne Vorlast (P,— 0) erweist sich als die wichtigste, weil sie am 
weitesten allgemein durchfiihrbar ist und aus dieser Theorie auch noch gewisse Schliisse auf den 
Stab mit Vorlast gezogen werden kénnen. Nimmt man daher zunichst P, = 0, so vereinfacht 
sich wegen « = 0 das Eigenwertproblem (20) betrachtlich, und die symmetrischen oder anti- 
metrischen Eigenfunktionen mit ihren Kigenwerten sind leicht anzugeben!. Die Normierungs- 
faktoren sind aus der Bedingung 
1/2 
[ Xpdé=1 
—i/2 
zu bestimmen. 
Die Berechnung der Kopplungskoeffizienten F’,,; in (25) bedingt eine weitere Fallunterschei- 
dung, je nachdem man das Resonanzverhaltnis 7 als verschwindend klein annehmen kann 
oder nicht. 


x) Vereinfachte Theorie; Resonanzverhialtnis klein. Verlangt man, daB® die An- 


_regungsfrequenz w sehr weit unter der tiefsten Eigenfrequenz der Langsschwingungen bleibt, 
so ist 7 — W/Woz sehr klein, und man kann in angebbarer Naherung die in den Kopplungs- 
_koeffizienten (24) vorkommenden cos-Funktionen durch | ersetzen. Das ist zulassig, wenn in 


ee ee 1.) 
cos = (8 nice | 7 


= 14 au le A) at + Ont 


cos is 
cml 


Potenzen 7” vernachlassigbar klein sind. Mechanisch verlangt diese Forderung, da die Langs- 
kraft im Stab keinen értlich nennenswert verschiedenen Lingsspannungszustand anregen darf. 
Diese Voraussetzung ist hinreichend gut nur bei sehr schlanken Staben verwirklicht. Dies er- 
sieht man unmittelbar aus einer anderen Deutung des Resonanzverhaltnisses mit Hilfe des 


Schlankheitsgrades ): 


EI 
a) mie a @ 
Meret hie a EF Chey had 
EV ia 


Da die wichtigsten Instabilitatsbereiche in der Nahe der doppelten Eigenfrequenzen liegen, wie 
spater ausfithrlich begriindet wird, hat die vereinfachte Theorie nur eingeschrankte zahlen- 


maBige Bedeutung. Fiir die 2. Kigenfrequenz findet man z. B. 26,—= 133,35 und die Forderung 


ead 
fiihrt bereits zu so hohen Schlankheitsgraden, wie sie in der Technik nicht vorkommen diirften. 
Trotz dieser Einschrankung ist die vereinfachte Theorie sehr wichtig, denn sie erméglicht einen 
fast vollstandigen Uberblick iiber die Erscheinungen und wird daher zunachst ausfithrlich ent- 
wickelt. 
Fiir 7 = 0 hat man die Kopplungskoeffizienten 


EG ns = A Yi ws (26) 
einfach mit Hilfe des bestimmten Integrales 
1/2 
Ymi= | XmXidé 
1/2 


zu berechnen. Fiir mi folgt aus den differenzierten Differentialgleichungen (20) nach ge- 


1 Vel. z. B. Lord Rayleigh, The Theory of Sound, 2, Aufl. London 1926. 
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eigneter Multiplikation und Addition 
Xi, Xi — Xi Xn = (Om — 61) Xm Xi 


Nach Integration iiber die Stablange unter Beachtung der Randbedingungen erhalt man daraus 


das gesuchte 
WI YIALI2 
[X; Xn 1/2 
OF aa, 


22 [x Bek 1/2 
™ 


at 


Ymi = (m F t) Y (27) 


Da in den Randausdriicken Produkte einer geraden und einer ungeraden Ableitung auftreten, 


Ingenieur-Archiy 7) | 


wird nz—= 0, wenn nicht beide Zeiger m und i zu einer symmetrischen oder antimetrischen | 


Eigenfunktion gehéren. Das Gleichungssystem (25) zerfallt also (fiir 7 = 0) in zwei getrennte 
Systeme: je eines fiir die symmetrischen und fiir die antimetrischen Ansatzfunktionen. Diese 
Higenschaft der y,,; kann man auch unmittelbar aus der geometrischen Deutung des bestimmten 


Integrales fiir y,,; entnehmen. Man erkennt dann weiter, da® fiir 7 ~ 0 eine Belegungsfunktion — 


unter dem Integralzeichen auftritt und infolgedessen der Zerfall in zwei getrennte Systeme 

nicht mehr erwartet werden kann. — Die Berechnung der y,,; fiir m= 1 kann durch unmittel- 

bare Integration oder durch Grenziibergang auf einem von Lord Rayleigh’ angegebenen Wege 
ausgefiithrt werden. Es ergibt sich 

if 4112 

Yi = 77 ( [§ X; 


1/2 
—1/2 


— [XX 


Pht ie4 


Das Ergebnis der Zahlenrechnung ist in Tab. 1 wiedergegeben. Alle dort gegebenen Zahlenwerte 
sind unabhangig von den speziellen Stababmessungen und daher kann man mit ihnen die Sta- | 
bilitatsgleichungen (25) ein fiir allemal lésen. Dies ist ein grofBer Vorzug der vereinfachten | 


Theorie. 
Tabelle 1. « = 0. 
= i 1 2 3 4 5 6 7 8 9 10 11 
1] 12,3069 0 9,7306 0 —71,6152 0 6,1079| 0 | —5,0664| 0 4,3295 
2 46,0501 0 17,1412 0 —15,2025 0 = 0 = 0 
3 98,9082 0 24,3492 0 —22,9838] 0 20,8527) 0 | —18,7853 
4 171,5923 0 SU lal? 0 = 0 > 0 
5 263,9739 0 38,0301] 0 |—37,9606| 0 36,3966 
(Ui) a6 376,1481| 0 S50 eno 
a 508.0461} 0 51,2309) 0 | —52,2028 
8 Pni = Vin +! -.0 = 0 
9 860,8818] 0 | 232.2902 
10 = 0 
11 1269,7329 
5; = | 22,3729] 61,6728/120,9032|199,8596|298,5569| 416,9893| 555,1678| -| 713,0769| | 890,7300 


— nicht berechnet. 


Fiir eine spatere Verwendung sollen noch einfache Naherungsformeln fiir y,,; und Wii an- 
gegeben werden. Durch unmittelbares Ausrechnen findet man fiir die Werte der normierten 
Kigenfunktionen am Rande mit der Abkiirzung 6;= 1,2 die (noch streng geltenden) Formeln 


rT 
—sen cos— , 
47 1 “st { 9 3 2 
X; ss = + X; 5 = — 27+ — 
. Le 
—sgn sin >, 
r; T; 
" i ' —tgh > Sgn coss , 
xi"(+ = 9 xt"(—})=— 28 
—ctgh “ sgn sin = , 
2 2 


in denen sich die oberen Zeilen jeweils auf die symmetrischen, die unteren auf die antimetrischen 


1 Vgl. FuBnote 1 von S. 167 (Nr. 164—165). 


(28) | 
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Funktionen beziehen. Da fitr die Wurzeln r; die Beziehung gilt 


lim Peer n)=0 (i> 0 ganz) 


i> co 


und schon fiir i= 1 7;> > ist, so kann man mit einem maximalen Fehler von 2° 


hey 


Un t 


setzen. Dann ergeben sich aus (26) und (27) sehr iibersichtliche Naherungsformeln 


T; Tin 
—sen cos = cos— , 
O72 2p 2 r 
ee a 7) (27 *) 
ht We Ey ; 
m i : ow aa 
—sgn sin ~ sin 


2 2 


und 


foms 2 
Wi, TE — 27; 


wobei wegen r; > 2 y;; stets positiv ist. Sehr wichtig fiir spitere SchluBfolgerungen sind die 
Grenzwerte: 

lim y;; = 00, lim Wmi = 9. 

i+ oc t+ co, m= const 

B) Verscharfte Theorie; Resonanzverhaltnis nicht mehr klein. In einer Ver- 
scharfung der bisherigen Theorie soll das Resonanzverhaltnis 7 nicht mehr auf extrem kleine 
Werte beschrankt bleiben. Jedoch ist die Be- 
schrankung auf Werte 7 <1, d.h. auf unter- 7 
kritische Anregungsfre quenzen bei technischen 
Fragestellungen von selbst gegeben. 

Zur Berechnung der Kopplungskoeffizien- 
ten hat man das bestimmte Integral (24) nun 
ohne Vernachlassigung auszuwerten. Dies ist 
zwar elementar moglich, fiihrt jedoch auf eine 
verwickelte Abhangigkeit der y,,; von 7(w). 
Da zur Auffindung der Stabilitaétsgrenzen 
spater nach w entwickelt werden muB, kommt 
man auf diesem Wege zu uniibersichtlichen 
Ausdriicken. Es ist daher einfacher und kein ——— Porabelniherung 
nennenswerter Verlust an Genauigkeit, wenn —— cos #7 g+3) 
man in dem Integral (24) die Belegungsfunk- 
tion durch die Parabelnaherung G7 


cos 5.0 (€+3) ~ (cos 37 — 1) (e-+4) + 1-45 =e 


Abb. 2. Ersatz der Belegungsfunktion durch eine Parabelniherung. 
(7 <1) 


ersetzt. Wie man der Abb. 2 entnehmen kann, macht man fiir die nur in Betracht kommenden 
Werte 7 < 1 durch diesen Ersatz keinen wesentlichen Fehler. 

Mit dieser Parabelnaherung ergeben sich die F,,; aus (24) mit einer einfachen Abhangigkeit 
von 7 : 


45 NY 
7-47 


Ef 
A cos 5 7] —I1 
Fi a Wai A 0? a Ymi * (29) 
cos 577 cos 5 7 


Die yn; haben die schon in Abschn. «) berechneten Zahlenwerte und die y,,; sind noch durch 


Integration 
1/2 


Yni = Yim = f (e434) Xm Xj; dé (30) 


Sap 


| 
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zu ermitteln. Auch die y,,; enthalten die speziellen Stababmessungen nicht. Das Ergebnis der; 
Zahlenrechnung gibt Tab. 2 wieder. Da die y,,; im Gegensatz zu den y,,; fiir kein Zeigerpaar : 
i, m verschwinden, sind alle Stabilitatsgleichungen (wenn auch verschieden stark) miteinander) 
gekoppelt. Fiir die y,,; lassen sich entsprechende Naherungsformeln wie fiir die y,,; aufstellen. 


Tabelle 2370 0) | 
; 1 3 5 q 9 11 
a 
I 4,,2269 1,5129 | —2,5633 DPA = FIPS} 1,8789 | 
3 30.8014 1,3834 —7,4688 8,1282 eo Lo 
(=) — 83,7302 —2,9409 —11,4094 13,9102 
7 Ve ha 162,9921 —11,4089 —14,1231 
9 268,5638 2,0355 
ill 4004542 | 


Insbesondere findet man wieder einfache Grenzwerte, welche hier ohne Zwischenrechnung mit- | 
geteilt werden: 


lim 
i— oo m= const 


Vai 10s 


lim 4 
perce tut 
] 


Spater werden von F,,; die Entwicklungskoeffizienten der Taylorreihe fir eine bestimmte | 
Stelle 69 bendtigt: 


F,; = FO1 (5 


80) Fmt (3 — 89)? Fit ++ 
Diese sollen schon hier angegeben werden. Dafiir ist es zweckmabig, (29) unter Verwendung 


2FOne 
von 7 = —; in der Form 
p teeny sf 
cal eee 
te 


cos — 7 
18) 5) ] 
zu schreiben. Dann lassen sich die gesuchten Koeffizienten leicht angeben: 


(0) 2 Vmi mi } 
SE a5 — TH Yi} » 


1 eee 


(jy aa \ sin d9/A @ *) 3] 
Fini F- (Ymi— Ymi) Ty eos Sofa aS 2 (31) 
(2) __ 4 30? Are 1 + sin? 5o/A 
Fai = 2? (Ymi — Ym) 2! cos? dp/A ” 


Wie man sieht, geht jetzt in die Koeffizienten F,,; und damit in die Stabilitatsgleichungen der 
Schlankheitsgrad A als zusitzlicher Parameter ein. Aus der fiir endliches A geltenden Theorie 
gewinnt man durch den Grenziibergang 14 —> co die Aussagen der vereinfachten Theorie zuriick. 
c) Berechnung der Gleichungskonstanten fiir den Stab mit Vorlast (P)<0). | 
Da jetzt 1 0 ist, hat man das Eigenwertproblem (20) abhangig von dem jeweiligen Parameter- 


wert « zu lésen. Daher soll hier nur ein Beispiel fiir symmetrische (normierte) Eigenfunktionen | 


durchgerechnet werden. In Tab. 4 sind die Zahlenwerte 6; fiir den speziellen Parameterwert 
a —= 3,1119 angegeben. In diesem Beispiel ist ~« und damit im wesentlichen die ruhende Vorlast 
P, etwa gleich der halben Knicklast gewahlt: P) = 0,4906 Py. Innerhalb der Rechengenauig- 
keit ist dann 0,— 70,. Dies entspricht einem interessanten Sonderfall [vgl. Abschn. 5 a) q)]. 

Die Berechnung der Kopplungskoeffizienten macht wieder die Fallunterscheidung 7 = 0 
und 7+ 0 erforderlich. 

a) Vereinfachte Theorie; Resonanzverhdltnis klein. Die Berechnung der Kopp- 
lungskoeffizienten F’,,; fiir 7 — 0 fithrt auf dem gleichen Wege, wie in Abschn. 3b) x) be- 
schrieben, zu denselben Endformeln, da die Randausdriicke mit dem Faktor « verschwinden. 
Zur Bestimmung der y,,; hat man nur in (27) die aus (20) zu berechnenden Eigenfunktionen X; 
einzusetzen. Da diese von dem Parameter « abhingen, werden auch die y,,; abhangig von a. 
Weil die y,,; auBerdem die gleichen Eigenschaften wie in der vereinfachten Theorie des Stabes 
ohne Vorlast aufweisen [vgl. Abschn. 3 b) «)], entkoppeln sich auch hier die Stabilitatsgleichungen 
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in zwei getrennte Systeme, welche jeweils von den symmetrischen oder den antimetrischen An- 
satzfunktionen herriihren. Die y;; sind nur durch unmittelbare Integration zu gewinnen und 
nicht durch Randausdriicke in der von Lord Rayleigh angegebenen Weise auszudriicken. Die 
Zahlenwerte der Tab. 3 gelten fiir das erwihnte Beispiel x = 3,1119. 


Mabellewirx—— onl) 


at 1 2 3 4 5 7 
™m < . 
1 12,6061 0 10,9147 0 —7,6896 5,9053 
2 46,5130 0 27,8200 0 0 
Cp es 99,0280 0 25,3528 | —23,4513 
4 Wni = Pani 173,6330 0 0 
5 264,0461 38,8954 
il 509,2421 
d;= | 16,0997 | 53,9217 | 112,6973 | 191,3692 | 289,8685 | 546,2238 


B) Verscharfte Theorie; Resonanzverhdltnis nicht mehr klein. Die Ausfithrungen 
des Abschn. 3 b) f) fiir den Stab ohne Vorlast lassen sich fast wértlich auf den vorliegenden Fall 
P,+# 0 itbertragen. Die Berechnung der F,,; geschieht wieder mit der Parabelniherung und 
fithrt auch hier auf den Ausdruck (29). Die darin vorkommenden yp, ; sind bereits in Abschnitt 3 b) 
&) berechnet. Die noch zu berechnenden Zahlenwerte der y,,; erhalt man, wenn man die Eigen- 
_funktionen X; von (20) in (30) einfiihrt und die Integrale ausfiithrt. Die wiederum von « abhan- 
_gigen Ergebnisse sind fiir das Beispiel « = 3,1119 in Tab. 4 niedergelegt. Da die y; fiir kein 
Zeigerpaar m,i verschwinden, ist wie beim Stab ohne Vorlast in der verscharften Theorie eine 
Entkopplung der Stabilitatsgleichungen nicht mehr méglich. 
Die Zahlenwerte von y,,; und Ym; 


fiir das Beispiel « = 3,1119 (vgl. Tab. 3 Tabelle 4. « = 3,119. 

und 4) unterscheiden sich kaum von Sere ; | F : ; 

denen des nicht vorgedriickten Stabes ae | 

ee ed us ot ‘o eee i 1 | 4.1180 |—3,3612 | 1,9082 | —2,2905 
er Einflu8 einer erheblichen Vorlast (9 14.3226 | 9.8501 12761 

(im Beispiel fast die Halfte der Knick- 3 30,9535 | —22,3675 

last) gering. Daher ist die ausfiihrliche 4 Meas = Vara: 56,6489 


Behandlung des Stabes ohne Vorlast 
besonders interessant und man kann erwarten, daf} die dabei gewonnenen Ergebnisse in guter 
Naherung auch fiir den gedriickten Stab gelten. 


4. Beantwortung der Stabilitatsfrage mit der Stérungstheorie. a) Allgemeine Betrach- 
tungen und Formeln. Nach den Vorbereitungen der Abschn. 2 und 3 kann man nun an die 
Beantwortung der Stabilitatsfrage gehen. Dabei soll eine Bewegung w(x, t) quer zur Stabachse 
wie iiblich nur dann stabil heiBen, wenn w(x, t) fiir eine beliebig gewahlte Anfangsstérung 


(xt == ty) We (t,t == 2) 


auf unbegrenzte Zeit endlich bleibt. Falls dagegen |w| trotz beliebig kleiner Anfangsstérung 
mit unbegrenzt wachsender Zeit iiber jede Grenze wachst, soll die Bewegung instabil heiSen. 
Fir die Funktion w(x, t), welche die Stabilitaitsfrage beantwortet, wurde der Naherungssatz 


N 
wee ! SX, (x). T(t) (32) 
==] 


gemacht. Da die Ortsfunktionen X; in (32) als Eigenfunktionen von (20) gewahlt sind, ist die 
Erfillung der Randbedingungen fiir die Lésung w des Gesamtproblems fiir alle Zeit gesichert. 
} Man hofft, in (32) die Lésung w mit einer Summe aus miglichst wenigen Ansatzgliedern annahern 
zu kinnen, denn fiir einen N-gliedrigen Ansatz ergibt sichi. a. das System (25) von N Differential- 
gleichungen zweiter Ordnung mit periodischen Koeffizienten zur Bestimmung der Zeitfunktionen 
! T,. Die allgemeine Lésung von (25) setzt sich additiv aus dem Fundamentalsystem der 2 N par- 
jtikularen Integrale T,,; mit ebensovielen Integrationskonstanten zusammen. Somit hat also 
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(32) die Gestalt 
B | 


= BEG 


eG ee SKS T m (2), 


und es kommt im weiteren auf die Diskussion der Zeitfunktionen 


jal) (itt 1 ere) 


an. Da die T;,, als Lisungen von (25) gefunden werden miissen, ist das Studium dieses Diffe- | 
rentialgleichungssystems das mathematische Hauptproblem fir eine Beantwortung der Stabili- 
titsfrage. Bei E. Mettler! Abschn. II und III findet man die Lésungstheorie in dem Umfange, 
wie er hier benétigt wird, und diese Theorie muf fiir das weitere als bekannt vorausgesetzt | 
werden. Zunachst entnimmt man der dort gegebenen Spezialisierung des Floquetschen Theorems | 
auf das vorliegende System (25), daB das Fundamentalsystem die Form 


V-lem 


Tim=e “im (8) » er ee (33) 


tm 


haben mu. Hierin sind die Ym (s+ 27) = Ym (s) periodische Funktionen von s mit der Periode | 
2%, welche sich in eine Fourierreihe entwickeln lassen. Die N Zahlen @Q,, sind im allgemeinen | 
komplex und werden als charakteristische Exponenten bezeichnet. Wenn es unter ihnen einen | 
mit nicht verschwindendem Imaginarteil gibt, so spaltet der betreffende Teil der Lésung (32) | 


einen reellen Exponentialfaktor 


pee | | 
ibe my t (34) 

ab. Der Klammerfaktor in (34) werde als der ,,Stabilitatsindex‘‘ bezeichnet, denn er gibt die 
Entscheidung itber Stabilitat oder Instabilitat. Gibt es in (32) auch nur einen positiven Stabili- | 
tatsindex, so herrscht keine Stabilitat. Die wesentliche mathematische Aufgabe ist im folgenden | 


exp 


daher die Berechnung der charakteristischen Exponenten, abhangig von der ,,Frequenz‘‘ 6 und 
von der relativen Amplitude ¢ der harmonisch pulsierenden Langslast. 

Dabei wird noch die weitergehende Frage gestellt, welche zusammengehérigen Werte von 6 
und ¢ die stabilen Lésungen von den instabilen trennen. Dazu ist zunachst zu sagen, daB fiir hin- | 
reichend kleine ¢ die Lésungen im allgemeinen stabil sind. Dies Verhalten ist mechanich plau- 
sibel, denn fiir ¢ = 0 erhalt man als Grenzfall die stabilen Bewegungsformen des freischwin- 
genden Stabes. In dem anderen Grenzfall groBer ¢ wird man als Folge der groBen pulsierenden 
Langskraft erwarten, daB die Bewegungen des Stabes im allgemeinen instabil sind und stabile 
Lésungsformen die Ausnahmen bilden. Dies kann man auch rechnerisch leicht verfolgen, indem 
man zeigt, dafs (25) fiir ¢> im allgemeinen instabile Lésungen hat, was dann wegen (32) auch 
fiir die Querbewegung w(x, t) selbst gilt. 

Schreibt man namlich das Differentialgleichungssystem (25) in der Form 


N 
ee eim 7 
eT. oh ot Fe — ©0088! Fin T= 0 (ae een 
mit 
0 (iA~m) 
cm = -| 
en Ceti) 


so kann man in den Koeffizienten von T,, bei unbegrenzt wachsendem ¢ in erster Naherung 


den Summanden 0} ¢;m/F;, vernachlassigen. Fiir die so definierten Naherungslésungen TH 
gilt dann ; 


= N ‘ss 
&T;° — ecoss 3’ Fim Tn = 0 (i=1,...,N) (25’) | 


i=1 
und es ist leicht, zu weiteren Aussagen iiber die Lésungen fiir e—> co zu kommen, indem man 
aus dem System (25’) eine Folge von einzelnen Mathieu-Gleichungen ableitet. Da das asymp- 


Welle FuBnore a von 5) 162: 
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totische Verhalten der Lésungen der Mathieu-Gleichung fiir gro®e ¢ bekannt ist, verspricht dieser 


Weg Erfolg. Die gekoppelten Differentialgleichungen (25’) lassen sich vorteilhaft als Eulersche 
Gleichungen des Variationsproblems 


ON BV 2 N N wats 
| [ot ST? — eos Sie Pala ths = Wxtr. (35) 
sh tL i=1 m=l 
auffassen, und man kann hierzu stets eine orthogonale Transformation 
x N : 
Tj= Sa; Ti  (j=1,.-.,0) (36) 
finden, welche die quadratische Form 
N N eS 
Sy > Fei Te T; 
t=1 m=1 


° . . . nN 7 
in eine rein quadratische Form iiberfiihrt und auferdem S’ T;? ungeadndert 148t. Denkt 


i=] 
man sich diese Transformation ausgefiihrt, so erhalt man an Stelle von (35) 


soi Ni N = 
i, ae T;? — ecoss = y 1571 ds = Extri, (35’) 


worin auch die ¥; reell sind. An Stelle von (25’) treten nun die Eulerschen Gleichungen 
oT; —ecoss-y,T,— 0 Cee SS yy (25.5) 


. Diese Differentialgleichungen sind im Gegensatz zu (25’) eine Folge von entkoppelten Mathieu- 
- Gleichungen, von denen bekannt ist, da sie fiir e—> 00 auf im allgemeinen instabile Lésungen 
| fihren. Infolge der Transformation (36) gilt dies auch fiir die Lésungen des Systems (25’). 
| Damit ist der Nachweis erbracht, daB die Bewegungen des Stabes bei sehr groBen pulsierenden 
Langskraften im allgemeinen instabil sind. 

Die Beantwortung der Stabilitatsfrage fiir endlich groBes ¢ soll auf zweierlei Wegen ange- 
| strebt werden. Zunichst gibt es! eine umfassende Lésung des Problems der weitergehenden 
Frage nach denjenigen zusammengehérigen Werten von 6 und ¢, welche stabile und instabile 
» Lésungen trennen. Diese Theorie ist hinsichtlich der Allgemeinheit und der numerischen Brauch- 
» barkeit die weitreichendste. Sie wird hier nach ihrem mathematischen Aufbau als Stérungs- 
i theorie bezeichnet. Allerdings ist das Dampfungsproblem in diese Theorie nicht eingearbeitet. 
Ferner sind einige Sonderfalle hinsichtlich der Eigenfrequenzen in der genannten Theorie aus 
* rein rechnerischen Griinden ausgeschlossen worden. Die Stérungstheorie wird daher zunachst 
» zur Aufsuchung der Grenzkurven der Instabilitats-Bereiche in der (6, ¢)-Ebene fiir den Stab ohne 
| Dampfung angewandt. Im nachfolgenden Abschn. 5 soll dann noch ein zweiter Weg zur Beant- 
) wortung der Stabilitatsfrage beschritten werden. Dabei soll mit der Theorie unendlicher Deter- 
* minanten der Stabilitatsindex, d. h. im wesentlichen der charakteristische Exponent, unmittel- 
* bar berechnet werden. Dieser Weg kann nicht ohne griéferen Aufwand die gleiche Allgemeinheit 
4 erreichen, bietet aber einen gewissen Einblick auch in das Problem der Dampfung und gestattet, 
i die genannten Sonderfalle zu behandeln. Zudem bietet diese andersartige Behandlungsweise die 
| Méglichkeit, wesentliche Ergebnisse der Stérungstheorie zu bestatigen. 
| Die Berechnung der Grenzkurven zwischen Stabilitaéts- und Instabilitatsbereichen in der 
é (6, c)-Ebene mit Hilfe der Stérungstheorie geht von der bekannten Tatsache aus, daB die In- 
3) stabilitatsbereiche auf der 6-Achse in den Punkten 


zl 


One - (p >0, ganz) 
* bei Bereichen erster Art der Ordnung p und in 
em 
P 


4 bei Bereichen zweiter Art der Ordnung p beginnen. Die Grenzkurven dg (¢) werden durch 


' do= 6, +06 + 24+... (37) 
I 1 Vel. FuBnote 5 von S. 162. 
12* 
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dargestellt. Die Berechnung der Entwicklungskoeffizienten gelingt dadurch, dai man in den 
Lésungen (33) die periodischen Funktionen @,,; in Form von Fourierreihen annimmt und dann 
sowohl deren Koeffizienten als auch die charakteristischen Exponenten selbst in Form von) 
Potenzreihen nach ¢ entwickelt. Nach Einsetzen dieser Reihen in die Stabilitatsgleichungen (25) ) 
kann man in einem Rekursionsverfahren nach Art der Stérungstheorie die Koeffizienten so be-> 
stimmen, daB der Imaginarteil von @,, gerade den Grenzwert Null bekommt. Das Verfahren. 


aa (j)) 
kann hier nicht noch einmal auseinandergesetzt werden. Es wurde fiir die Berechnung der 0G 
in unverdnderter Form verwendet. Dabei wurde die Theorie insofern erganzt, als die Koeffi-: 


zienten 6) bis j= 2 und bis zur Ordnung p = 4 als geschlossene Formeln explizit in ihrer Ab: 
hangigkeit von den Gleichungskonstanten ausgerechnet wurden. Diese Formeln sind in Tab. 5) 


zusammengestellt, Auf die Wiedergabe des Ausdrucks fiir éf? fiir p= 1 mvB seiner Lange: 
wegen jedoch verzichtet werden. 


Tabelle 5. Grenzkurvenentwicklungen fiir Bereiche zweiter Art: dg = 09 + 6) et 6?) a + Oe. 


2 
P | of? og 
\ N (0)\2 
1 : 1 pee | ] | 1 | | 
Fo FO) 
4 Or Ok 
ieee EY) 1 


= | 
~ 2Vorde Sy 2 2 (8-1)? — 48, 


N (0) \2 
: 1 oes el 1 1 i 
cao Soe D \@,—1P— oe, | +1 — | 


n=T,k m=1 


ont (0) \2 
ieee lal 1 1 | 
i . 5 Bi, Sal 2 \G— 1p —& |G, + 1a 


n=r,k m=1 


Anmerkungen: 1. 6) = 


2, > : Summation nur tiber ry und k. 


m=r,k 


3. Grenzkurvenentwicklungen fiir Bereiche erster Art in den Formeln enthalten, wenn } 
r =k gesetzt wird. 


4. Bezeichnungen wie in Ing.-Arch. 17 (1949), S.418 mit Ausnahme von: 6, statt,, 
Cie INE SUBIE I). 


Man entnimmt dem Rechenverfahren und den Formeln der Tabelle 5, da®B sich fiir einen’ 
Bereich p-ter Ordnung die Entwicklungen (37) seiner beiden Grenzkurven erst in den Koeffi-; 
zienten 6?) unterscheiden. Demnach haben die Bereiche erster Ordnung schon zwei verschie-} 


dene 6), die Bereiche zweiter Ordnung zwei verschiedene 6), was in der Tabelle 5 durch die 
Doppelvorzeichen kenntlich ist. Die Bereiche dritter Ordnung unterscheiden sich erst in den} 
Koeffizienten 6) und so fort. 


Ferner deutet sich in der Formelzusammenstellung der Tabelle 5 eine weitere Tatsache an, 
welche in der allgemeinen Theorie eine groBe Rolle spielt und worauf auch in Abschnitt 5 noch 
besonders eingegangen werden wird. Die Entwicklungen dieser Stérungstheorie versagen nam- 


lich, wenn zwischen den Biegeeigenfrequenzen 5, und 6, und dem FuBpunkt eines Instabilitats- 
bereiches 6, eine Beziehung 


(0, + pd)? — oF = 0 (pS 0, ganz) 


besteht. Es mu8 also ausdriicklich verlangt werden, daB keine Eigenfrequenz ein ganzzahlig Viel- 
faches einer anderen ist. Dieser Sonderfall kann beim Stab ohne ruhende Vorlast nicht eintreten. 
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Beim Stab mit Vorlast kann dagegen bei geeigneter Wahl der ruhenden Last dieser Ausnahme- 

fall vorkommen und soll dann mit anderen Methoden in gewissem Umfange gelist werden. 
Kin weiteres Ergebnis aus der allgemeinen Theorie gestattet die Berechnung des Maximums 

der Imaginarteile der charakteristischen Exponenten. Die Formel soll fiir Bereiche erster Art 

erster Ordnung (6) = 2 6,) schon hier angegeben werden, da sie spiiter wichtig wird: 

| as 

8 2 E 4. Oe. (38) 


b) Stab ohne Vorlast. Resonanzverhaltnis klein (vereinfachte Theorie). Fir 
extrem schlanke Stabe (A—> oo) erhalt man durch den Wegfall der Resonanzbedingung sehr 


_vereinfachte Formeln. Zufolge (31) verschwinden fiir A—> co samtliche FY) fiir j > 1 und 
F®) ist durch (26) gegeben. Die numerische Rechnung ist daher vollstandig unabhangig von 
den speziellen Daten des Stabes und 1a8t sich in den durch (12) und (13) eingefithrten Parametern 
6 und ¢ ganz allgemein durchfihren. 

Da iiber die Konvergenz des Rechenverfahrens' hinsichtlich der fiir eine brauchbare Naherung 
notwendigen Zahl der Ansatzglieder N und der noch zu beriicksichtigenden Potenzen von « all- 
gemein keine Aussagen gemacht werden kénnen, ist zunachst das Konvergenzverhalten in einem 
Zahlenbeispiel erprobt worden. Hierfiir wurde der Instabilitatsbereich erster Art erster Ordnung 
gewahlt, fiir welchen also die Anregungsfrequenz in der Nahe der doppelten ersten EKigenfrequenz 
der Biegeschwingung liegt. Dabei wurden bis zu 6 Differentialgleichungen entsprechend einem 
| 6gliedrigen Ansatz aufgelést. Da das Gleichungssystem allein den symmetrischen Schwingungs- 
formen zugeordnet ist, hat man so die 11. Eigenschwingungsform beriicksichtigt. Wahrend die 
_Hinzunahme weiterer Gleichungen keinen grofen Aufwand erfordert, ist die Ausdehnung der 
- Rechnung auf héhere Potenzen von ¢ mit vergleichsweise viel gréBeren rechnerischen Schwierig- 
keiten verbunden, so daB die Rechnung nur bis zu Gliedern ¢® getrieben wurde. Das Ergebnis 
ist der Tabelle 6 zu entnehmen. Man erkennt, da man sich im Rahmen einer durch die verein- 


| Im 0,| max —— 


Tabelle 6. 


Bereich 1, Art Somes N +6())/2 64-8 +6?)/2 6,:€ +62)/2 o,:e| +0 é| 
Deo. 
1 ii —0,01472 «2 | +0,00090 «° 
Vereinfachte Theorie: 2 —0,02097 «2 | +0,00053 «? 
n=0(A=co) 3 | 6¢ = 44,7458- 1 | <0,24266¢) | —0,02149 «2 | £0,00046 <3 | | + 9 e 

Stab ohne Vorlast: 4 —0,02158 <2 | £0,00046 <* 

, P, == (a —— (0) 5 —0,02160 fx +0,00046 ee 
6 | —0,02161 <? | £0,00046 «® | 


Tabelle 7. Zu 6, = 22,3729 und 6, = 120.9032 (symmetrische Schwingungsformen) gehérend. 


5¢ 13. | 3p. |Se/5o| 58/5, | d/6, | 6&6, 


1 44,7458 | 1 | =0,24266} —0,02097 | +0,00053 


1 143,2761 | 1 | +0,02577) —0,00497 
3 241,8664| 1 | =0,06678) —0,00103 
—0,20424 o 
il 2253129 | 1 0 0.03633 0,00000 
2 Tl O3 Sisal 0 
—0,01989 L 
3 120,9032| J 0 +-0.00299 +0,00176 
—0,07156 | £0,07476 : 
ea ee nS ea 47,1587 | 1 0 +.0,06136 | $-0,00775 
3 80,6025} 1 0 —0,00387 | +0,00162 
~— 0 
A d ange eae a ONGhee 35,8190 | 1 0 —0,03015 0 
3 60,4516} 1 0 —0,00432 0 ; 
1. Art 2. Art 
9 
Grenzkurvenentwicklungen d¢ = 0) | 1 4 ete +——e8 +>" (= Os IN =2)- 
do do 09 | 


1 Vgl. FuBnote 5 von S. 162, IIL. Abschnitt. 
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fachte Theorie gerechtfertigten Naherungslisung durchweg auf zwei Ansatzglieder und auf in ¢ 
quadratische Glieder beschranken kann. In dieser Naherung wurden die Resonanzbereiche erster | 
und zweiter Art fiir die Ordnungszahlen p = 1 bis 4 fiir die ersten beiden symmetrischen und anti- 

metrischen Schwingungsformen berechnet. Vermutlich sind diese Formeln bis etwa ¢= 0,5 
brauchbar, wahrend in der Praxis wohl kaum Werte iiber 0,1 vorkommen diirften. Das zahlen- 
miaBige Ergebnis ist in den Tabellen 7 und 8 enthalten und die Grenzkurven sind in ‘Abb. 3 in| 


Tabelle 8. 
Zu do = 61,6728 und 5, = 199,8596 (antimetrische Schwingungsformen) gehérend. 
P oa Do 59/99 6/5, ot? [59 6) /5g bo 59/99 68) /5, df?) 15, 62) /6, 
2 123,3455 | 1 | =0,11949) —0,00659 
1 261,5323 | 1 | +=0,01165| —0,00147 
4 399,7191 | 1 | +0,04240) —0,00040 
S| Weiterey | ee on eee +-0,03765 
+0,00891 
2 130,7622 | 1 0 
—0,02818 
4, 199,8596 | 1 0 +0,00498 +0,01815 
2 41,1154} 1 0 —0,01614 | = 0,00888 
3 87,1766 | 1 0 —0,00833 | = 0,00015 
4 133,2404) 1 0 —0,00090 | - 0,00028 
2 30,8364 | 1 0 —0,01437 0 
A 65,3831 | 1 0 —0,00624 0 
4 99,9298 | I 0 —0,00167 0 
1, Art 2. Art 
a) 6?) 60) 
Grenzkurvenentwicklungen 6g = 0) | 1 + mC See tte ps eas (“== 0 Ni 2): 


3 3p 


0 


rt 


der (6,¢)-Ebene aufgetragen. Man entnimmt dieser Auftragung anschaulich, daB in der Technik 
wohl nur die Bereiche erster Ordnung Bedeutung haben werden, da die Bereiche héherer Ord- 
nung (p > 1) infolge ihrer geringen Frequenzbreite durch die nicht erfaBte Dampfung vermutlich 
ausgeléscht werden. Weiter entnimmt man der Anschauung, daf die Bereiche erster Ordnung 
beider Arten hinreichend genau durch das lineare Glied in der Grenzkurvenentwicklung 6¢(e) 
formelmaBig dargestellt werden kénnen. Ferner ist auffallig, da® offenbar die Frequenzbreite der 
Bereiche erster Art bei Ubergang zu Bereichen in der Nahe der zweifachen héheren Eigenfre- 
quenzen zunimmt. Im gleichen Sinne nehmen die Bereiche zweiter Art offenbar ab. Dies Ver- 
halten erscheint einer néheren Betrachtung wert. 
Die Beschrankung auf das lineare Glied ist fiir p == 1 und fiir ¢ < ausreichend: 


dg ~ Oot OMe. 
Dann kann man fiir d¢ eine sehr einfache Naherungsformel finden. Mit 


27 w.; 
iH) eee ie Lt 
= 


und den in Abschnitt 3b)qa) cingefiihrten Naherungen 


Uf ~ = (21+ 1) 


wird gendhert fiir Bereiche erster Art erster Ordnung der i-ten Eigenfrequenz 


le BGM eh 


i 
2 
Dit —— 2 ry, 
8% 


7 a> Il 


Damit ist auch eine Aussage itber den halben Offnungswinkel x; des Bereiches in dieser Naherung 


do = by 4 [Bat — 


gewonnen, denn der Faktor von ¢ ist gleich tg = (Abb.4). Demnach wachsen die Offnungswinkel 


der Bereiche erster Art und erster Ordnung mit héheren Eigenfrequenzen an, um dann aber 
schnell dem Grenzwert 2 7? zuzustreben, wie man auch anschaulich aus der Abb. 3 ersehen 
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kann. Fiir die Bereiche zweiter Art und erster Ordnung (6, + 6, = 69) ist es leicht, eine ent- 
sprechende Betrachtung durchzufiihren. Mit 


wird fiir e—>0 


2 er} 
4% -8rir. 


Fp lazes Agee |) ee ee 
Flee an ee greater craeurry 


Lor 7, 


(Fre) 2 Er) © 


Og = Oo “Eb 


Die Tangenswerte des halben Offnungswinkels der Bereiche zweiter Art nehmen also rasch auf 
den Grenzwert Null hin ab, wenn man zu hoheren Schwebungsfrequenzen iibergeht, was man 


auch der Abb. 3 entnehmen 
kann. Von denzuVielfachen 
der Schwebungsfrequenz 

| 69 = 6, + Ox 
gehérenden Bereichenhaben 
also nur diejenigen prak- 
tische Bedeutung, welche 
den tiefsten Eigenfrequen- 
zen zugeordnet sind, da 
ihre Frequenzbreite schnell 
abnimmt. 


dg (e) 


20; 


7 Abb, 4. 
Offnungswinkel der Bereiche erster 
Art und erster Ordnung, 


Aber auch die Bereiche 
erster Art erster Ordnung 
in der Nihe von 6, = 2 6; 
diirften nur fiir die tiefsten 
Frequenzen oder sogar nur 
fiir die erste Eigenfre quenz 
in der Technik wichtig sein, 
denn die Instabilitat wird 
in den Bereichen héherer 
Ordnung sehr viel schwa- 
cher. Das kann rechnerisch 
verfolgt werden, denn das 
MaB der Instabilitat ist 
durch den Imaginarteil des 
charakteristischen Expo- 
'nenten bestimmt. Mit (38) 
findet man fiir den Betrag 
des hier in Betracht kom- 
-menden charakteristischen 
Exponenten in _linearer 


Naherung fiir ¢e< 


Am? Wij 

8 0? 

qe 1; — 2r,; a 9 
beats in (2a-Elyt 


4 
(1 ana 62 (39) 
Erregerkraft 
»(e=konst.) wird also 
[Im @:| max bei Ubergang 
| auf hohere Kigenfrequenzen 
» schnell kleiner. Sieht man 
von dem Einflu8 einer etwa 
. vorhandenen Dampfung ab, 
so folgt aus (39), daB die Aufschaukelung der Querschwingungen in héheren Instabilitatsbereichen 
-langsamer vor sich geht als in tieferen Bereichen, da [Im g;| mit wachsendem 7 rasch abnimmt. 


— 


[Im Bleak o— 


gleiche 


400 


fe 
20, 
d= 
ob 
& 
ch 
50 700 750 


0, 
Q7 
0 200 250 300 550 


Abb. 3. Instabilitatsbereiche beider Arten bis zur vierten Ordnung 
fiir die ersteu vier Frequenzen. 
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Wenn man den DampfungseinfluB mit beriicksichtigt, so ist der Stabilitaétsindex (34) maf- 
gebend. Dieser wird fiir gréRere i kleiner und schlieBlich sogar Null, so daB damit die Dampfung | 
die Instabilitat verringert und schlieBlich ganz beseitigt. Dies Ergebnis entspricht genau den 
Eigenschaften gleichermafen belasteter Stabe bei freiaufgestittzter Lagerung’. 

In Zusammenfassung der Ergebnisse der vereinfachten Theorie kann man sagen, da tech- | 
nische Bedeutung wohl nur Bereiche erster Ordnung haben werden. Und zwar wahrscheinlich | 
nur diejenigen, welche zu den tiefsten Eigenfrequenzen gehéren. In guter Naherung lassen sich 
die Grenzkurven dg (e) dieser Bereiche durch in ¢ lineare Funktionen angeben. Die Bereiche (erster | 
Art) iiber den doppelten Eigenfrequenzen nehmen bei Ubergang zu héheren Eigenfrequenzen an 
Frequenzbreite schnell auf einen Grenzwert hin zu. | 
Von den Bereichen zweiter Art haben nur diejenigen | 
Bedeutung, welche zu den ersten Schwebungsfrequen- _ 
zen gehéren, denn ihre Frequenzbreite nimmt bei | 
Ubergang zu héheren Eigenfrequenzen schnell auf 
Null hin ab. 

Die Grenzen der Anwendbarkeit dieser vereinfach- | 
ten Theorie kénnen erst mit Hilfe der Ergebnisse der _ 
verscharften Theorie angegeben werden. 
Gg c) Stab ohne Vorlast. Resonanzverhaltnis | 
ii Oo WE a nicht mehr klein (verscharfte Theorie). In | 
09-20; T Res der verscharften Theorie mu®B man sich unter Ver- 
pease Genii ieh eae eae unterhalb der zicht auf die volle Allgemeinheit auf einen endlichen | 

Schlankheitsgrad / festlegen. Es wird daher ein Zahlen- — 
beispiel fiir einen sehr gedrungenen Stab mit A= 40 gewahlt. Fiir den wichtigsten Bereich 


A= | /n/inearer 
=4g{ Naherung 
yyy 


Vl 


7 


erster Art erster Ordnungist nach Tabelle 1 6) = 44,746 und wegen > p= ee ist dann in Bereichs- 
mitte das Resonanzverhaltnis 7 = 0,712. Der Resonanzfaktor wird 


fe 9980 


cos— 0 


2 


und ist bereits wesentlich von 1 verschieden. Die Berechnung der oly) hat nun unter Beriick- 


sichtigung der Frequenzabhangigkeit der Kopplungskoeffizienten mit Hilfe der FY) nach (31) 
zu geschehen. Es ist also bereits F°) gegeniiber der vereinfachten Theorie gedndert, wie man durch 
Vergleich von (31) mit (26) sieht. 

Mit einem zweigliedrigen Ansatz sind einige Resonanzbereiche berechnet worden. Besonders 
interessant ist die Grenzkurvenentwicklung fiir den Bereich in der Nahe der doppelten tiefsten 
Kigenfre quenz 


Og = 44,7458 (1 F 0,44799 ¢ + 0,29874 6? =F 0,71341 23 -L---). 


Diese Grenzkurven in der (0, ¢)-Ebene sind in Abb.5 dargestellt und die Ausweitung des Bereiches 
zur nahe benachbarten Resonanzstelle hin ist in die Augen fallend. Allerdings ist die Konvergenz- 
frage fiir den Fall einer derartigen Resonanznihe sehr kritisch zu beurteilen, wie man aus dem 
Anwachsen der Koeffizienten sieht. Man kann daher den Ergebnissen nur eingeschrankt quanti- 
tative Bedeutung beimessen. 


Ferner zeigen zwei weitere Zahlenbeispiele, namlich fiir den Bereich erster Art zweiter Ord- 


nung 
== (NYO) 206 
—_ Papaegss » il 2 “ cies — fee all 
96 : 9 + 0,04520° (> 2 
und den Bereich zweiter Art zweiter Ordnung: 
== (DOSS 6,6 
Oye == TGS fT ? 2 516.6 Ea reese 
7 ( + 0,00763° (> 2 iF 


daB auch fiir 7A 0 die Bereiche héherer Ordnung keine nennenswerte Frequenzbreite besitzen 
und daher auch kaum von praktischem Interesse sind. 


1 Vel. Fufnote 1 von 5S. 162. 
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Da man bei technischen Fragestellungen die pulsierende Lastamplitude hichstens auf etwa 
10% der Knicklast anwachsen lassen wird, so kann man sich auf ¢ < 0,1 beschranken!. Dann 
kénnen die Grenzkurven der praktisch wichtigsten Bereiche erster Art erster Ordnung wieder 
durch eine in ¢ lineare Funktion hinreichend genau angegeben werden: 


dc ~ bo + Of) « (= 01)e 


Hierin ist 6{) aus der Tabelle 5 zu entnehmen. Bei nicht zu hohen Genauigkeitsanspriichen darf 
man fiir y,;/y;; bereits fiir kleine i den Grenzwert 1/3 nehmen und erhalt damit 


6) xe 2 Vi eye ; 1 
: 6; \3cos do/A 3.) 


UL 
In dieser Naherung 1a4Bt sich die vereinfachte Theorie mit der verscharften vergleichen. Dabei 
bezeichne der angefiigte Stern * im weiteren die Ergebnisse der verschiarften Theorie. Dann ist 


1 2; 
aed Weaegae (j= 1), 
ec ul 


und fiir unterkritische Anregungsfrequenzen (7, < 1) und e—>0 ist stets: 
| 6" | =| 6e?| 


Die Abgrenzung beider Theorien gegeneinander ]a8t sich auch zahlenmaGig weiter verfolgen. 
Lat man zu, daB sich die Resultate der verscharften und der vereinfachten Theorie um hichstens 


on) 


ol) 


x%% von 6) unterscheiden, so ist damit fiir das Resonanzverhiltnis eine obere Schranke gegeben: 
= 1 5 1087, 
jeg 0.110 0,246 0,349 . 
Da 51 = 3 ist und bei Bereichen erster Art erster Ordnung auBerdem 6, = 2 6; festliegt, kann 
man angeben, bis zu welchen Werten des Schlankheitsgrades / es zulassig ist, nach der verein- 
fachten Theorie zu rechnen. Ist z. B. fiir 6 eine Abweichung bis zu 10% gegeniiber der stren- 


geren Theorie zulassig, so mu der Schlankheitsgrad oberhalb einer gewissen Grenze liegen: 


hare, ce ae 


LE 
30,349 


Da die Eigenfrequenzen des Stabes ohne Vorlast festliegen, sind die zulassigen Schlankheits- 
grade A; abhangig von der iten Eigenfrequenz angebbar. So mu fiir die Berechnung des 
-Resonanzbereiches erster Art erster Ordnung in der Nahe von 6, = 2 6; nach der vereinfachten 
Theorie 

Ae. 0L), Ag > 224,96 , Ay ASI,OL os 


sein, wenn 10% Abweichung zugelassen sind. 

Will man also auch nur die zweite Eigenfrequenz 6, in die Rechnung einbeziehen, so ist die 
vereinfachte Theorie schon auf extrem schlanke Stabe beschrankt. Um fiir technische Stabe die 
Bereiche auch zahlenmabig richtig zu erfassen, ist man gezwungen, die verscharfte Theorie heran- 
zuziehen. Da jedoch in der Technik ¢ wie erwahnt sehr klein sein wird und wohl auch nur die 
Bereiche Bedeutung haben, welche zu den tiefsten Kigenfrequenzen gehéren, kann man sich mit 
Naherungsformeln schnell einen Uberblick verschaffen. Die Naherungsgleichungen des Ab- 
schnitt b) lassen sich leicht auf den vorliegenden Fall hin erganzen. Man findet so fitr Bereiche 
erster Art erster Ordnung bei sehr kleinem € 


8a 2 ih 4 0) 1 
—~ + 2 =O) ee 
96 45 + [2a eallsecan ie RT aS, 
Fiir Bereiche zweiter Art kann man mit 


20? mi ymi 
6) == 55 (" Ze | Vm ) 


Vom 6; \ 8 Oo/A 


eine entsprechende Naherungsbeziehung aufstellen. Die mi, Ym; und 0; sind zwar unabhangig 


1 Vgl. die FufBnote 1 und 2 von S. 162. 
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von den speziellen Stababmessungen, sind aber nicht so einfach durch die Ordnungszahlen m, 1 |} 

auszudriicken, sondern ihre Zahlwerte miissen jeweils aus den Tabellen 1 und 2 entnommen 

werden. Beispielsweise gilt fiir den wohl wichtigsten Bereich zweiter Art in der Nahe der Schwe- 

bungsfrequenz 6, = 6, + 0, 

20° / 8,22 : 

: |. == 22,313; Og =— 1205903. 

5,6, a | SL ieee nits 07 5 3 . 

i) 


Diese Formel ist allerdings wie alle itbrigen der verscharften Theorie nur giiltig fiir 5 =< 


wla 


Das bedeutet wegen 6) = 143,27, daB 1 > 100 sein muB, womit hier immer noch verhaltnismaDBig | 
schlanke Stabe vorausgesetzt werden. Daraus kann iiberhaupt der bemerkenswerte Schluf gezogen | 
werden, daB Instabilitatsbereiche zweiter Art nur bei schlanken Staben in die Betrachtungen ein- 
geschlossen werden miissen, solange jedenfalls die Anregungsfrequenz unter der ersten Reso- | 
nanzfrequenz der Langsschwingungen bleibt. In Abschnitt 5 wird auSerdem noch an Zahlen- 
beispielen gezeigt, da der fiir die Starke der Instabilitat maBgebende charakteristische Expo-_ 
nent in Bereichen zweiter Art wesentlich kleiner ist als in Bereichen erster Art. Daher haben auch 
aus diesem Grunde die Bereiche zweiter Art wohl sehr viel geringere Bedeutung. 

In der verscharften Theorie entkoppeln sich die Stabilitatsgleichungen, welche von geraden | 
und ungeraden Ansatzfunktionen herrithren, nicht mehr. Als Folge davon treten Instabilitats- 
bereiche iiber den Schwebungsfrequenzen zwischen Eigenfrequenzen von geraden und ungeraden | 
Schwingungsformen auf. Hierfiir wird in Abschnitt 5 ein Zahlenbeispiel gegeben, woraus man _ 
ersehen kann, dafS Besonderheiten gegeniiber den bisher betrachteten Bereichen zweiter Art 
nicht auftreten. | 

d) Bemerkungen zu dem Stab mit Vorlast. Wenn eine ruhende Vorlast (P)+ 0) 
vorhanden ist, so 4ndern sich die Gleichungskoeffizienten von (25) mit dem Parameter a. Wie 
jedoch ein Vergleich der Tabelle 1 und 2 fiir « = 0 mit den Tabellen 3 und 4 fiir ~« = 3,1119 zeigt, 
ist der EinfluB von « auf y,,; und y,,; und damit auf F,,,; ziemlich gering. Dies ist um so be- 
merkenswerter, als zufolge 

Ast? Py 
ast. male Re 
die Vorlast P, mit fast der Halfte der Knicklast sehr hoch gewahlt ist. Lediglich die Kigen- 
frequenzen 0; des gedriickten Stabes sinken grundsatzlich mit wachsender Vorlast (wachsen- 
dem «), wie man auch aus Tabelle 1 und 4 ablesen kann. 

Beschrankt man sich fiir eine Naherungsbetrachtung wieder auf das in ¢ lineare Glied von d¢ 

und auBerdem auf Bereiche erster Ordnung, so andert sich mit wachsendem « der Koeffizient 


F°, 
OO = + a im wesentlichen nur mit dem Nenner Vm 6;. Daraus folgt z. B. fiir den 


ersten Bereich erster Art erster Ordnung nach Tabelle 2 und 3 


oy? | a=3,119 6, | a=0 22,4 
1) ib A 6, | «=3,119 16,1 


Allgemein geht aus dieser Betrachtung hervor, dab die Instabilitatsbereiche erster Ordnung eines — 
Stabes mit Vorlast gréfBer sind als die eines nicht gedriickten Stabes, was vermutlich auch fiir die 
Bereiche hiéherer Ordnung gilt. 


5, Beantwortung der Stabilitatsfrage mit der Determinantentheorie, Man kann die Stabilitats- 
frage aus den Gleichungen (25) heraus noch auf etwas anderem Wege als bisher beantworten und 
dabei auch die Dampfung in die Rechnung einbeziehen. Man mu dann den Stabilitatsindex (34) 
unmittelbar berechnen, um iiber Stabilitat oder Instabilitat entscheiden zu kénnen. Die wesent- 
liche Aufgabe ist also die Bestimmung der charakteristischen Exponenten @,, abhangig von den 
Parametern 6 und ¢. Im weiteren soll daher ein Rechenverfahren entwickelt werden, mit welchem 
man unmittelbar den Zahlenwert der charakteristischen Exponenten @,, (m= 1,..., N) fiir 
jeden Punkt der (6, e)-Ebene finden kann. Die Kleinheit des Parameters ¢ wird dabei nach wie 
vor vorausgesetzt. Im weiteren wird jedoch der in der Stérungstheorie hinsichtlich der Eigen- 
frequenzen ausgeschlossene Sonderfall 6, = p06, (p ganz) in die Rechnung aufgenommen. SchlieB- 
lich scheint es ferner lohnenswert, einige wesentliche Ergebnisse der allgemeinen Theorie auf 
diesem Wege zu bestatigen. 
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a) Ausgangsgleichungen der Determinantentheorie. Zur Berechnung der charakte- 
ristischen Exponenten von (25) soll die Theorie der unendlichen Determinanten von H. v. Koch} 
angewandt werden. H. v. Koch verwendet seine Theorie zur Integration gewéhnlicher Differential- 
gleichungen und behandelt als Beispiel eine Differentialgleichung, welche die Mathieu- Gleichung 
als Sonderfall enthalt. An dieses Beispiel wird in ahnlicher Weise angekniipft, wie es G. Lemaitre- 
4 ai E, Lamla? und J. 0. Fleckenstein* bei ihren Arbeiten iiber Hillsche Systeme getan 

aben. 


Dividiert man (25) durch 6? und setzt zur Abkitrzung 


+ 0; 
ee 
so lautet das Differentialgleichungssystem 
N 
Ti + @T:—ecoss >’ T, = (§=1,...,N). 


Si 


Macht man ferner fiir eine partikulare Lisung nach (33) den Ansatz 


f bere Ss Cis; Pa aia a) 


v=—o 


und fiihrt diesen in das Differentialgleichungssystem ein, so findet man fiir die Entwicklungs- 
koeffizienten c;,, das unendliche Gleichungssystem: 


N 
Fi i, é 1 
fe aah 262 q?—(v+ e)? (Cr,o-1 + 041) = 0 


v=) 


io otis CALNE aes) 


Die unendliche Determinante dieser homogenen Gleichungen ist eine Normaldeterminante. Sie 
mu8 fiir eine nichttriviale Lésung verschwinden und liefert so die Bestimmungsgleichung fiir N 
charakteristische Exponenten, wie sich zeigen wird. Die Elemente der Determinante werden 
zweckmaBig im AnschluB an G. Lemaitre-O. Godart? wie folgt bezeichnet: 


By.y+in-o+i =1 Nidamesete ae) 
& Fi i, 
9 g2 “1, | v—07 | 


(40) 


ieneee=o, ip ees | 


B - , >, —— 
N-v+i, N-v’ +i q? — (v-+e)? 


) 

; 

Zur Erlauterung dieser Bezeichnungen ist in Abb. 6 ein Ausschnitt der Determinante fiir N — 2 
schematisch wiedergegeben. Auber den Haupt- 
diagonalelementen vom Betrage 1 sind nur noch 

z Elemente vorhanden (durch x bezeichnet). 


a) Aufsuchung der analytischen Eigenschaften 
der Determinante. Wegen 


yy +i,No’ +i (O-+ 1)—= Bye_r+i,ne’—n +H (0)> 

und weil es gleichgiltig ist, welches Element man 

als Ursprungselement einer unendlichen Deter- 

minante ansieht, ist die Determinante eine peri- 

odische Funktion der Periode 1: 
A(e+ 1)= Ale). 

Die Funktion A(g) hat Pole an den Stellen 


aoe aD eeere ed Seal ea LaatN 2 


wie man aus den verschwindenden Nennern von 


(40) ersieht. Je nach der Vielfachheit dieser Pole 


ats Fh Abb, 6. Determinantenausschnitt zur Erlauterung der 
werden einige Falle getrennt untersucht. Bezsichnuipen. 


1 
| 1 H. v. Koch, Acta mathematica 16 (1892—93), S. 217, lehrbuchmafig dargestellt bei A. R. Forsyth, 
| Theory of Differential Equations, Part. III, Vol. 1V, 5. 348. Cambridge 1902. 

2 G. Lemaitre-O. Godart, Bull. Classe Sc. 5e Ser. Tome 24. Acad. Roy. Belgique Bruxelles 1938, 5S. 19. 
| 3 E. Lamla, J. reine angew. Math. (Crelle) 179 (1938), S. 134. 

4 J. 0. Fleckenstein, Comment. Math. Helv. 15 (1943), S. 367. 
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1. Einfache Pole (Residuen). Wenn alle g; voneinander verschieden sind und sich nicht um | 
ganze Zahlen unterscheiden, so sind alle Pole einfach. Fiir jedes beliebige Ziffernpaar r, k mu | 


dann 


G+ KF P p20 ganz 
sein. 
2. Zweifache Pole. Wenn 
ILqg-+qg=p oder Il.q¢.—4%=—pP (41) | 
fiir ein Zeigerpaar r,k ist, so fallen die Pole | 
I.+q-+n, mitden Polen + q+ ™ Guage) 
oder II. + q,-+ 7, mit den Polen + q,-+ ™% 
zusammen und es treten also jeweils zweifache Pole auf. Falls in I. r= k, so ist q,= p/2. Auch) 
hier treten Pole zweiter Ordnung auf, weil die Pole + q,-|- n, = F + n, und — q,-+ n= ‘ +n, | 
zusammenfallen. 
3. Vierfache Pole. Falls zugleich 
q=% ud ginu=P (ppp Pp!’ Z % ganz), (42) 


so ist auch q,— p’’/2. Mithin fallen die Pole g,+ p mit den Polen — q,-+ p’, + q+ p’ und || 
—q,-+ p’”’ zusammen. Es tritt daher ein vierfacher Pol auf. 

Es sind weitere Sonderfalle mit Polen noch héherer Ordnung 
denkbar. Diese sollen jedoch im weiteren ausgeschlossen werden. 
Dann kann man aus Abb. 7 fiir ein belicbiges Zeigerpaar r, k die 
Vielfachheit der Pole von A(g) entnehmen. Ist 6 und damit q,, q; 
d so gewahlt, da®B man sich nicht auf einer der Geraden des Dia- 

eramms befindet, so hat J(g) nur einfache Pole [vgl. 1.]. Auf 
einer der Geraden treten Pole zweiter Ordnung auf [vgl. 2.] und 
,, im Schnittpunkt zweier Geraden hat A(o) Pole vierter Ordnung 
[vgl. 3.]. 
B) Umformung der Determinante in eine rationale Funktion von 
Oe eee runs trigonometrischen Funktionen. 
Nach einem Satz der Analysis! kann man eine periodische 
Funktion der Periode 1 durch ctg-Funktionen und deren Ableitungen ausdriicken. 
1. Mit diesem Satz kann man in dem Fall einfacher Pole fiir die Determinante 


N 
A (0) — = [R;, x ctg x(o — 4) + Raactgu(o+ q:)|=K (43) 


Tk 


NIN 


3 
7 Zi 


NIN 


mit der noch zu bestimmenden Konstanten K ansetzen. Die Funktion der linken Seite hat keine | 
Pole mehr, wenn die Residuen R,,, Rj, entsprechend bestimmt sind. Die rechte Seite ist daher | 
nach dem Satz von Liouville eine Konstante. Da A (Q) eine gerade Funktion in 9 ist, so verlangt 
A(e) — A(—e) = 0 das identische Verschwinden jedes einzelnen Summanden in (43): 


[ete 2(@—— 9;). += ctg.a(o-- 4,)) Ry +> Rpg) == Ole (hel etn ie 


Setzt man also einfacher R;, = — R;,—= R;, so wird: 


N 
A(e)= K+ SAR, [etg x(o — gi) — ctg a(o + 9)]=0. 


t=1 


Zur Bestimmung der Konstanten K wird 9 = a+ if gesetzt und der Imaginirteil unbegrenzt 
vergréBert. Dann verschwinden in der Determinante alle Nichthauptdiagonalelemente und man 
erhalt: 
lim A(a + if) =1. (44) 
prc | 
Da andererseits lim ctg (~-+ i B) = — i, so ist damit K bestimmt: K — 1. Nach einigen trigono- 


B00 


* W. F. Osgood, Lehrbuch der Funktionentheorie, S. 466. Leipzig 1920. 
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metrischen Umformungen erhalt man endlich 


N : 
A(o) = 14 2a R; sin 2 7 q; (45) 


7 
ja C08 2 1g; — cos 27%@ 


_worin die Koeffizienten R; unmittelbar die Residuen an den Stellen 9 = q; sind. Fiir die Einzel- 
-gleichung ist N= 1 und in (45) ist das Ergebnis von H. v. Koch! enthalten. 


2. Die Einbeziehung der Sonderfille (41) ist fiir die Beantwortung der Stabilitiitsfrage not- 
wendig. Da A(e) nun Pole zweiter Ordnung hat, ist ein entsprechend erweiterter Ansatz er- 


forderlich N 
A(o) = K+ ek [ctg x(@ — q;) — ctg z(o + 4) ]+ 


ta 4 . 
= paca ldo ctg %(@— 4.) + Go te mer Mk \ + 


| Hierin ist bereits beriicksichtigt, daB A(g) eine gerade Funktion ist. Durch den Grenziibergang 


(44) ergibt sich auch hier K — 1, wegen lim sin-?(a + if) = 0. Nach dhnlichen Umformungen 


wie in (45) erhalt man Ja 
N 22R,sin2274q, 1 —cos 2 2q, cos 2 2 
Se esa oct 2p qk @ ; 
A(e) r= cos 2% q; — cos 27@ ve ¥ (cos 270g), — cos 270)? (46) 


=r 
3. Der Sonderfall (42) soll nur in der speziellen Form behandelt werden, wie er in einem 


spateren Zahlenbeispiel auftritt. Dort ist 


ea nO nt cin Gd, 


und man schlieBt ganz entsprechend wie in den vorigen Abschnitten auf 


cos 2 %q; —cos 2 7@ 


N 
ee d 1 , B 1 
A(e)=1- > a4 a Py a erg (0 *} t+ 7 P, jactex(o—4), 
it 


xr 

+k 
worin beriicksichtigt ist, daB 4 (0) eine gerade Funktion ist und die Pole zweiter und vierter Ord- 
nung an den Stellen + 1/2-+ n und — 1/2 +n liegen. Nach einigen Umformungen wird 


2at P| se 2 cos? x (o—-)| 


N 
2 aR; sin 2 74q; Geode 
A(e)=1+ > eae : 


47 
cos 2 %q; — cos 220 , 1 a, 1 (47) 
: sin? a Cia sin''2\@ — 5 


oar, k 
Offenbar sind die charakteristischen Exponenten notwendig aus den Gleichungen 4(0) = 0 zu 


-bestimmen. Nach (45), (46) und (47) sind dies Gleichungen N-ten Grades in cos 270, und man 
-erhalt mithin N Lésungspaare fiir die charakteristischen Exponenten 


+ Om> — Om CE ica SP aa AND 


Wenn ein oder mehrere 9,, einen nicht verschwindenden Imaginarteil haben, besteht nach (34) 
die Méglichkeit einer unbegrenzt wachsenden Partikularlésung. Ferner treten ebensolche 


Lésungstypen i. a. auch dann auf, wenn es mehrfache Wurzeln Q,, gibt. Dieser Fall kann im 
weiteren jedoch beiseite gelassen werden?. 

y) Bestimmung der Koeffizienten. Zur zahlenmafigen Auswertung von (45), (46) und (47) 
miissen noch die Koeffizienten berechnet werden. Ihre Berechnung ist die Hauptschwierigkeit, 
und man ist gezwungen, dabei zu einer Naherungsrechnung iiberzugehen. Zu dem Zweck werden 
die R;, P, und P”, nach Potenzen des kleinen Parameters ¢ entwickelt. Diese Entwicklung soll 
auf Potenzen ¢° einschlieBlich beschrankt werden. 

Nach H. v. Koch? gilt fiir die hier zu betrachtenden Normaldeterminanten ein ahnlicher Ent- 
wicklungssatz wie fiir endliche Determinanten: 


. 0 pee 1s | pr 
Bat oe SY Bee On Belt Bee. (48) 
Pq B oT, 


ye Pye 


mit p< gc rc ... ganz. 


1 Vgl. S. 292 der in Fufnote 1 von 5.181 genannten Arbeit. 
2 Vgl. Nr.6 der in FuGBnote 5 von S. 162 genannten Arbeit. 
8 Vgl. FuBnote 1 von S, 181. 
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Da die Liésungen von (25) unabhingig vom Vorzeichen von ¢ sein miissen, darf die Entwick- 
lung (48) keine ungeraden Potenzen von ¢ enthalten. Tatsichlich iiberzeugt man sich leicht, 
da® auch die Hauptminoren dritter Ordnung in (48) verschwinden. Bei Beschrankung auf Gig 
der ¢® einschlieBlich, hat man also nur Hauptminoren zweiter Ordnung zu betrachten. 

1. Falls keine der Beziehungen (41) und (42) eintritt, gilt die Gleichung (45) und die R; sind] 
unmittelbar die Residuen: 

R; = lim (9 — qi) A(@)- 
o> 45 
Zu ihrer Berechnung hat man also nur diejenigen zweireihigen Hauptminoren aufzusuchen, 
welche fiir 90—>q; unendlich werden. Diese sind alle in der aus der Gesamtreihe der Haupt- 
minoren herausgegriffenen Teilsumme 


NN 
er Ne au ae 1 ; 1 ; 
— (rx) 2k UT elt —@— 04 | eg — OF | 
t=1 j=1 a J ; d (49) 
1 1 
l ae ie : 
» Jae (e= Dla = (oe 25) laf ee 1 ge 2)? | j 
enthalten. Die Summanden, welche den Nenner (9 — q;) enthalten, liefern 
(Fey | Memeo 1 ; 
R=? 4 2 ! Oct, 50 
"4, (4q7—1) . PAF saw ' @— (a+ 1? me ( )| 


Ai 
worin das Glied j = i seiner iiberwiegenden GréBe wegen aus der Summe herausgezogen wurde. 
Fiir die Mathieu-Einzelgleichung (N= 1) fallt die Summe in (50) weg, und es bleibt nur dieses 
erste Glied stehen. Es findet sich schon bei H. v. Koch}. 

2. Die Berechnung der Koeffizienten fiir die Falle (41) ist fiir Pp >1 sehr weitlaufig. Daher 
werden im weiteren nur die praktisch wichtigsten Falle g,-+- q,—= 1 und q,— q,—1 durch- 
gerechnet. Wahrend auch in (46) die R; wieder die Residuen sind, findet man fiir die Pole 
zweiter Ordnung aus (46) 

lim (0 ~ q)?A(o) = — Py. 

eg, 
In (49) liefert dann fiir Ge g,— 1 der erste Summand (2 N— 2) Pole erster Ordnung (fiir 
i= koderj=r) und einen Pol zweiter Ordnung (fiir i= k und j =r), der zweite Summand 
N Pole erster Ordnung (fiir 1 = k), der dritte Summand N Pole erster Ordnung (fiir i= r) und 
der vierte Summand, wie alle weiteren nicht hingeschriebenen Summanden der zweireihigen 
Hauptminoren, ce keine Pole. Also gibt es (4. N— 2) Summanden in? 


hep oes ctls a) 
2 2 4, (4; ently i ! 2; 2 4, ee ese ! 


i=] 
k 


L 


| (51) 
N we ‘ N jue : 
my ore ! bs (nF) Lei One | 
et 2a lG—G@ FT" <y 2a f@G—@— 2 | } 
und einen Summanden 
2 2 
i Sa NN aoa 
OS a) 4g, 4, Deadhne (2) 
Wenn dagegen q, — q,= 1 ist, so findet man entsprechend 
R, = const. e? + Oet (51’) 
und 
(aap Fe) 
P,=4 —-+ Oe. (52’) 


Aq, 4, 


1 Vgl. FuBnote 1 von S. 181. 
 Wegen (1 — q,)? —(@ — 1)? = —(e —q,) (0 + a, — 2). 
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Wie sich spater zeigen wird, ist es bedeutungsvoll, daB (52’) sich im Vorzeichen von (52) unter- 
scheidet. Die genaue Formel von (51’) ist fiir das Weitere uninteressant und daher nicht an- 
gegeben. 


3. In einem spater zu behandelnden Zahlenbeispiel fiir den Sonderfall (42) mu8 nur noch 
der Koeffizient P, fiir den Pol zweiter Ordnung bekannt sein. Den Pol vierter Ordnung kann 
man auBer Ansatz lassen, da seine e-Potenzentwicklung von P; mit ¢4 beginnt und hier nur 
Potenzen bis ¢° beriicksichtigt werden. Den Koeffizienten P, findet man aus (49) durch Betrach- 
tung des Summanden 


€ Fx) ume i ee 
2 52° **) Tq? — 0%) [a2 —(@— 1") ’ 
aus welchem man den Nenner (@ — q)° abspalten kann. Mit q, = 1/2 ist dann 
2 
P,=—(s85 Fa) + 004. (52’’) 


Fiir die Residuen erhalt man ebenso wieder R; = const. ¢2 + 0e4, worin die Konstante fiir das 
Zahlenbeispiel jedoch nicht weiter benétigt wird. 


b) Einige allgemeine Schliisse. 

Mit den Ergebnissen des Abschn. a) lassen sich einige Resultate der allgemeinen Theorie 
mit Hilfe der Determinantentheorie bestatigen. 

Fiihrt man die Abkiirzungen 


z(0) = cos 270, 2; == cos 27q;, m;(e) == 27 R; sin 2704; (53) 


ein, so wird im allgemeinen Falle, wo keine der Beziehungen (41) oder (42) eintritt, die Deter- 
minantengleichung (45): 


N 
m; 
A aise Dae 
oder 
N N N* 
IT (%,— 2) + 3S) m, II (z— 2) = 0. (54) 
i= i=l Ee 


Damit der Fall der Mathieu-Gleichung in (54) enthalten ist, werde festgesetzt, daB fiir N— 1 
das letzte Produkt gleich 1 ist, was durch den Stern am Produktzeichen kenntlich sein soll. 
Far N= 1 ist dann also z,— z-++ m, = 0 und eine komplexe Wurzel 2) kann hierin nicht auf- 
treten, denn z, und m, sind nur reeller Werte fahig. Die durch komplexes @ bedingten Instabili- 
‘taten sind wegen (53) bei der Mathieu-Gleichung also nur méglich, wenn eine Wurzel |z,)| >1 
eintritt. Derartige Instabilitéten heiBen solche erster Art. 

Ist allgemeiner in (54) N >1, so kénnen die Wurzeln 2, den Betrag | fiir e— 0 nur in 
der Nahe ausgezeichneter Frequenzwerte erreichen oder iibersteigen. Da nimlich fiir «= 0 
die m; verschwinden, hat man 


[|= 1 fiir 20g, p’u (p’20, ganz) 
oder mithin 


é. 
g(e= N= f= F (pise0 ganz). 
0 


Die Méglichkeit solcher Instabilitaten erster Art besteht also fiir nur in der Nahe der Frequenzen 
6) =26,/p. DaB diese Instabilitaéten tatsachlich eintreten, soll wenigstens fiir den wichtigsten 
Fall (p= 1) 6, = 20; nachgewiesen werden. Es sei also g,— 1/2 und keine der Beziehungen (41) 
und (42) erfiillt. Dann findet man mit den gleichen Uberlegungen wie in Abschn. 5 a) f) 


ay ALI zt P, ieee It een 
AG ia a) Fag ( 2} ee) 


In (55) ist A(@) bereits als gerade Funktion angesetzt und der letzte Summand liefert die Pole 


sowohl an den Stellen 0 = + + nm und, 0= —+ +n. Nach kurzer Zwischenrechnung wird 


186 Weidenhammer: Der eingespannte, achsial pulsierend belastete Stab usw. _Ingenieur-Archiy 


aus (55) mit den Abkiirzungen (53) 
Al) = 14 yr (56) 


2, — & z+ Zp, 


i—1 
=k 


worin z= 1 ist. Aus (56) folgt durch pean ete der Nenner 


N N 
F= (4+ 2) Hf (%— 2) + (& + ae m; Ut —2)— 20 PHT (2) = 0. (56’)) 
x uh 4: Fk 


= 


Fiir ¢— 0 verschwinden die m,; und P, und es gibt eine Wurzel 24) = = — 2% = — 1 vom Betrage} 
eins. Diese Wurzel hat fiir ¢ 0 einen Betrag >1, wie mit Hilfe einer Taylorreihenentwicklung ; 
fiir z) (e) durch implizite Differentiation von (56’) gezeigt werden wird. Zunachst findet man, , 
daB in der gesuchten Entwicklung 


Zu (e) = 2) + en) + otf) +... (57)) 


kein lineares Glied enthalten ist, denn in! 


dz 
Fe F,7 = 90 (58)) 


ist (F.),-»= 0, (F)enp% 0 und mithin (Z) ete 0. 


e=0 
Demnach hat man in 


FP 2 FR + EF, gay =* (59) 


Wd? : ; : ; = | 
nur F’., und F, zu berechnen, um 2(?) = (zea) zufinden. Es ergibt sich mit P,—= P,, e?+ 0¢*: | 
e=0 


2 \de 
pedipes) 
(Ue wey = — 4a? P, IT (%;+ %) 
Z=—Z), a 
N 
(Feng =1+ Hf (+4) 
2=—Zz se 
und daraus z(2) = 27? P, 
Da nun (vgl. 52) P, = — (38) + Oe? 
ist, hat man (in den urspriinglichen Bezeichnungen) 
es Fie’ | 
Cosa Et 2a (3) eit Ost. (60). 


Da die rechte Seite den Betrag 1 iibersteigt, ist der Nachweis der Instabilitat fiir Frequenzen in | 
der Nihe der doppelten Biegeeigenfrequenzen durchgefiihrt. | 
Wenn die Anregungsfrequenz 6, ein ganzzahliger Bruchteil der Schwebungsfrequenzen 


Or + Ok Or — Ok 
1p 


(= (41), oder’ Yog— 


(41) 


wird, so besteht die Méglichkeit von Instabilitaten. Nach (46) ist nimlich mit den Abkurgungag ' 
(53), wenn man die Nenner beseitigt, 


4 N N N 

P(e, 3 oT 4) 31 (% —- 2)? Sm, IT (z-—2) — 
i= t=1 j=1 
ar Fr “a 


(61) | 
— 4a? Pi (1— m2) HT (3) = 0. | 


= 


oe 


* Tiefgestellte Indizes ¢ und z bezeichnen partielle Ableitungen, 
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Fiir ¢ = 0 gibt es die Doppelwurzel 2) = %. Fir e 0 besteht hier die Méglichkeit eines kon- 
jugiert komplexen Wurzelpaares und damit ist Instabilitat moéglich. Aus der Stérungstheorie! 
ist bekannt, da dieser Fall der Instabilitat tatsachlich nur in der Nihe der Frequenzen (41’) 
eintritt. Diese merkwiirdige Eigenschaft der Schwebungsfrequenzen soll nochmals mit der 
Determinantentheorie fiir den speziellen Fall p = 1 mit Hilfe der Entwicklung (57) nachgewiesen 
werden. Da in (61) (F,),.)—=0 und auch (F.).e=9 = 9 ist, kann man 21) nicht aus (58) be- 
rechnen, sondern mu dazu (59) heranziehen. Da auch (F,.)e=9 = 0 wird, findet man 


2(1) = |. a G ‘) . 


~ worin 
biy N 
(Feo = — 4 27 P, 2(1— 2) I (2; — 2) 
i=1 
und ca 
N 
( peo a ett (2; — %,) 
a; 


ist. Damit ist 
cos 27 0,,, = cos 27q, + 2m sin 27 q, /P, e+ 0e2, 
wenn man das Ergebnis wieder in den ursprimglichen Gréfen ausdriickt. Da nach (52) 
fir g,+ q,=1 P, <a 
/und nach (52’) 
| fir ¢.—q,= 1 PS 0 


ist, hat also cos 27o,,,, fiir e—>0 nur iiber den Schwebungsfrequenzen 6) = 6, +- 6, einen nicht 
_ verschwindenden Imaginarteil und fithrt so zu den bekannten Instabilitaten zweiter Art. Die 
| Nichtgleichwertigkeit der Schwebungsfrequenzen hat sich also auch auf diesem Wege bestiatigt. 
In dem Sonderfall q,—= 1/2, q¢,—= 7/2 1aBt sich durch eine ¢-Potenzreihenentwicklung auf glei- 
chem Wege zeigen, daB gegeniiber (60) keine Besonderheiten eintreten. 


c) Zahlenmafige Auswertung in einigen Beispielen. 


| Mit der Determinantentheorie lassen sich einige interessante Stabilitatsbereiche berechnen, 
) welche mit der Stérungstheorie nicht ohne weiteres aufzufinden sind. Insbesondere wird auch 
} die Dampfung mit in die Rechnung einbezogen. 

Durch eine geeignet gewahlte Vorlast von etwa der Halfte der Eulerschen Knicklast hatten 
‘sich die Zahlenwerte der Tab. 3 und 4 ergeben und innerhalb der Rechengenauigkeit ist der in 
| der Stérungstheorie ausgeschlossene Sonderfall 6,—= 76, verwirklicht. Das Stabilitatsverhalten 
‘eines Stabes mit diesen Daten und mit Dampfung wird in den nachfolgenden Zahlenbeispielen 
) aus den Stabilitatsindizes (34) abgelesen. Zu dem Zweck hat man zunachst die charakteristischen 
Exponenten @,, als im allgemeinen komplexe Wurzeln von Gleichungen N-ten Grades fiir 
}cos 270 zu berechnen. Da man aufberdem fiir jedes Wertepaar 6, ¢ eine gesonderte Rechnung 
| durchzufiihren hat, ist es wesentlich, daB man sich im Ansatz (32) in guter Naherung auf zwei 
Glieder (N= 2) beschranken kann. Dann sind die Gleichungen A(0)= 0 nach (45), (46) und 
(47) nur noch quadratische Gleichungen fiir cos 279. Zur Rechtfertigung dieser Beschrankung 
‘soll die Konvergenz hinsichtlich der Zahl der Gleichungen wenigstens in einem speziellen Falle 
‘numerisch verfolgt werden. 

Fiir eine Anregungsfrequenz (6 = 30,000) dicht unterhalb der doppelten ersten Kigenfrequenz 
‘und fiir einen Betrag der pulsierenden Lastamplitude von der Hialfte der ruhenden Vorlast 


fe 1d URS ae = P;) ergeben sich die charakteristischen Exponenten mit den Zahlenwerten 


‘der Tab. 9. Da es sich in diesem Beispiel um ein Wertepaar 6, ¢ in einem Bereich erster Art 
-handelt, liefert sogar ein eingliedriger Ansatz schon einen brauchbaren Wert. Um jedoch auch 
, die Bereiche zweiter Art zu erfassen, wird stets mit einem zweigliedrigen Ansatz weitergerechnet. 


1 Vgl. Fufinote 5 von S. 162. 


; : al 
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Tabelle 9. 
6 = 30.0000 6 = 0.2453 Glieder- 1. Wurzel 2. Wurzl | 3. Wurzel | 4. Wurzel 
; fap Z, + 2m Qy +27Q, +273 + 2% Q4 
—: Naherung liefert Sn] Se aa 
keinen Wert 1 —1,05388 tot eset as ; 
2 —1,05518 | +0,50000 1,53633 ae = 
-+1i+ 0.33070 
3 —1,05567 | +0,50000 1,53596 | 2,12547 = 
+i+ 0,33215 
4 —1,05565 , +0,50000 1,53596 | 2,12547 | 1,30269 
+i+ 0,33209 


Da in der naherungsweisen Berechnung der Gleichungskoeffizienten nur Potenzen é* einschlieB- 
lich beriicksichtigt wurden, und die Konvergenz dieser Reihen nicht gepriift ist, beschranken| 
sich die nachfolgenden Zahlenrechnungen auf ¢< 0,25. | 


Das Ziel der Zahlenrechnung ist jeweils, die Imaginarteile der @,, abhangig von den Para-, 
metern 6 und ¢ zu finden, um dann mit Hilfe von (34) die Stabilitatsfrage zu entscheiden. Fiir> 
die Stabilitatsgrenze folgt insbesondere aus (34) 

= = 1,9) 5, aN 62), 

[Fm Om (es 9)| = x 5F Gen 9 onal e0 9 ihe ( ) 

Die Daimpfungskonstante f soll durch das mit der tiefsten Biegeeigenfrequenz w,—= 6,f ge-+ 
messene Dekrement ausgedriickt werden: 


aliba 
Ss ERT 
Dann lautet die Grenzbedingung (62) 
[Lm Qm (€6, 9)| = ou 2 (62) 


und die Grenzkurven der Instabilitatsbereiche d¢ () sind aus einer Kurvenschar | Im @,,| ab- 
hangig von 6 mit dem Parameter ¢ leicht numerisch anzugeben. Die rechte Seite von (62’) er- 


ame, 
46 é 
G20 
G3 410 
£0 
ar Fes g 
JO aS 
es g 20; 
28 30 Se BY 36 ersterArt  ersterOrdnung 
d,=20; A=00  =G015 
Abb. 8. |Imo,| abhingig von 6 mit dem Parameter ¢ Abb. 9. Bereich erster Art erster 
in der Nahe von 6, =20,. Resonanzverhiltnis 7 = 0. Ordnung in der vereinfachten Theorie. 


gibt in dieser Auftragung eine gleichseitige Hyperbel, welche fiir verschwindende Dampfung 
zur Abszissenachse ausartet (vgl. Abb. 8, 10, 12). 


Im einzelnen werden die nachstehenden Instabilitatsbereiche durchgerechnet. 


«) Bereich erster Art erster Ordnung; Vereinfachte Theorie 7 = 0 (sehr schlanke Stéibe). Ini 
Abb. 8 ist |Im 0,| und die Dampfungshyperbel fiir diesen Bereich in der Nahe der doppelte 
ersten Biegeeigenfrequenz aufgetragen. Etwa in Bereichsmitte tritt der Sonderfall (42) ein 
Hier ist (45) fiir die Zahlenrechnung unbrauchbar, und man mu dann auf (47) itbergehen 
Dabei liegt in der Berechnung benachbarter Werte von Om nach verschiedenen Determinanten 
gleichungen eine willkommene Méglichkeit, die numerische Rechnung zu proben. Wie man sieht 
liegt dasMaximum von Im Q,, ziemlich genau in der Mitte des Bereiches, so daB die vereinfacht 
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Berechnung (38) der charakteristischen Exponenten? i. a. eine gute Naherung gibt. Die Grenz- 
kurven sind in Abb.9 eingetragen, wobei die Dimpfung mit einem Delaoment von 3} = 0,075 
sehr groB erwahlt ist. Zur Anregung der Instabilitat ist bei Vorhandensein dieser Dampione 
mindestens ein gewisser Schwellwert ¢, der pulsierenden Amplitude erforderlich. 


B) Bereich zweiter Art erster Ordnung; Vereinfachte Theorie 7 = 0 (sehr schlanke Stabe). Die 
doppelte Schwebungsfrequenz 6, = 6, + 6g ist in der vereinfachten Theorie die tiefste Schwe- 


tm er ne,| 


a6 : 
€=G25 
420 
Q3 Q10 
a7 gb-- 
ke ln, s died 
128 129 730 zwelter Art erster Ordnung 
d=, +d, A=-~ — B=GOT5 
Abb. 10. | Im 0, abhangig von 6 mit dem Parameter ¢ Abb. 11. Bereich zweiter Art 
in der Nahe von 6,=—6,+ 63. erster Ordnung in der verein- 
Resonazverhiltnis 7 = 0. fachten Theorie. 


bungsfrequenz, in deren Nahe Instabilitat der Querbewegung méglich ist, denn nach Abschnitt 
3c) a) zerfallen die Stabilitatsgleichungen (25) in der vereinfachten Theorie in zwei getrennte 
Systeme, welche von den geraden Eigenschwingungsformen (Index 1, 3, 5,...) und den un- 
geraden (Index 2, 4, 6,...) herrithren. In Abb.10 ist wiederum | Im @,,| aufgezeichnet, und 
daraus wurden die Bereichsgrenzen gefunden und in Abb. 1] eingetragen. 

Ein Vergleich der beiden Bereiche erster Ordnung ist aufschluBreich (Abb. 9, 11). Man er- 
sieht, dai die Bereiche zweiter Art eine sehr viel geringere 
praktische Bedeutung haben werden, als die Bereiche erster € Y 
Art, denn ihre Frequenzbreite ist sehr viel geringer. Auch Uj 
die Dampfung hat auf die Bereiche zweiter Art einen sehr G20 Y, 


19°? \/me-/me,| 


40 €=G25 
G10 
€9 
OS 
oS 
68 70 Wh 
5; +0; 
a7 EN Jb zwelterArt — ersterOrdnung 
; 69 70 1 A=30 B=0075 
Oy + Oz Abb. 13. Bereich zweiter Art erster 
Abb. 12. | Imo,'| abhangig von 6 mit dem Parameter ¢ Ordnung in der verschiarften 
in der Nahe von 6,= 06, +;. Theorie fiir den Schlankheitsgrad 
Resonanzverhaltnis 7 0. A = 30. 


viel starkeren EinfluB, so da ein gréBerer Schwellwert ¢, zur Anregung einer Instabilitat er- 
forderlich sein wird, als in Bereichen erster Art. 

Das MaB der Instabilitat in den beiden Bereichen ist nach den Abb. 8 und 10 allerdings nicht 
‘ohne weiteres zu vergleichen, denn in Bereichen erster Art gibt es nur jeweils einen charakte- 
ristischen Exponenten, welcher mit seinem Imaginarteil fiir die Instabilitat maBgebend ist. In 
Bereichen zweiter Art gibt es dagegen zwei derartige Exponenten, denn fiir cos 270 findet man 


ein konjugiert komplexes Wurzelpaar. AuBerdem gehen in die Naherungsliésung (32) fiir w 


1 Vel. S$. 441 der in FuBnote 5 von S. 162 genannten Arbeit. 
13* 
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| 


nach (33) von den Anfangshedingungen abhangige Integrationskonstanten ein, tiber die keine) 


allgemeinen Aussagen gemacht werden kénnen. 


y) Bereich zweiter Art erster Ordnung; Verschirfte Theorie 1 ~ 0 (Schlankheitsgrad 2 = 30)., 
Betrachtet man auch sehr gedrungene Stabe, so sind Instabilitaten zweiter Art in der Nahe» 
simtlicher Schwebungsfrequenzen im Gegensatz zur vereinfachten Theorie méglich [vgl. Ab-- 
schnitt 3c) B)]. Die niedrigste Frequenz, welche dann zu einem Instabilitatsbereich fiihrt, ist) 
5, = 6,-+ 6,. In Abb. 12 und 13 sind die zahlenmaBigen Ergebnisse wiedergegeben. Besonder-- 
heiten gegeniiber den Bereichen zweiter Art der vereinfachten Theorie sind jedoch nicht zu be-} 
merken. 

6) Bereich erster Art in der Nahe der Lingsresonanz. Die verscharfte Theorie gestattet einen 
Ausblick auf sehr merkwiirdige Erscheinungen, die eintreten, wenn die Anregungsfrequenz die} 
Langsresonanzfrequenz erreicht. Wie in Abschn. 4a) mit Hilfe der Transformation (36) gezeigt | 
wurde, ist der in der Nahe der Langsresonanz schwingende Stab gegen eine Querauslenkung} 


cos 2¢@,| 


12Z- instabi/ 


10 “ Vf, 
470 €-G0- G25 Ul ] 
stabil 4 Z| 4 
h dé Y d 
ts 30 35 40 45 50 O76 30 | 36 3b 42 «+46 5D 
0-26, ORes Op=20; Shes 
Abb. 14. | cos2 0,| abhangig yon 6 mit dem Parameter ¢ in Lingsresonanznihe. Abb. 15. Bereich erster Art in Langsresonanzniahe 


instabil. Dies Verhalten ist auch mechanisch unmittelbar einleuchtend. Demnach muB es bei} 
den Langsresonanzfrequenzen 6; xs, zusatzliche Instabilitatsbereiche geben. Es wurde versucht, 
einen Instabilitatsbereich erster Art (69= 20,= 32,20) unmittelbar unterhalb der ersten Langs4 
resonanzfrequenz (Oz res. == 47,12) mit Hilfe der Determinantentheorie zu berechnen, obgleict 
die Grundlagen der Rechnung in Resonanznahe véllig unsicher sind, weil eF,,,, nicht mehr kleiz 
ist. Aus der Auftragung in Abb. 14 wurden die Grenzkurven der Abb.15 gewonnen. Dem Er 
gebnis kommt keine zahlenmafige Bedeutung zu. Es soll lediglich behauptet werden, da Er 
scheinungen der skizzierten Art méglich sind. Hierauf deutet auch schon der mit der Stérungs+ 
theorie berechnete Bereich der Abb. 5 in Abschnitt 4c) beim Stab ohne Vorlast hin. Die dort 
angewandte Methode der Grenzkurvenberechnung durch eine Potenzreihenentwicklung kansi 
jedoch den iiber der Langsresonanzfrequenz liegenden Teil des Bereiches nicht erfassen. 


Eine auch zahlenmabig befriedigende Beschreibung der Instabilitat in der Nahe der Langs 
resonanz wird erst méglich sein, wenn die asymptotische Integration der Stabilitatsgleichungen 
(25) fir ¢ > etwa in dem Umfange bekannt ist wie bei der Mathieu-Einzelgleichung. 


6, Zusammenfassung, Die vorliegende Arbeit untersucht die Stabilitat der Querbewegunger 


eines eingespannten, achsial pulsierend belasteten Stabes. Zu dem Zweck wird abhangig vom 
Frequenz und Amplitude untersucht, ob bei einer beliebigen Anfangsstérung die Querbewegung, 
unbegrenzte Amplituden annehmen kann oder nicht. In Ermangelung eines geeigneten Ansatzeq 
wird eine geschwindigkeitsproportionale Dampfung zugrunde gelegt. Die rechnerische Behand4 
lung wird iiber ein Variationsproblem durchgefiihrt, wobei das Variationsproblem durch einen 
Ritz-Ansatz mit zeitabhangigen Koeffizienten genahert gelést wird. Das auf diesem Wege ent! 
stehende gewoéhnliche Differentialgleichungssystem mit periodischen Koeffizienten kann als Ver 
allgemeinerung der Mathieuschen Differentialgleichung angesehen werden und dient zur Be+ 
urteilung der Stabilitat. Sieht man von der Dampfung zunichst ab, so kann mit der allgemeinen 
Theorie von E. Mettler+ Grenzkurven in der (0, ¢)-Ebene berechnen, welche diejenigen zusammen} 
gehérigen Werte von Frequenz und Amplitude einschlieBen, die zu instabilen Bewegungsformen 
fiihren kénnen. Um auch die Dimpfung und einige bisher ausgeschlossene Sonderfalle in di¢ 


1 Vgl. FuBnote 5 von S. 162. 
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Betrachtung einzubeziehen, werden die zur Beurteilung der Stabilitat dienenden charakte- 
ristischen Exponenten des Differentialgleichungssystems mit Hilfe unendlicher Determinanten 
unmittelbar berechnet, und daraus die Grenzkurven abgeleitet. 

Die Behandlung des Problems wird auf zwei verschiedenen Stufen der Naherung durchgefiihrt. 
Die bisher allein gebrauchliche erste Naherung unterdriickt den Resonanzfall véllig, gestattet 
aber einen vollstandigen Uberblick iiber die Erscheinungen. An Einzelbeispielen wird gezeigt, 

“wie die verscharfte Theorie die Resonanzfrequenzen der Lingsschwingungen beriicksichtigt und 
zu weiteren Instabilitatsbereichen fiihrt. Die Abgrenzung der Naherungen gegeneinander er- 
gibt, da die vereinfachte Theorie nur fiir extrem schlanke Stibe zahlenmaBige Bedeutung hat. 

Fir die technisch wichtigsten Instabilitatsbereiche des Stabes mit pulsierender Langslast 
ohne ruhenden Anteil werden sehr einfache Naiherungsformeln angegeben, welche einen raschen 
Uberblick iiber das Stabilitatsverhalten technischer Stabe gestatten. 


Aus dem Institut fiir Mathematik und Mechanik der Bergakademie Clausthal. 


(Eingegangen am 7. Oktober 1950.) 
Anschrift des Verfassers: Dr. F. Weidenhammer, Clausthal-Zellerfeld, Schulstr. 2. 
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Uber den Mechanismus der Wirbelkerne einer Karmanschen Wirbelstrabe. 


Von W. Kaufmann, 


1, Einleitung, Wird ein prismatischer Kérper senkrecht zu seiner Achse gleichformig mit der 
Geschwindigkeit U in einer ruhenden, wenig zahen Fliissigkeit bewegt, so bildet sich bekanntlich 
hinter dem Kérper ein mehr oder weniger stark ausgepragtes Wirbelsystem aus. Bei geeigneten 


Werten von U und entsprechenden Kérperabmessungen ordnen sich dabei die einzelnen Wirbel | 


in bestimmten Abstinden auf zwei nahezu parallelen Geraden an. Sie bilden, wie man sagt, 


eine WirbelstraBe, die dem Korper, als Ganzes betrachtet, mit der kleineren Geschwindigkeit 
u< Ufolgt. Zwischen den beiden Wirbelreihen stellt sich eine pendelnde Bewegung der Fliissig- | 


keit ein. 
Diese in Wasserrinnen beobachtbare Erscheinung veranlafte Th. v. Karman, die Stabilitat 


einer derartigen Wirbelkonfiguration genauer zu studieren’. Dabei werden folgende ideali- | 


sierende Annahmen gemacht: 
1. Die Fliissigkeit soll reibungsfrei und raumbestandig sein. 


2. Der prismatische Kérper (rechteckige Platte, Zylinder od. dgl.) wird senkrecht zur 


Strémung als unendlich lang vorausgesetzt (ebenes Problem). 


3. Die einzelnen Wirbelfaden stehen senkrecht zur Strémungsrichtung und werden als Po- | 


tentialwirbel eingefiihrt, deren Drehsinn auf beiden Seiten der Strabe entgegengesetzt ist. 

4. AuBerhalb dieser Wirbelfaiden (die sich in der Bildebene als ,,Wirbelpunkte‘‘ projizieren) 
herrscht Potentialstrémung. 

5. Die WirbelstraBe wird als beiderseits unendlich lang angenommen. Die Frage nach der 
Entstehung der Wirbel wird nicht diskutiert, die Existenz der Wirbel vielmehr als eine physi- 
kalische Realitat angesehen. 

Unter den vorstehenden Voraussetzungen findet v. Karman, daB nur die aus Abb. 1 ersicht- 
liche Wirbelanordnung gegen kleine Stérungen der Anfangslage stabil ist, bei welcher die Wirbel- 
faden der beiden Reihen um die Strecke 1/2 gegeneinander verschoben sind?. Dabei ergibt sich 
fiir das Verhaltnis zwischen dem Abstand h der beiden Wirbelreihen und dem seitlichen Wirbel- 


abstand / (Wirbelteilung) der von der speziellen Querschnittsform des Kérpers unabhangige 
Wert 


h 
7 = 0,281. (1) 


Demgegeniiber kann man bei Laboratoriumsversuchen haufig trotz sorgfaltigen Experimen- 
tierens bereits vom dritten oder vierten Wirbelpaar ab ein Zerflattern der WirbelstraBe beob- 
achten, was darauf schlieBen 1aBt, daB irgendwelche Einfliisse vorhanden sind, welche die fiir 
die ,,ideale‘‘ Strémung nachgewiesene Stabilitat stéren. Die Ursachen hierzu sind noch nicht 


hinreichend geklart; im allgemeinen wird die Fliissigkeitsreibung dafiir verantwortlich gemacht. | 


Man kann nun zeigen, daB diese auch bei wenig zihen Fliissigkeiten nicht in Erscheinung tritt, 


sofern man, wie oben geschehen, die Bewegung der WirbelstraBe als Potentialstrémung ansieht. | 


Fiir die raumbestandige Fliissigkeit lautet die Kontinuitatsbedingung 
Ovx , Ovy _ 


0x aie Oy Of (2) | 
AuBerdem ist wegen der im Gebiet auBerhalb der Wirbelfaden vorausgesetzten Wirbelfreiheit | 
Ovx Ovy 5 aoe 
eo Ve (3) 


wenn v, und vy die Geschwindigkeitskomponenten nach den Koordinatenachsen bezeichnen. | 


: 1 Th, v. Karman, Nachr. Ges. Wiss. Gottingen 1911/12 und v. Karman u. Rubach, Physik. Z. 13 (1912) 
49 


> Auf die Méglichkeit, daB bei ganz speziellen Stérungen auch die Wirbelanordnung nach Abb. 1 instabil | 


sein kann, hat C. Schmieden hingewiesen, Ing.-Arch. 7 (1936), S. 215. 


° Vel. hierzu die von A. W. Maue, Z. angew. Math. Mech. 20 (1940), S. 130, angegebene Verallgemeine- 
rung des Kdrmdnschen Stabilitatskriteriums. 


| 
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Differentiiert man (2) nach x und (3) nach y, und addiert die so entstehenden Ausdriicke, dann 
ergibt sich , 


=. Avy 0. (4) 


Kine entsprechende Beziehung gilt fiir die y-Richtung. Die Werte 4v, und Av, sind aber gerade 

die GréBen, in denen sich die Navier-Stokesschen Gleichungen der zaihen Flissigkeit von den 
Euler-Gleichungen der idealen unterscheiden. Im Innern der Flissigkeit, d.h. fern von allen 
Wanden, treten demnach die Reibungskrafte dynamisch iiberhaupt nicht in Erscheinung. Die 
bei ,,natiirlichen“ Flissigkeiten tatsachlich beobachtete Stérung der Stabilitat 14Bt also darauf 
schlieBen, da die vorausgesetzte Potentialstrémung sich auf die Dauer gar nicht einstellt, daB 
vielmehr entweder von den Randern der Fliissigkeit oder von den Wirbelfaden — fiir welche ja 
Gleichung (4) nicht gilt — Stérungen ausgehen, welche das ,,ideale‘* Strémungsbild allmablich 
verandern. 


Den Einflu8 von Wanden, die zu den Wirbelachsen senkrecht stehen, hat K. Schlayer! unter- 
sucht; er findet, daB eine schwache Instabilitat auftritt, wenn der Abstand der zu den Wirbel- 
achsen senkrechten Begrenzungsebenen gréfer ist als 1/100 der StraBenbreite h. Die quanti- 
tativen Ergebnisse seiner Theorie kénnen jedoch nicht als endgiiltig angeschen werden, da die 
von ihm angestellte Betrachtung iiber das Energiefeld der WirbelstraBe fiir eine ,,natiirliche*, 
wenn auch nur wenig zahe Fliissigkeit den wirklichen Verhdltnissen nicht gerecht wird. 


Die vorliegende Untersuchung beschaftigt sich nicht mit dem eventuellen EinfluB der Wande 
auf die Stabilitat der StraBe (Einflu®B der Gerinnetiefe und Gerinnebreite). Sie soll vielmehr 
der Frage gelten, ob das theoretische Strémungsbild etwa dadurch wesentlich gestért wird, daB 
in einer ,,natiirlichen’ Fliissigkeit keine ,,Potentialwirbel‘* auftreten, sondern Wirbel von end- 
lichem Querschnitt. 


Nach der Prandtlschen Grenzschichttheorie sind derartige ,,Wirbel‘+ das Endprodukt einer 
sich aufrojlenden Trennungsschicht, die sich infolge der Reibungseinfliisse in einer schmalen 
Zone (Grenz- oder Reibungs- 
schicht) in unmittelbarer Nahe 
des umstrémten Kérpers ausbil- 
det?. Eserhebt sichjetztdie Frage ©. 
nach der Gréfe dieser Wirbel, 
um dadurch beurteilen zukénnen, 
wie weit gegebenenfalls dasStrom- 
linienbild der Karmanschen Ideal- 
strémung von dem der wirklichen Abb. 1. Schema der Wirbelstrabe. 

Strémung abweichen wird. Nach- 
stehend soll der Versuch gemacht werden, mit Hilfe des Energiesatzes etwas ither die GréBe 
dieser Wirbel auszusagen. 


2. Die kinetische Energie des Geschwindigkeitsfeldes der WirbelstraBe, Die durch die Wirbel- 
straBe erzeugte Potentialstrémung ist instationaér. Ihr komplexes Strémungspotential, bezogen 
auf das aus Abb. 1 ersichtliche Koordinatensystem, lautet mit 


I AIO 
n=i+it (5) 


w= ptip= in — (6) 


worin (x, y) das Geschwindigkeitspotential, w(x, y) die Stromfunktion und J die Zirkulation 
der einzelnen Wirbel darstellen, wahrend z— x-} iy einen beliebigen Punkt des Feldes an- 
gibt 3, Spaltet man diesen Ausdruck in seinen reellen und imaginadren Teil auf, so liefert der 


K. Schlayer, Z. angew. Math. Mech. 8 (1928), S. 352. 

L. Prandtl, Fihrer durch die Strémungslehre, 3. Aufl., $.127. Vgl. hierzu auch A. Betz, Naturwiss. 
1950), S. 193. 
Vel. 


37 
1. etwa W. Kaufmann, Angew. Hydromechanik, 1. Bd. (1931), 5.193. 
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imaginére Teil als Stromfunktion der WirbelstraBe den Wert 
; l\ a ' h\ a Nall as h\ x 
r, sin? (. Fr ‘| 4 of? (y de 5) i + cos? (: = i i Sin? (y+ 4 d 


ae esi leor|o—2) 9] ome Dae oD) 


Abb. 2 zeigt einen Ausschnitt des Stromlinienbildes einer stabilen WirbelstraBe nach v. Karman. 
Daraus ist ersichtlich, daB das Strémungsfeld in zwei verschiedene Bereiche zerfallt: das Gebiet _ 
um die Wirbelpunkte (kreisende Bewegung) und das zwischen ihnen liegende Gebiet der pen- | 
delnden Bewegung. Die Stromlinien, welche beide Bereiche begrenzen, sollen hinfort als Grenz- 
stromlinien bezeichnet werden. 

Fir den innerhalb einer oberen Grenzstromlinie liegenden Wirbelpunkt mit den Koordinaten | 


x=1/4+ nl und y=h/2 (n = ganze Zahl) folgt aus Gleichung (7) 
Vd ae 0.8 (8) 


(7) 


Die kinetische Energie der kreisenden Bewegung um einen Wirbelpunkt lat sich unter Be- 
achtung der Abb. 2, bezogen auf die Tiefe ,,eins‘‘, wie folgt darstellen: 


VW o=2n 
Ey = > | vednds, (9) | 
Per o=0 


wobei v die skalare Geschwindigkeit, ds das Lingenelement einer Stromlinie und dn den Ab- } 
stand zweier benachbarter Stromlinien bezeichnen. Das Doppelintegral erstreckt sich iiber den 
ganzen Bereich von der Grenzstromlinie (y,,) bis zum Wirbelpunkt (py). Der Ausdruck 


vdn = dy 
stellt das sekundliche DurchfluBvolumen | 
zwischen zwei Stromlinien dar, das aus 
Kontinuitatsgriinden konstant ist, weshalb 
(9) auch in der Form 
YW b=2n 


B= £ |  { vds (10) 


Ver b=0 


geschrieben werden kann. Beachtet man 
noch, dah 


O=22 
Abb. 2. Stromlinienbild der WirbelstraBe. | vds= 1, 
9=0 
die Zirkulation langs einer Stromlinie um den Wirbelpunkt angibt, welche fiir alle Stromlinien 
des betrachteten Bereiches den gleichen Wert besitzt, so folgt aus (10) 


Ey = 1's (yw — Yer) » (11) 


und dieser Ausdruck wird wegen (8) unendlich gro8. Zu dem gleichen Ergebnis gelangt man 
fiir das Feld jedes der ittbrigen Wirbel. 

Nach dem Energiesatz ist die kinetische Energie der hinter dem bewegten Kérper entstehen- 
den Strémung das Aquivalent der Arbeit, die bei der Vorwartsbewegung des Kérpers geleistet 
wird. Da diese Arbeit durchaus endlich ist, widerspricht das gewonnene Ergebnis den physi- 
kalischen Tatsachen. Die Ursache ist nach dem weiter oben itber die Wirbel Gesagten leicht 
verstandlich. In der Umgebung des theoretischen Wirbelpunktes bildet sich ein Strémungs- 
gebiet aus, in dem die bisher angenommene Potentialstrémung nicht mehr gilt, und dessen Energie 
einen endlichen Wert besitzt. In der Folge soll dieser Bereich als Wirbelkern bezeichnet wer- 
den. Seine GréBe ist zundachst unbekannt. ; 

Es sei nun angenommen, daf$ auBerhalb dieser Wirbelkerne auch weiterhin Potential- 
strémung vorhanden sei, und daB der Rand des Wirbelkernes zugleich Randstromlinie der 
duBeren Potentialstrémung ist. Bezeichnet man die Stromfunktion dieser Randstromlinie mit 
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Yr, 80 kann die Energie der kreisrunden Bewegung auBerhalb eines Wirbelkernes nach (10) 


in der Form 
E,, = ww f ee 


angeschrieben werden, wobei jetzt y, die Wie rae parallel darstellt. Man erhalt dann 
, Q 
E,, = . Ty (y — Wor) « (12) 


Der Wert y, ist zunichst unbekannt. Um die gesamte Energie der kreisenden Bewegung 
um einen Wirbelpunkt zu erhalten, hat man zu (12) noch die Energie Ey des Wirbelkernes 
hinzuzufiigen, so daf 


Ee = 5 Dy (Ys. — Yer) + Ew (13) 

wird. 
Zur Berechnung der Energie der pendelnden Bewegung soll ein Gebiet von der Lange 1 
betrachtet werden, dessen seitliche Rander die Lotrechten x= %/,/-+-nl und x’=#/,I4+-(n+-1)1 


bilden, und das im iibrigen durch die obere und untere Grenzstromlinie eingeschlossen ist. Die 
kinetische Energie dieses Strémungsbereiches ist 


Por B 


E,= 2 | vdnds. (14) 
ae 
Yer 


Darin bezeichnen Wer die Stromfunktion der oberen, Wor diejenige der unteren Grenzstromlinie, 
A und B Anfangs- und Endpunkt einer beliebigen Stromlinie zwischen den beiden Lotrechten, 
die das betrachtete Gebiet seitlich begrenzen (Abb. 2), wahrend ds und dn die gleiche Bedeutung 
haben wie oben. Setzt man wieder vdn= dy, so kann (14) in der Form geschrieben Werdan 


Por B 
E,=2{ ay| vas. (15) 
Var 


B 
Das Integral fv ds hat eine einfache physikalische Bedeutung. Um sie zu erlautern, betrachte 
A 


man die Zirkulation langs des geschlossenen Linienzuges, der von folgenden Linien gebildet 
wird: Stromlinie A B, Lotrechte durch B bis y = oo, Gerade y = © zwischen den Lotrechten 
durch A und B, Lotrechte durch A von y = co bis zum Punkte A. Da fiir y = oc die Geschwin- 
digkeit v= 0 ist, und da die beiden Lotrechten durch A und B alle Stromlinien senkrecht 
schneiden, kann nur das Stromlinienstiick yon A nach B einen Beitrag zur Zirkulation liefern. 
Es mu8 also sein 


B 
fvds=Iy, 
A 


da innerhalb des betrachteten Linienzuges nur ein Wirbel von der Zirkulation /’) liegt. Somit 
erhalt man aus (15) 


E, = "5-15 (Yer — Ver) - 


Wie man leicht feststellt, ist y,— — p,, weshalb 
tse a dad Ai pee (16) 
Fat man schlieBlich die gesamte Energie der kreisenden und der ape Bewegung in dem 
Gebiet zwischen den beiden Geraden x= 3/,1/-+ nl und x’ = 3/,1-++ (n-+-1)l zusammen und 
beachtet, daB darin die Energie von z w ei Wirbeln enthalten ist, so erhalt man aus (13) und (16) 
E=2E,+ £E,= ol oy,+ 2 Ey. (17) 


Zur zahlenmiBigen Berechnung dieser Energie bedarf es also nur der Kenntnis der Energie des 


196 Kaufmann: Uber den Mechanismus d. Wirbelkerne einer Karmanschen Wirbelstrabe. Ingenieur-Archiv — 


Wirbelkernes sowie der Stromfunktion des Kernrandes. Beide Gré8en hangen miteinander | 


zusammen. 


| 


3. Der Wirbelkern. Wie bereits unter Ziffer 1 bemerkt wurde, enthalten die Wirbelkerne 
die spiralférmig aufgespulten Unstetigkeitsflachen, die infolge Grenzschichtbildung an dem um- | 
strémten Hindernis (Platte, Zylinder od. dgl.) standig neu entstehen. Nimmt man zunachst | 
an, daB die Querschnitte der Wirbelkerne unserer WirbelstraBe klein gegeniiber dem Gesamt- | 
bereich der kreisenden Bewegung um die theoretischen Wirbelpunkte sind, so kénnen die Kern- 


rander angenahert als Kreise aufgefaBt werden. Unter dieser Voraussetzung laBt sich die Energie 
eines Wirbelkernes in der Form darstellen + 


omar | 
£55, 8 ? 18 
in | iaade oe 
r=0 


wobei ry den Halbmesser des Kernrandes, r denjenigen einer Stromlinie innerhalb des Kerns | 


und J’, die mit r veranderliche Zirkulation lings dieses Kreises bezeichnen. Fiir r= 0 ist thee! 
fiir r= r, dagegen gleich J), d.h. gleich der Gesamtzirkulation eines Wirbels (s. oben). Der 
Verlauf I, = T°,(r) ist nicht bekannt. Man kann daritber auBer den vorstehenden Randbedin- 
gungen nur vorschreiben, daB [dJ/dr],.,, = 0 sein muB, da fiir r= r, die Zirkulation nicht 


mehr wachst. Aus der Untersuchung der Wirbelkerne, die sich hinter Tragfliigeln einstellen, | 


ist bekannt, daB [dI)/dr],,.9—> co geht®. Sieht man diese Bedingung auch im vorliegenden | 
Falle als giiltig an, so ist der einfachste Ansatz, den man fiir /,(r) machen kann, eine tiber den | 


Halbmesser r elliptisch ansteigende Zirkulationsverteilung. Es mag also dafiir der Ansatz 


P=Tol/27, =leih (19) 


gewahlt werden, welcher jedenfalls die vier oben gestellten Bedingungen erfiillt. Damit lautet (18) | 


oder 


ae (20) 


Das bedeutet, daB die Energie eines derartigen Wirbelkernes nicht von seiner GroBe, sondern 
nur von der Zirkulation J’, abhangt. Mit diesem Ausdruck fiir Ey liefert Gleichung (17) als 
gesamte Energie des betrachteten Bereiches 


E=0D(¥+73): (21) 


Die Relativgeschwindigkeit des durch dic Fliissigkeit bewegten Korpers gegen die WirbelstraBe | 
ist U—u. Demnach betragt die Zeit, welche verstreicht bis die WirbelstraBe sich um die | 


Lange / vergréBert hat, 
l 


T= _—_. (22) 


U—u 
Wahrend dieser Zeit bewegt sich der Kérper absolut um die Strecke 


Ul 
U—u 


i,= UT. = 


vorwarts. Dabei leistet die zu seiner Bewegung erforderliche Kraft P die mechanische Arbeit Ploml 
welche, sofern vorerst von allen Reibungseinfliissen abgesehen wird, das Aquivalent der Energie | 


(21) darstellt. Man erhalt somit 
Tk 
Feelin ta (24) 


1 W. Kaufmann, Ing.-Arch, 17 (1949), S. 190. 
2 W. Kaufmann, Sitzungsber. Bayer. Akad. Wiss., Mathem.-Naturwiss. Abt. 1946, S. 109. 


(28) | 
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Bezeichnet ¢,, die ,,Widerstandsziffer‘‘ des bewegten Kérpers und b seine Hihe (senkrecht zur 
Strémung), so kann P, bezogen auf die Tiefe ,,cins‘‘, bekanntlich in der Form angeschrieben 
werden 


se Ue 
Pog. Uh, (25) 
wahrend /, mit der Abkiirzung 
u 
Go (26) 
die Form 
l 
lb={=n (27) 
annimmt. Mit diesen Werten lautet Gleichung (24) 
a) l 3 I, 
of 02 b= oly. 4g). (28) 


Es mége nun x, die Abszisse des Punktes C bezeichnen, in dem der linke Rand des Wirbelkernes 
die durch den theoretischen Wirbelpunkt gehende Horizontale y = h/2 schneidet. Dann erhalt 
man als Stromfunktion der durch C gehenden Randstromlinie aus Gleichung (7) mit x= x, 
und y= h/2 den Wert 


ra, lect Dgfeor tts) + son ant Deol 


» == —2] 29 
Y 4m : “oe ( ; 1 ( ) 
sin x, = 4 7 
Setzt man noch zur Abkiirzung 
h : h 
Coj? ; n) ore Sin? (3 n| 16, = 65 (30) 
und 
l\ x V\ x I 
letgi=a, (ue —gf=a—3 a 
und fiihrt (29) in Gleichung (28) ein, so wird 
@ re OIE c, sin? « + ¢,cos*a , 
Wy el arenas Ax (In cos? & el () 


Nach v. Karman? besteht zwischen der Fortschrittsgeschwindigkeit u der WirbelstraBe und der 
Zirkulation J, die einfache Beziehung 


10) 


g . 
1/8 


Damit wird wegen (26) 


8 n2l2’ 
so da (32) iibergeht in 
Cig I b as 2 
TE ea 3 == In) (e717 0-1-6) < (33) 
Setzt man die Konstante : 
Cry TE b 
Ww =. _ 34 
4n?(1—n) 1 2 ie ee) 


so folgt aus (33) 
e* = ¢, tg? a+ ¢, 
und somit 


Paacg No 
tga = ee (35) 


Cy 


1 Th. v. Karman, a. a. O. 
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Aus dieser Gleichung kann zunachst « und wegen (31) x, berechnet werden, wenn ¢,, n, b und 1 
bekannt, bzw. durch Messung bestimmt sind. 
Bei den Karmdnschen Versuchen wurden eine senkrecht zu ihrer Ebene bewegte Platte von 


der Hihe b— 1,75 cm und ein Kreiszylinder vom Durchmesser b= 1,5 cm verwendet. Die | 


entsprechenden Versuchswerte waren: 
Platte: n= 0.2... be= 1,15 em, Wes. 9,8 emese eos 
Zylinder: n= 0,14; b= 1,50 cm, 1 64cm, c= 0,92; 
Damit erhalt man aus (34) fiir die Platte x, = 4,02 und fiir den Zylinder x, = 7,04. Nach (30) 


ist unter Beachtung von Gleichung (1) c,= 2,0 und c,= 1,0. Aus (35) folgt also fiir die Platte | 
tga —= + 5,22. Diesem Werte entsprechen die beiden Winkel «, = 1,38 und ag =a — o,= 1,168) 


Aus (31) erhalt man damit als Abszisse des Punktes C (Abb. 2) x, = 1,86 cm und als Abszisse 
des rechts vom Wirbelzentrum liegenden Punktes C’ des Kernrandes x, = 3,04cm. Demnach 
wird fiir. die Plattenstrémung der horizontale Kerndurchmesser dp = x, 38, == 1 Bie eee 
entsprechenden Zahlenwerte beim Zylinder ergeben sich wie folgt: 


teo + 23,88, 0, = 1,5287, o«,—1,6129, x,=—1,51, «%=—1,69, dz=0,18cem. 


Da nun die Breite der Wirbelstrafe fiir die Platte theoretisch nach (1) hp = 0,281 1 = 2,75 cm 
und fiir den Zylinder 1,8 cm ist, so wird dp = 0,43 hp, dz = 0,1 hz. 

Die vorstehende Theorie liefert also fiir die Platte einen wesentlich gréBeren Wert d/h als 
fiir den Zylinder. Unter diesen Umstinden erscheint es fraglich, ob fiir die Platte die oben 
gemachte Annahme kreisférmiger Kernquerschnitte und der darauf beruhende Ansatz (19) fiir 
die Zirkulationsverteilung noch zulassig sind. 

Von Interesse ist in diesem Zusammenhang die Feststellung, da die von v. Karman und 
Rubach experimentell gefundene Breite h der WirbelstraBe beim Zylinder sehr gut mit dem 
theoretischen Wert nach Gleichung (1) iibereinstimmt, wahrend fiir die Platte als Mittelwert 
hp = 3,0 cm, unmittelbar hinter dem Kérper sogar 3,5 cm, festgestellt wurde, im Gegensatz 
zu dem mit / == 9,8 cm aus (1) berechneten Wert von 2,75 cm. Es scheint danach, dai die Theorie 
die wirklichen Verhaltnisse fiir den Zylinder besser wiedergibt als fiir die Platte. Und das kénnte 
damit zusammenhiangen, daf eben beim Zylinder der kleinere Kerndurchmesscr eher dem Ideal- 
bild des Potentialwirbels entspricht als bei der Platte. Der Unterschied in den Verhaltniswerten 
d/h fiir Platte und Zylinder wird iibrigens erheblich gemildert, wenn man als c,,-Wert der Platte 
nicht den aus der Kdrmdnschen Widerstandsformel berechneten, sondern den fiir ebene 
Strémung giiltigen Wert 2,01 einfiithrt’. Mit diesem Wert liefert die obige Rechnung xp = 5,82 
und damit dp = 0,48 cm, oder dp/h = 0,175. 

Behalt man trotz der vorstehend geaduBerten Bedenken zunachst einmal Gleichung (20) bei, 
wonach die Kernenergie lediglich von der Zirkulation I’) abhangt, so kann die Energiebedingung 
(32) wie folgt geschrieben werden: Die linke Seite dieser Gleichung stellt, wie aus ihrer Herleitung 
folgt, die Arbeit PJ, dar, wahrend der erste Summand der Klammer wegen (33) und (34) gleich x 
ist. Somit wird 


el 
Pl, = 4 (x+ 3). 


Mit den oben errechneten x%-Werten ergibt sich daraus fiir die Platte 
els 


Pl, = £—* (4,02 + 3) 


bzw., wenn mit ¢, == 2,01 gerechnet wird, 


Pl, = 228 (5,22 4 3), 


und fiir den Zylinder 


Pi, = £22 (7,044 3), 


wobei der zweite Summand jeweils den Beitrag der Kernenergie darstellt. Man erkennt daraus, 
da®B bei der Platte von der geleisteten Arbeit 43% (34%) in Kernenergie umgesetzt wird, beim 
Zylinder dagegen nur 30%. 


> Velsiuttee Bde 2ieAutlemonso3. 


| 
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Um jetzt noch einmal auf die in Ziffer 1 aufgeworfene Frage nach der Stabilitit einer ,,physi- 
kalischen* WirbelstraBe zuriickzukommen, kann beziiglich der Wirbelkerne folgendes gesagt 
werden: Bei ebenen, reibungsfreien Strémungen in konzentrischen Kreisen kénnen die Zentri- 
fugalkrifte je nach ihrer radialen Anderung stabilisierend oder destabilisierend wirken. Bezeich- 
net u den Mittelwert der tangentialen Geschwindigkeit, so hat man es nach Prandtl? mit Stabili- 

tat zu tun, wenn das Produkt ur nach auen hin zunimmt, im anderen Falle mit Instabilitat. 
_ Bei Potentialstrémung mit Zirkulation herrscht indifferentes Gleichgewicht. In den hier be- 
trachteten Wirbelkernen ist ur—TJ’./2, und dieses Produkt nimmt nach aufen hin bis zum 
Kernrand zu. Der Wirbelkern stellt demnach ein st abiles Gebilde dar. Die Tatsache allein, 
da in einer physikalischen WirbelstraBe derartige Wirbelkerne an Stelle der Potentialwirbel 
auftreten, kann also kein Grund fiir das haufig beobachtete Zerflattern der WirbelstraBe sein. 

Andererseits mu8 festgestellt werden, da der in Ziffer 1 gezogene SchluB beziiglich des 
Wegfalls aller Reibungseinfliisse in der 4uSeren Potentialstr6mung fiir die Kerne nicht mehr zu 
trifft, da in den Kernen Gleichung (3) und somit auch (4) keine Giiltigkeit haben. 

Damit erfahrt aber auch die weiter oben zunachst fiir reibungsfreie Bewegung aufgestellte 
Energiebilanz eine gewisse Anderung, und zwar dergestalt, da% durch die in den Kernen wirk- 
same Fliissigkeitsreibung standig Bewegungscnergie aufgezehrt wird. Die Folge davon diirfte 
eine allmahliche Veranderung des theoretischen Stromlinienbildes sein, die schlieBlich zu einer 
vollstandigen Auflésung der Wirbelstrafe fiihrt. Die Méglichkeit, daB dieser ProzeB durch Sté- 
rungen, die von den Begrenzungsflachen der Fliissigkeit ausgehen, mehr oder weniger stark be- 
schleunigt wird, soll durch die vorstehende Betrachtung natiirlich nicht ausgeschlossen werden. 


(Eingegangen am 19. Oktober 1950.) 
Anschrift des Verfassers: Professor Dr. W. Kaufmann, Miinchen, Technische Hochschule. 


1 L, Prandtl, EinfluB der stabilisierenden Krafte auf die Turbulenz. Vortr. aus dem Gebiete der 
Aerodynamik und verwandter Gebiete, 5S. 5. Aachen 1929. 
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Uber die Beulsicherheit von Reckteckplatten mit querverschieblichen 
Randern. 


Von S. Woinowsky-Krieger. 


1, Einleitung. Die bisher durchgefithrten Stabilitatsuntersuchungen von Rechteckplatten 
beziehen sich vorwiegend auf Fille, wo die dem AuSendruck ausgesetzten Plattenrander entweder | 
als gelenkig gestiitzt oder als eingespannt vorausgesetzt werden. In der technischen Praxis | 
kommt es aber vor, daf ein in der Ebene der Platte belasteter Rand keine Stiitzung senkrecht zu | 
dieser Ebene besitzt. Die Stabilitat einer Rechteckplatte wird nun durch die fehlende Stiitzung 
an einem oder gar an mehreren ihrer Rander so stark beeintrachtigt, daB eine genauere Unter- 
suchung ihrer Beulsicherheit unter diesen besonderen Grenzbedingungen wohl angebracht er- | 
scheint. } 

Wir nehmen im folgenden eine isotrope, diinne Platte an; ihre kritische Belastung mége 
Druckspannungen verursachen, die unterhalb der Elastizitatsgrenze des Materials liegen. Wir | 
wollen uns auf Stiitzungsfalle beschranken, die der strengen Behandlung zuganglich sind und 
beginnen gleich mit einem fiir das vorliegende Aufgabengebiet besonders aufschluBreichen Fall. | 


2. Rechteckplatte mit zwei gelenkig gestiitzten und zwei querverschieblichen Randern. Die 
Rander x= 0 und x— a der Platte (Abb. 1) mégen gelenkig gestiitzt sein, die Rander y = + 6/2 
seien querverschieblich. Ist p die parallel zur y-Achse wirkende Druckbelastung je Langeeinheit 

des Randes und bezeichnet weiter D die Biegesteifigkeit der Platte, so 


a lautet die partielle Differentialgleichung ihrer verbeulten Mittelflache | 
paw 
pe Ee ee ge (1) 


2 2 
worin unter /\ = fate der Laplacesche Operator zu verstehen 


Sle 
ist. Wir befriedigen die Grenzbedingungen an den Randern x= 0, 
a, indem wir von dem iiblichen Ansatz 
mmx 
et 1) SS, Sia) ae (i ei oa) (2) 
Gebrauch machen, worin unter Y eine Funktion von y allein zu ver- 
je stehen ist. Fithren wir den Ausdruck (2) in die Gleichung (1) und 
sate eae setzen noch zur Abkiirzung k,,—= pa?/m? x? D, so erhalten wir fiir Y 


die totale Differentialgleichung 
Abb. 1. Platte mit zwei querver- ates rie 
schieblichen, einem gleichmaGigen IV m* It m* It 
Yn + (lm—2) Yad Ye 08 (3) 


Druck ausgesetzten Randern. 2 m —| 4 
a 


Ist einmal die kritische Belastung fiir den ersten Eigenwert m= 1 ermittelt, so kénnen wir 
spater, wie aus der Gleichung (3) hervorgeht, zu m= 2, 3, ... tibergehen, indem wir einfach 
den Wert a durch a/m ersetzen. Wir diirfen daher, ohne die Allgemeingiiltigkeit der Ergebnisse 
einzuschranken, zunachst m — | annehmen und setzen dementsprechend auch 


2 
k= =. (4) 
Den Lésungsansatz fiir die Differentialgleichung (3) wahlen wir in der Form 
i Tb Hi y 
Y= eC en” (5) 


i=1 
mit den vier Integrationskonstanten C; und erhalten zur Bestimmung der zugehérigen Werte 
des Faktors x; die Gleichung 


wit (k— 2)47+1=0. (6) 
Um den vier Freiwerten C; entsprechend, ebensoviele Gleichungen anzuschreiben, bendtigen 
wir zweiGrenzbedingungen langs jedem der beiden Rander y= + b/2. Da das Biegemoment m, 
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hier jedenfalls verschwinden muB, hei®t die erste dieser Grenzbedingungen 
0? w aw 
viene gta Mais 79 (7) 
worin y die Querdehnungszahl bedeuten soll. 

Um die zweite Randbedingung zu formulieren, benétigen wir eine Vereinbarung iiber das 
Verhalten der am Rande angreifenden Krafte im Augenblick, wo sie den kritischen Wert erreicht 
haben und wo das Ausbeulen der Platte gerade beginnt. Nehmen wir z. B. an, da diese Krafte 
die infinitesimale Drehung der Platte um die beiden Geraden y = + 6/2 mitmachen, so miissen 


wir die Grenzbedingungen des freien Randes anschreiben, d.h. es wird dann die Auflager- 
reaktion 


3 
ate 
Man kénnte unschwer nachweisen, dafi die eben angeschriebene Randbedingung des freien 
Randes in Verbindung mit der Bedingung (7) auf keine reellen Werte fiir die kritische Belastung 
unserer Platte fiihren wiirde +. Die vom technischen Standpunkt aus sinnvollere Annahme ist 
die, daB die 4uBeren Krafte ihre urspriingliche Richtung auch an der in die benachbarte Gleich- 
gewichtslage iibergegangenen Platte beibehalten. Unter dieser Voraussetzung laBt sich die an 
dem Plattenrande angreifende Belastung p nach den Richtungen der Tangente und der Normale 
zu der deformierten Mittelflache der Platte zerlegen. Die letzte der beiden erwahnten Kompo- 


ee 


—0. (8) 


nenten spielt dabei die Rolle einer Auflagerkraft. Sie hat in y= — 6/2 die GréBe r,= Ps 5 


kann aber auch andererseits durch die Forme! (8) ausgedriickt werden, wahrend am Rande 
y = - h/2 die beiden obigen Ausdriicke fiir die negative Auflagerkraft gelten. Setzt man diese 
Ausdriicke einander gleich und beachtet noch die frither eingefiihrte Abkiirzung (4), so erhalt 
man die gesuchte Grenzbedingung in der Form 

ew , ew , Bkow _ 
Goda ACES oy a Pade 


On (9) 


Indem wir jetzt auf die charakteristische Gleichung (6) zuriickgreifen, stellen wir ihre Wurzel 
durch die Formel 


(eke 1-44 ia=m: (= y—]) (10) 


dar. Man sieht, daB die Form der Funktionen Y und somit auch die Form der Gesamtlisung (2) 

wesentlich von der GréBe k abhangt. Wir miissen hier drei Méglichkeiten, namlich k < 4, 

k— 4 und k > 4 unterscheiden und beginnen gleich mit dem fiir das weitere wichtigsten 
Fallk< 4. Die vier Wurzeln (10) sind in diesem Fall paarweise konjugiert komplex und 


man kann ihnen die Form 


i 


i 


mit 


+ (A + tH) | 


ce: Vk | ey) 


A ES ——— 


2 ‘ 9 
geben. Die allgemeine Liésung (2) 1aBt sich somit in folgender Weise darstellen: 
A 7 
w= (4 Gin “22 cos eae B oj = sin ery C Cof — wee ay 4 | 
; (12) 
+ D Sint 2% sin me) Sig Oo | 
a a a 


Unterwirft man nun die elastische Fliche (12) denGrenzbedingungen (7) und (9) auf y= + 6/2, 
so erhalt man zur Bestimmung der vier Freiwerte A bis D ein System von ebenso vielen homogenen 
Gleichungen. Infolge der Symmetrie des Spannungszustandes hinsichtlich y= 0 spaltet sich 
dieses System in zwei voneinander unabhangige Systeme auf, von denen das erste nur die Kon- 
stanten 4 und B, das zweite nur C und D enthalt. Der ersten Gleichungsgruppe ist eine in bezug 
auf die Gerade y= 0 antimetrische, der zweiten Gruppe eine symmetrische 


1 Ahnliches gilt z. B. fiir einen Stab unter einer tangential gerichteten Druckkraft an seinem freien 
Ende, siehe A. Pfliiger, Stabilitatsprobleme der Elastostatik, S. 218. Berlin 1950. 


| 
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elastische Flache (12) zugeordnet. Wir gehen jetzt dazu iiber, diese beiden Beulformen im ein- 


zelnen zu betrachten. 
a) Antimetrisches Beulen. Auf dem soeben beschriebenen Wege erhalt man fiir | 


die Konstanten A und B das Gleichungssystem 
A [(#2— w?— v) Ge— 2Aw Cs] +B [(M#—p?—») s+ 2Au Se]—0, | 

A [(kA-+/3 — 3Ap?—2A+ vd) Cc — (ku— pF +32 — 2u-+ rp) Ss] 
+ B[(ku — w2+32u—2u+yp) Ce+ (kA +2 — 3Ap? — 24+ 9) Ss] =0, | | 
worin zur Abkiirzung | 


. Amb 
S= Cin, 
a 


ab Ds eab 
er C= COR 


Gime Cotes Mamaia 
2a 


gesetzt ist. | 

Da die Gleichungen (13) homogen sind, kénnen die beiden Konstanten nur dann von Null | 
verschieden sein, wenn die Nennerdeterminante des Systems verschwindet. Die Erfiillung dieser 
Forderung fiihrt auf die Beulbedingung 


~~: (ab y4—k aut Brida ee | 

[(3-+ ») (l—») k] VE Gin (*PA—E =F [k— —8])4— Fsin( WF = 05 (14) | 
worin im vorliegenden Falle der Antimetrie das obere Vorzeichen zu nehmen ist. Wie 
aus (4) hervorgeht, ist die Beulbelastung bis auf einen konstanten Faktor gleich k, welche GréBe | 
durch (14) vollstandig bestimmt ist. Wir wollen nun den Verlauf des Beulfaktors kin Abhangig- | 


keit vom Seitenverhaltnis der Platte b/a verfolgen. 
Eine offensichtliche Wurzel der Gleichung (14) ist 


ky = (3+) (1—»), (15) 


i 
zee verschwindet, welche letzte Bedingung die Beziehungen — 


sofern zugleich auch sin 
ab yk, = ance (n—— |, 2,3,...) “und bia == nV ky nach sich zieht. 

Unsere k-Kurve schneidet also die Gerade k= k, regelmafig in Abstanden von 2/Viy. Mit 
b/a—co wachst auch die Gin-Funktion in der Beulformel (14) unbegrenzt. Der erste Klam- 
merwert in dieser Formel 
strebt somit gegen Null 
und ergibt k > k,. Mit zu- | 
nehmendem Verhaltniswert 
b/a schmiegt sich also die 
k-Kurve immer mehr an die 
Gerade k= k, an. Im ande- 
ren Extremfall b/a > 0 er- 
halt man aus der Beul- 
formel (14) den Grenzwert 
kj=2(1—v). (16) 
y=0 Die k-Kurven sind fiir 
die beiden extremen Werte 

sty-g3 der Querdehnungszahl 

Y=) ny! Pe (0.3 

berechnet worden und in 
Abb. 2 samt den zugehori- 
gen k,-Geraden dargestellt. 
S=symmetrisch Man erhalt insbesondere im 
A=antimetrisch | Falle v= 0 die Werte k,= 2, 
ofa k,= 3 und fiir die Schnitt- 


0 676 7 1316 2 6030 3 
punkte mit der k,-Geraden 
Abb. 2. Beulbelastung der Platte nach Abb. 1 in Abhangigkeit vom Seitenverhiltnis di . 
b/a fir symmetrische und antimetrische Beulformen, sowie fiir » = 0 und » = 0,3. ie Abszissen 


bla 2 nfy3 = 1,155 n; 


fiir y—=0,3 ergibt sich entsprechend ky= 1,40, kj= 2,31 und 6/a= 2 n/j/2,31 = 1,316 n 
(ae RS ee) 
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auftretende Potenz von M,, und unterdriickt die Ausdriicke mit rat if , so erhalt man 
1 1 
“u—1 
LRM ernest WPOE eee 
ae = -+ — (2% — 1) sin os 
R, Os | R | ) “Eh : i 
—.—] —x cos? 
2 
3 : *x—l1 
1 Aa y (2% + 1) cos Pe ear Ly EOE 
fe a Loa 
: ge ee 5 
R,, On R cos p ¢ He a ; i 
Seat) aw Cony 
Das Verhaltnis der Kriimmungshalbmesser ist hier an jeder Stelle der StoBfront 
33 
eg 2 (2 % + 1) cos*y — 5 (% + 1) 
R, Te Dee 4] sin 29 ‘ 


v 


Ausgezeichnete Stellen sind hier ferner: 


d) Diejenigen Orte an der StoBfront, wo der Ablenkungswinkel 6 sein Maximum 
erreicht. Hierbei gilt 


ae 
R, 6 max R 


1 = 
2tgy Mi cos? p (2+ M+ cos* gp —x M} cost p —> = ™ Mi} cost p +1} ( 


“x+1 
D 


“—l1 
2 


M? cos? y + 1) 


(“+ mt cosy +1) +( 


BF one 7 (tet), = 


1 [ (x + 1) tg? Mj} cost p (“+t M3 cos? p— x M} cost p ——™ My cos? p +1) 


at; ae x—] 2 “+1 

(BS M? cos? p + 1) HAG 
Die 5,,..7 Werte fiir die bei vorgegebener Machzahl M, der Anstrémung der Ablenkungswinkel 6 
gerade sein Maximum erreicht, erhalt man aus Gleichung (61) des Abschnittes 11 durch Null- 
setzen des Zahlers, d.h. 


° 


(M? cos? py — 1) 


\2]3/2 
M? cosy sing) %5 


max 


5 
| 


(M? cos® g — 1) ?6 


max 


23/2 
M? cos ¢ sin 0) | 


L 3 


Lola 1 1 ee | Sahel 


Die fiir verschiedene Machzahlen der Anstrémung M, errechneten Werte 9;,,,, sind in Tabelle 3 
zusammengestellt. 


(30) 


e) Diejenigen Orte an der StoBfront, wo die Strémungsgeschwindigkeit w gleich 
der értlichen Schallgeschwindigkeit a ist. Bei vorgegebener Machzahl der Anstrémung 
M, entspricht der Polarwinkel gy_;=90°—ym-—1 derjenigen StoBfrontneigung yy—1, bei der 
hinter der StoBfront im Strémungsfeld gerade Schallgeschwindigkeit herrscht. Man erhalt die 
@m-i-Werte aus Gleichung (63) des Abschnittes 11, wenn man den Zahler der rechten Seite 
gleich Null setzt, d. h. 


cos? Yy—1 = | 
1 1 ioe ill Vo ee : he oe 
a 8 an + | an Bl al GT a caram te nilspialiecan. [-| 


Die fiir verschiedene Machzahlen der Anstrémung M, errechneten Werte py—; und yy _; sind 
in Tabelle 3 zusammengestellt. Hierbei ist fiir das Verhaltnis der spezifischen Warmen c,/c, der 
Wert x = 1,405 der Rechnung zugrunde gelegt. 

15 
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Tabelle 3. Stofrontneigung fiir gréfte Ablenkung, Schallgeschwindigkeit hinter dem Stofs 
und grépte Stromdichte. 
Pd max ?M =1 - ~* 
M, C0319) = ™, £8 P5 mas [Neugrad] SAIN hea [Neugrad] eee [Neugrad] 
1,0 1,00 1,0 0 1,0 0 1,0 0 
1,05 0,95238 0.98462 11,1811 0,97698 13,6853 
Ueslt 0,90909 0,97153 1552273 0,95757 18,6118 
1S 0,86957 0,96039 17,9775 0,94127 21,9273 
1,2 0,83333 0.95091 20,0298 0,92766 24,3631 0,92269 | 25,1977 
1,3 0,76923 0,93599 22,9008 0,90711 27,6562 
1,4 0,71429 0,92525 24,7707 0.83350 29,6484 0,87376 | 32,3352 
ils) 0,66667 0,91759 26,0264 (0,884.94 30,8394. | 
1,6 0,62500 0,91220 26,8759 0,88000 31,5092 0,84095 | 36,3989 | 
Usa 0,58824 0,90850 27,4459 0,87757 31,8334 
1,8 0,55556 0,90605 27,8161 0,87686 31,9263 0,81678 | 39,1516 
1,9 0,92632 0,904.48 28,0522 0,87727 31,8725 
2,0 0,50000 0,90364 28,1770 0,87842 31,7199 0,79760 | 41,2209 | 
2,20 0,444.44 0,93229 28,2294. 0,88283 SIRL279 } 
25D 0.40000 0,90434 28,0724 0,88790 30,4335 0,76060 | 44,9815 
Ze 0,36364 0,90593 27,8351 0,89244 29,7990 
3,0 0,33333 0,90764 27,9765 0,89651 29,2184 | 0,73092 | 47,8183 
350 0,28571 0,91083 27,0892 0,90297 28,2765 
4,0 0,25000 0,91343 26,6852 0,90762 27,5797 0,68141 | 52,2740 
5,0 0,20000 0,91706 20,1115 0,91393 26,6695 0.64003 | 55.7840 } 
ee) 0,13333 0,92128 25,4289 0,91981 25,6689 
10,0 0,10000 0,92291 25,1603 0,92211 25,2934 0,50358 | 66,4033 
20,0 0,05000 0,92457 24,8856 0,92436 24,9184 0,37684 | 75,4022 
100.0 0,01000 0,92511 24,7944. 0,92510 24,7957 
oo 0 0,92513 24,7906 0,92513 24,7906 
Fiir die Kriimmungshalbmesser der Strom- und Orthogonallinien folgt aus (28a) und (29a}/ 
mite 1 
l 1 ~M} cost p —" = * My cost p + 2 MY cost —* + * My —1 
Bes 2 21 
(a : * 4 T Mz cos p sing |“ Mz cos?» +1) + (22+ mp cos p sing) |” a 
— — =i 


Go =e: | 


f) Stromdichtenverlauf hinter der StoB®front. Fiir das auf die Dichte 9, und aie! 
Schallgeschwindigkeit a, des AusgangsstrO6mungszustandes bezogene Produkt ow hinter der StoB- } 
front kann mit den Gleichungen (46) und (51) des Abschnittes 11 geschrieben werden 


ow M, cos @ = By ie y *x#+1 aro ; yi 
Tea 5— Mi cos*p-- 1 + | 5 M?cospsing) . 
aie kk cos*y + 1 


(32) 


Die Stelle hinter der StoBfront, wo fiir einen gegebenen Gaszustand p,, 0,, M, das Produkt ow) 
sein Maximum erreicht, erhalt man in iiblicher Weise durch Nullsetzen des nach dem Parameter Q? 
differenzierten Ausdruckes (32), d.h. es gilt hier 


d /ow eo d /w w d/o 
d > oa Berroa 
P\Q; % OQ, 4P\qy a, €P\ 


; *x*— 1 ; 2 
My sin g |x" Mi cost p + 3% Mi cost —M32(*41 M? 4 > 1) cost p —1 | 
bye || 2 === || 2 / 2711/2 == (1). 
( 5 M? cos? @ + 1] (2 5 Mj cos? p + 1) +(- M? cos @ sing) i 


Die Auflésung des im Zahler in den geschweiften Klammern stehenden Ausdruckes nach cos? ) 
ergibt p*-Werte, die im gesamten M,-Bereich zwischen gy_, und @o liegen. Hieraus folgt, daB 
in einem isoenergetischen Strimungsfeld die Stelle hinter der StoBfront, wo das Produkt i 


sein Maximum erreicht, nicht identisch mit dem Ort ist, wo an der StoBfront gerade Schall- 
geschwindigkeit herrscht. An dieser Stelle erreicht lediglich der 0w-Wert, die Stromdichte dieser 
Stromlinie an der StoB®front das Maximum. Fiir einige Machzahlen der Anstrémung M, sind die 
berechneten y*-Werte in Tabelle 3 angegeben. 
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9. Erster und zweiter Differentialquotient der Geschwindigkeit in Stromlinienrichtung, Aus 


(19a) und (29a) erhalt man fiir den Geschwindigkeitsgradienten in Stromlinienrichtung 0w/ds 
bzw. fiir d1n w/ds die Beziehungen 


aw 1 w 
9s RMP Se} My). | 
33) 
dlnw 1 1 ( 
Fe = RTM). | 


Hierbei ist (vgl. Gl. [63] des Abschnittes 11) 


aa) as 
1 9 eas M? cos? pm + 1) (> Mi cost —“—*) 


VS 2S 2S 
2; 


pias —. 
Mi cos* p — x Mj cost — ~~ M$ cos? + 1 


Bezieht man die Strémungsgeschwindigkeit w auf die Schallgeschwindigkeit der Anstrémung a,, 
so erhalt man mit Gleichung (51) des Abschnittes 11 fiir den Geschwindigkeitsgradienten in Strom- 
linienrichtung folgenden Ausdruck: 


Ow Cd eee 
aoe OR (; Sa (34) 
= M? cos?p+ 1) (~My cost —“—*} (2%+-1) Mj cos*p —3 nye cos*p-+2 M? cost — M?—1 
= 2 2 
M, cos p (Mj cos? p — 1) (ee M} cos? + 1) oP = : M3 cospsing| | 
Fiir die Symmetrieachse vereinfacht sich dieser Ausdruck zu 
ene aa 
CRAM CR ge asl (35) 
as oo Rs x+] M, 


Die Geschwindigkeit hinter der StoBfront w ist durch die Gleichung (51) des Abschnittes 11 
festgelegt. Fiir eine konstante Machzahl der Anstrémung M, ist sie lediglich eine Funktion des 
Neigungswinkels der Stoffront an dieser Stelle oder, anders geschrieben, des Komplement- 
winkels zur StoSfrontneigung gy (p = 90° — y). Mithin ist auch der zweite Differentialquotient 
d?w/dq? an dieser Stelle der StoBfront bekannt. Andererseits gilt nach den allgemeinen Regeln 
der partiellen Differentiation 

Pw  dw/ds\? Pw ds dn , Pw /dn\? dw es , dw On 
dg? ds? ae dsdndp dp | On (55) ' ds 0g? ' On Og?” 


Aus den in Abschnitt 11 angegebenen Zusammenhangen erhalt man [vgl. (58) und (59) | 


di : dilnw 
FERS (op + 6) = sin (p+ 8) (ete » — de ). 


(36) 


(38 a) 


d d\n w 
Ju 008 (p+ 8) = — 008 (9 + 8) (4 tpt ) 
mit 
Gee aor 1 
2 
A= *x—1 a 
5 Mj cos*p + 1 
so daB mit (24a) und (24b) 
d’s : ( d\nw dR 
= 3 
ian sin (y + 0) [R( ete ~ Ie ae ; (37a) 
dn d\nw dR| 
)~‘e 


faa — 008 (P+ 0)|R(4te 9 bE 
gilt 


Fir den Kriimmungshalbmesser Q,, einer analytischen Kurve im cartesischen Koordinaten- 
system kann geschrieben werden (Hiitte I, S. 133) 
peel wey Re 


C., Via 


15* 
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Die Anderung von g,, auf der Strecke dx lautet 


dOp. Sy Ceey Me ye va ¥(2) 
Ges Vie é 
=R 
gleichzeitig ist mit y= are ctg y’ ils 
dy = yO z 
dx 12a y2 2 
woraus mit 9, = — R AE: 
dR 2. R 3 y/2y’—_(l+ y?)y’” Kriimmung der StoBfront. 
dp os ; 


Aus Schlierenaufnahmen von Verdichtungsstéen geht hervor, daB der beim abgelésten Ver- | 
dichtungssto8 auftretende StoSfrontverlauf im wesentlichen einen hyperbelahnlichen Charakter _ 


aufweist. Fiir eine Hyperbel = == i erhalt man aber nach kurzer Rechnung fiir den 
Kriimmungshalbmesser Ry an einer Stelle x, y den Ausdruck 


poe CR Uh 3 4 oll 
oder mit R, = h2/o als Scheitelkriimmungshalbmesser | 
pes Rolie G cs 1) + 1". (39a) 
Fiir die Anderung des Kriimmungshalbmessers Ry bei einer Winkeldrehung dg gilt somit 
i= 3 Ru “y(t a) (40) 
T= 3 R, ae ia 1) ! il etal (2 ! 2 (40a) 


Hierbei sind x, y die laufenden Koordinaten eines Hyperbelpunktes, wobei der Koordinaten- 


anfangspunkt im Kreuzungspunkt der Asymptoten liegt, und a, bdie in x- bzw. y-Richtung | 


fallenden Halbachsen der Hyperbel. 

Da die StoBfrontneigung y in groBer Entfernung vom Kérper mit dem Machschen Winkel 
der Anstrémung zusammenfallt, kann man bei einer hyperbelahnlichen StoBfrontform das 
Achsenverhaltnis b/a durch die Machzahl der Anstrémung ausdriicken. Es gilt! 

1 b SIN Yo J 


sin Yy> = — > SS i = — 
Yo M, = Ora a5 vo yM?—1 


2 
Mit tgy = ctgg= 2 = kann fiir die Kriimmungsanderung “ae einer Hyperbel auch 


dR y? b2 y2 7 
Gp = 8Buteo(se+ 1)(Bt 1) = 3Ruteo (Fe +1) (gay tI) 
oder 
dRy tgp (y¥* ie tgp 4s 
Fe Raga & te 1) = 3Ry SE Ai (40b) | 
mit 


geschrieben werden (vgl. Tabelle 4). 
Die Gleichung (40b) wiirde exakt gelten, wenn die StoBfront eine Hyperbel mit dem festen 


| 


durch die Machzahl der Anstrémung festgelegten Achsenverhiltnis ware. In praktischen Fallen 


ist jedoch mit Abweichungen zu rechnen, da die StoBfront sich der durch die Kérperform be- | 


dingten Geschwindigkeitsverteilung im Strémungsfeld anpassen muB. Die Schwierigkeit liegt 
dabei in der genauen Erfassung des Klammerausdruckes (y?/h? + 1), da ja dieser bei einer durch 
die Kérperumstrémung bedingten Anderung der Stoffrontform nicht mehr dem fiir die Hyperbel 
ermittelten Zahlenwert gleich sein kann. 


1 Vg]. FufSnote 1 von §. 211. 
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Tabelle 4, 4}, der Schallgeschwindigkeitsstromlinie als Funktion von M,. 


(tg?@) m=1 | 0,04767 | 0,09063 | 0,12867 | 0,16198 | 0,21521 | 0,27698 | 0,29600 

(y?/B) w=1 | 0,86941 | 0,75924 | 0,63126 | 0,58265 | 0,45327 | 0,23466 | 0,10947 

(Apo ea 9) 1,c604 || 157592 We6sis. 1.5827 | 14533 | 1.2847 | 1.1085 
M, 2,5 3,0 4,0 5,0 1,5 10,0 

(tg?~) m—1 | 0,26839 | 0,24425 | 0,21388 | 0,19832 | 0,18203 | 0,17606 

(y?/b2) m—1 | 0,05388 | 0,03149 | 0,01446 | 0,00833 | 0,00331 | 0,00178 

(42) 47-1 | 1.0539 | 1,0315 | 1,0145 | 1,0083 | 1,0033 | 1,0018 


Die GréBe des hierbei méglichen Fehlers kann abgeschatzt werden. Wie Tabelle 4 zeigt, 
nahert sich der Ausdruck: 4 = y?/b?+- 1 nur im ,,Unterschallbereich*‘ der StoBfront bei sehr 
kleiner Machzah] M, der Anstrémung im Schallgeschwindigkeitspunkt dem Zahlenwert 2,0. 
Sowohl mit zunehmender Machzahl M, wie auch mit abnehmendem @ fallt er aber schnell auf 
den Wert 1,0 ab. Ein kleiner Fehler in der Abschaétzung von y/b wirkt sich daher auf den Wert 
Ay nicht sehr stark aus. Wiirde man z. B. bei der Abschatzung des (y/b)y_-Wertes im Schall- 
geschwindigkeitspunkt bei der Anstrémmachzahl M,= 2,0 einen Fehler von 10° machen, so 
wiirde der Wert Jj nur um 2,07 bis 1,90% gefalscht werden. Aus diesen Erwagungen heraus 
wird man also fiir den ,,Unterschallbereich‘‘ der StoBfront bei der Ermittlung der d?s/d y- 
und d?n/dq?-Werte die Gleichung (40b) verwenden dirfen. 

Mit (40b) nehmen die Gleichungen (37a) und (38a) folgende Form an: 


d : d\n w 3t 
ig = Rsin (p+ 6) (cts p ae 5 sais din) 5 (37b) 
a dinw 3tgq@ 
T= Res (p+ 8)(4 tee 7 Sa (38b) 
mit den Abkiirzungen 
= M? cos? — 1 y2 a2 
Ae iy ; Ai epee bas (eT) 


grt M? cos?p + 1 


Durch Einsetzen von (37b) und (38b) in (36) unter Beriicksichtigung von (33) erhalt man nach 
umfangreicher Zwischenrechnung 

1 ow 
w Os? 


= 73 L(y, M,). (41) 


Auf die Wiedergabe der umfangreichen Beziehung L(y, M,) mu an dieser Stelle verzichtet 


werden |. 
Fir die Geschwindigkeit w, auf der Symmetrieachse nimmt (41) folgende einfachere Form an: 


: ei a ( : ) if ° | (42) 
w, Os IRE Neg p= (m3 — (#54 i + 1) 
, ae 
hg ITS 2 (emt 7) (43) 
sz rales 7] M, (M? — 1) 


10, Relative Lage des Sto8frontscheitels zum Kérper, Unter der Voraussetzung, dafi die 
- Geschwindigkeitsverteilung auf der Symmetrieachse lediglich durch die Geschwindigkeitsverhalt- 
nisse und den Gaszustand im Anstrémgebiet sowie durch den Scheitelkriimmungsradius der StoB- 
front bestimmt ist, kann die Entfernung des StoBfrontscheitels vom Kérper auf der Symmetrie- 
achse wie folgt abgeschatzt werden: 


1 In der Originalarbeit ist die Funktion L(y, M,) unter Beriicksichtigung der Stofverluste an- 
gegeben. Jedoch ist die Funktion nicht tabuliert worden. 
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Fiir die Geschwindigkeit w, auf der Symmetrieachse in der nachsten Umgebung des StoB- | 
frontscheitels gilt die Entwicklung 


2 f 
Ow. % 02 w, x2 


w,(x) = w,(0) 4 aoc | an or nay 


wenn man den Koordinatenursprung in den Schnittpunkt der StoBfront mit der Symmetrieachse 
verlegt. 

Bezieht man die Strémungsgeschwindigkeit w, auf die Schallgeschwindigkeit a, des Anstrém- 
gebietes, so ergibt sich aus Gleichung (51) des Abschnittes 11 auf der Symmetrieachse 


x—l,,. 
ty \ ial Laie he 
( e+) M, ; 
und man erhalt mit den Gleichungen (35) und (43) 


a; 


(51a) | 


\2 
x—l,,. RL aed py os 
a oe es Be ( 2 ae A 2 « 2 ca 2 a 
x 


Th vec lilels aM x+1 M, R, %+-1} M, (M?—1) 2R? | 


ay 


Da dem Schnittpunkt der Kérperkontur mit der Symmetrieachse ein Staupunkt in der Stré- | 
mung entspricht, gilt fiir die Entfernung x, der Sto8front vom Kérper 


ee Ne 
‘ wm —*>3) . 
sey es 2 of gingge kia Ee penne? | a: z A EV 
2 ie x+1 ; 2 7R, \e+1 Mi —1 2 R? 
Man erhalt nach kurzer Rechnung 
M2+)1 
Soar, Mi eeal | ; 2 1 
nes 2 ngs Bea : : M?—1 I}. (44) 
Ab 
24 


Die Ergebnisse der Auswertung sind in der Tabelle 5 und in Abb. 5 wiedergegeben!. 


ay 11, Zusammenstellung einiger 
<< —$$ {> = 1 aus den StoBgleichungen folgen- 

den Beziehungen fiir den Ge- 
| schwindigkeitsvektor und die Gas- 
zustandsgréBen, Die Zustands- 
anderungen im Gas beim Durch- 
schreiten einer Sto®front kénnen 
sehr anschaulich mittels des 
Busemannschen StoBpolaren- 
diagramms und des Druckberges? — 
dargestellt werden. Der formel- | 
maBige Zusammenhang zwischen 
den einzelnen physikalischen und 
geometrischen Gréfen ist in iiber- | 


sichtlicher Weise bei F. Schubert? 


G3 


Oe 


1M, 


10 09 08 G7 06 O5 OY O3 G2 g7 0 M, 
angegeben. 
10 105 11 1512 «13 «14 «15 16 = 18 = 20 80535 80 40 50 1510 20 ~ Man erhilt mit d : d | 
; en in er 
Abb. 5. Verhiltnis x,/R der Entfernung’x," des’ StoBfrontscheitels vom Staupunkt Arbo: Ss b 
des umstrémten Kérpers zum Scheitelkriimmungsradius Rs der StoBfront. Die ge- rbeit von chu ert angegebenen 
strichelte Kurve entspricht den abweichenden Ergebnissen von Dugundji. Beziehungen, wenn man an Stelle 


der StoBfrontneigung y den Komplementwinkel y — 90 — y einfiihrt, folgende Ausdriicke fiir 


' Bei Umrechnung der Ergebnisse von J. Dugundji erhalt man fiir den zweiten Differentialquotienten 
der Geschwindigkeit auf der Symmetrieachse ein gegeniiber hier abweichendes Ergebnis. Demzufolge er- 
geben sich auch fiir die Entfernung x, des StoBfrontscheitels vom Staupunkt andere Werte, fiir groBe Mach- 
sche Zahlen der Anstrémung sogar iiberhaupt keine reellen Werte (vgl. Auftragung in Abb. 6). Eine Nach- | 


priifung der von Dugundji abgeleiteten Bezichung ist schwierig, da ein Teil der Zwischenrechnung weg- 
gelassen ist. 


2 Vel. FuBnote 1 von S. 211. 
3 Vel. FuBnote 1 von S, 208, 
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1. den auf den Druck im ea ay P, bezogenen Druck p: 


Dre Se aS ll 2 *%—I1 
a = Mi cos? ~ — eT = gpa (*Micos# 9 —*5"), (45) 
2. die auf die Dichte im Anstrimgebiet 0, bezogene Dichte 0: 
0 i eee M? cos” YP 
Rand. (46) 


@ tee! ae a 


3. die auf die Temperatur im Rasen T, bezogene ea hee E: 


—] 
M? cos? yp — —.— 
lt x—l1 ¢ P 9 
Se gy | he oon 2 2 | 
T, ( 5 M? cos? y 4 1) ae ie ) (47) 
= M? cos? p 
4. den Ablenkungswinkel : der Stromlinien hinter der StoBfront: 
[m3 — —- OG cos os Mi eee M? 
Cr Ol —— ct 
g i es fay a 8° \ Mpcotg@=ia +)” (48) 
eee ae aes M, &? 
5. den Drosselfaktor K: 
x 1 2 poe x +1 1 
i T= Mteostp \ 75 Le a 
se ecea beeer le (49) 
1 5 M? cos?y + 1 x M? cos? o—- 
Tabelle 5. Relative Lage des Stoffrontscheitels zum Kérper (Gl. 44) Mit diesen Gleichungen kénnen 
ay 
sowie Zahlenwerte iad Gleichungen (51a), (35), (42), (44). also bei vorgegebener Machzahl 
Gl. (51a) Gl. (35) Gl. (42) (Gl. 44) der Anstrémung M, die Ab- 
x Gonchuindies Geschwindig- Gradient der Relative Lage lenkung 0) sowie die einzelnen 
Anstrém- ese 8 keits- Geschwindig- | des StoBSfront- Zustandsgrifien des Gases un- 
Machzahl gradient keitsinderung scheitels 4 
M, — mittelbar hinter der Stoffront 
w/a, See eae eee Re Xo fiir jede Stoffrontneigung 
w, Os De Ge 2 R, 
= 90 — 
1,0 1,00 —0,83160 — —_ vy 7) 
1,05 0,96882 —0,91539 | —10,0002 0,36374 levoht Geetimintverdent 
1,10 0,94124 —1,00027 — 5,9319 0,43602 
aS 0,91627 —1,08598 — 4,6362 0.46308 
ie2 0,89508 —1,17226 — 4,0339 0.47114 
11683 0,85861 —1,34567 — 3,5225 0,46280 
1,4 0,82976 —1,51884 — 3,3567 0,44226 
iS 0,80700 —1,69034 — 3,3273 0,41805 
1,6 0.78919 —1,85900 — 3,3637 0,39603 
1,8 0,76512 —2,18418 — 3,5271 0,35569 
2,0 0,75260 —2,48906 — 3,7379 0.32384 
2520 0.74850 —2,83762 — 4,1140 0,29207 
BAS 0,75364 —3,14884 — 4,2788 0,26857 
2,75 0.76550 —3,42403 — 4,5224 0,25059 
3,0 0,78240 —3,66597 — 4,7416 0,23658 
3,5 0,82700 —4,06378 — 5,1093 0,21659 
4,0 0,88150 —4,36934 — 5,3964 0,20334 
5,0 1,00832 —4,79036 — ,9,1967 0,18749 
eo 1,37387 —5,29062 — 6,2773 0,17155 
10,0 Oli —5,43043 — 6,4705 0,16592 
20,0 3,40956 —5,69746 — 6,6712 
oo 0,158616 
* Bei der Ermittlung der Zahlenwerte wurde mit x = 1,405 BUGS ete polare mech. Busemony: 
gerechnet. 


Weiterhin erhalt man aus den Gleichungen (45) und (46) fiir die auf die Schallgeschwindigkeit 
im Anstrémgebiet a, bezogene lokale Schallgeschwindigkeit a die Beziehung 


a 1 s¢—— | “x—1] 
= M? cos2 1) (4 M2 cos? mw — ; 50 
a, eee \( 9 1 cos e+ ) [x 1 COs" Pp 5) . (50) 
2 
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Beriicksichtigt man, daB fiir eine durch die StoBfront hindurchgehende Stromlinie die normal | 
zur StoBfront verlaufenden Komponenten Wp und w, der Geschwindigkeitsvektoren vor und 
hinter der StoBfront der Beziehung! 

women Os 

Mi Q 
geniigen miissen, so folgt mit Gleichung (46) fiir die Strémungsgeschwindigkeit w unmittelbar | 
hinter der StoBfront 


— 2 
— M2? cos? y + 1) 
w? —= Ww; | Wy = Ww? sin? yp — w? cos? p fa) = wi sin? y — wy cos? p-—_ ary ats | 
e | 5 M? cos?* 0) | 


Bezieht man die Strémungsgeschwindigkeit w auf die Schallgeschwindigkeit im Anstrém-_ 
gebiet a, und beriicksichtigt dabei, daB w,/a,— M, ist, so erhalt man fiir die Strémungsgeschwin- | 
digkeit w/a, 


w 1 jee 20 m+1 Dale 
Pera pe ae M? cos? 1) + M? cos sin (51) | 
Zs E d M, cos Pp | ( 2 ‘ aie ( 2 u a 9) 
sowie fiir die beiden w-Komponenten in Normal- bzw. Tangentialrichtung 
“= * Mi cost +1 
iy i : (52) 
“1 z 2 M, cos 
ww; . 
ail M, sing (53) | 
bzw. in x- und y-Richtung 
w, Micos*o—1 ., alga | 2 : 2 ] 
apn ebay sin p = 77 te Y (Micos Pit ares (54) 
2 4 YCOS Pp 
PES Y pa ee lee (M, cos? p— 1) Eee (55) 
a, aerate pata i M,~ 


Beriicksichtigt man, da fiir den Leitstrahl r der Busemannschen StofBpolaren (vgl. Abb. 6) 
mit (46) und. M, = w,/a, | 


x— 1 


M? cos* p +1 


Bee cos? p (56) 


r= Wp, 


— Wy, = W, COS YP — W, COS | 


2 ety ae. pe Poe ee a? 
x+ 11 “ x+1M?cstp x+1 01 COSP —~ Tt cos P 


geschrieben werden kann, so gilt fiir die Strémungsgeschwindigkeit w 


w = ]/(r sin y)? + (w, — rcos g)? = \w? =-\1? — 2710; cos (57) | 
AnschlieBend werden einige bei der Rechnung oft benutzten Ausdriicke, die sich aus den 


yes a (45) bis (55) leicht herleiten lassen, iibersichtlich zusammengestellt. Es ist (vgl. 


x 5 d cos p sin 


sin (p + 6) = retail 


Ww = 2 lane 
VC 5 1 cos? p +1) le = M} cos? 9 sin? 


2 (58) 


x—l1 
cos (p+ 6) = “»— 2 
es (a M? cos2 1) (* 1)" Mj 2 +5 
~~ 9 1 COS* Pp -- a 5) 1 COS? Pp sin? p 
1 Vgl. FuBnofe 1 S. 211. 


M? cos* p+ 1 


- (59) 
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= 2 M? cos p sin y 
tg (p+ 6) = — | 
5 M? cos? y + 1 
Der Ablenkungswinkel 6 wird 
2“ M? — M? cos? 9 +1 
6 = are ctg acne ctg @|, 
daraus 
dé x Mj cost p — “= * My cost p+ 2 M cost —“** Mp —1 
Tap 2 
ie eae M}%, cos? p — # Mj cos* gm + (% — 1) Mi cos? y+ 1 
Die értliche Machsche Zahl M hinter der Stoffront ist 
= 2 2 
ee Pe (“5+ my cos? y + 1) + (4+ My cos psin og} 
i (2 Mi cos? p + 1) (x M3 cos? p —“—*} 
daraus 


uw 


Rees 
: a tt My cost —x Mj cosy —"—™ M3 cost y +1 
M?— |= 


2 x— 1 
IS 
Aus Gleichung (45) erhalt man ferner 


== Il 
M3 cos? p +1] (. IE cos? p —*—*| 


az 
dinp Py Di a5 2 % M? cos? p 
oak pee % MY cost p —“—* 


Die Winkel sind in Neugrad gemessen, daher ist y= 100 — ¢8. 


Eingegangen am 4. November 1950. 


Anschrift des Verfassers: Harald Melkus, Wright-Field MCJ-AXF Dayton Ohio USA. 
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Membrantheorie verallgemeinerter Rotationsschalen. 
Von W. Zerna. 


1. Einleitung. Die Bestimmung des Membranspannungszustandes beliebiger Schalen laBt sich 
auf die Lésung einer partiellen Differentialgleichung zweiter Ordnung fiir eine Spannungsfunktion 
zuriickzufithren, aus der dann die Schnittkrafte als zweite Differentialquotienten folgen’. In 
vielen Fallen ist die Berechnung von Membranschalen auf diesem Wege jedoch nur schwerlich 
durchfiihrbar, da sich die Randbedingungen fiir die Spannungsfunktion nicht immer ohne weiteres 
angeben lassen, und aufSferdem erscheint der Umweg iiber eine Spannungsfunktion, die selbst ja 
gar nicht interessiert, sondern nur ihre zweiten Ableitungen, unzweckmaBig, insbesondere im 
Hinblick auf die Anwendung von Naherungsrechnungen, wie z. B. das Differenzenverfahren, auf 
die bei allen nicht ganz einfachen Problemen zuriickgegriffen werden muB. 

Fiir eine spezielle Flachenklasse, die aber eine groBe Mannigfaltigkeit verschiedenster Flachen 
enthalt und verallgemeinerte Rotationsflachen genannt werden midge, soll im folgenden ein 
Lésungsweg angegeben werden, der die Einfithrung einer Spannungsfunktion vermeidet und fiir 
die praktische Rechnung besonders geeignet erscheint. 


2. Geometrische Beziehungen. Es werde ausgegangen von einer ebenen Kurve K, deren Orts- | 


vektor a sich in einem durch dieEinheitsvektoren i,, i, festgelegten kartesischen Koordinaten- 
system x, y mittels 
a= x(9)i+ ¥(P)t, (1) 


darstellen 14Bt, wenn ? einen 


dem Punkt von K werde nun 
ein Vektor 

ri=r(8)[E(~)ts+7(P)te] (2) 
abgetragen, worin r(@) eine 
beliebige Funktion des Para- 
meters #, sowie &(p), 1 (p) be- 


ten Parameters p und i, einen 
Abb. 1. Verallgemeinerte Rotationsfliche mit Breitenkreisschnitt. auf a und te senkrecht stehen- 
den Einheitsvektor bezeichnet. 


t= a(d) + r(d) Ey) (3) 


ist dann eine Flache bestimmt (Abb. 1). 


Durch 


Wird die Kurve K als Gerade gewahlt und beschreibt rf Kreise, deren Ebenen senkrecht zu | 
K stehen, so geht die Flache in eine Rotationsflache iiber. Es werde deshalb die durch (3) dar- 


gestellte Flache verallgemeinerte Rotationsflache genannt. Die Kurven # = const. entsprechen 
den Breitenkreisen und die Kurven y = konst. den Meridianen. 
Fiir r= const. ergibt sich eine Translationsflache mit K als Erzeugende. Ist K eine beliebige 


Gerade, so lassen sich Kegelflachen gewinnen, wie iiberhaupt die meisten praktisch inter- | 


essierenden Flachen erfaBt sind, was sich leicht einsehen laBt. 

Ableitungen nach # werden durch einen Strich, Ableitungen nach y durch einen Punkt an- 
gezeigt. 

Die Basisvektoren der Flache folgen aus (1) bis (3) zu 


gv tty iet rt, (4) 

a Se (5) | 
wenn 

PN Sy ie (6) 


1 Vgl. A. Pucher, Beton u. Eisen. 33 (1934), S. 298. 


Parameter bezeichnet. In je- | 


liebige Funktionen eines zwei- | 
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Die Einheitsvektoren i,, i, lassen sich durch i und f ausdriicken: es ist 
: i ‘ 
b= > Gi ep, | 
eae ye) | 


mit 
D= én — &y (8) 
Mit (7) wird aus (4) 
Wi / § / se 
a= wit y 5 it(r yp)! (9) 


/ iy Be 8 uv ee 
aaeity pit(s vp)t | 


e=e=r'i, : (10) 
es= >, (Di+ Ht), | 
mit 
H= En — Ey 7 (11) 


Die Koeffizienten der ersten Fundamentalform der Flache (die Komponenten des MaBtensors) 
errechnen sich aus (5) und (9), wie folgt: 


1 — 8° = x2 1 yy’? 1 r/2(€2 1 2) 1 2 x’ yn , | 


Boo Cg'Cg—= 7(E2-+ 77), (12) 
Sig CyCg—= rr’ (EE + ny) +7ry'7 . | 
Die Determinante des MaBtensors ergibt sich damit zu 
& = [Sin] = 811 822 — Bip » | 
$= (E29?) [x? 4 2 2 r'y'g + (E+ 9?)] (13) 


ee a Ape a ES 2H = 7)? 
3. Gleichgewichtsbedingungen. Die Langskrafte mégen mit Ny und N,, die Schubkrafte mit 
T bezeichnet werden. Die Belastung der Schale sei durch den Belastungsvektor in der Form 
P= Petit Ppit prt (14) 


gegeben. Es werden reduzierte Laingskrafte und Belastungskomponenten eingefiihrt, indem 
gesetzt wird 


No | Ne ie NA (15) 
ll ! $20 
Pe=J8Pz> Pyo=VEPp> Pr=VEPr- (16) 


Mit Verwendung der definierten reduzierten GréBen lassen sich die allgemeinen Gleichgewichts- 
bedingungen der Membrantheorie in der Form darstellen! 


ex(No + T)+ eo(T’ + Nz) + e\No+ 20e,T+e,N, +p g=0. 
Einfiihrung von (5), (9), (10), (16) liefert die Gleichgewichtsbedingungen 
Mt+T+5N4+2=0, 


x 


PN = 


N,+6N,+ P,=0. J 


(17) 


D = = 
PN cee Ny 0, 


1 Vgl. W. Zerna, Ing.-Arch, 17 (1949), 5. 149. 
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Darin bedeuten 


| ; 
J 


(19) | 


In den Gleichungen (17) sind als Sonderfall die Gleichgewichtsbedingungen der Rotations- | 
schale in der vom Verfasser! angegebenen Form enthalten. 


4. Differentialgleichung des Problems. Ahnlich wie bei den Rotationsschalen lassen sich die _ 
drei Gleichgewichtsbedingungen (17) auf eine einzige Differentialgleichung zuriickfiihren. Zu 
diesem Zweck werden aus den ersten beiden Gleichungen die Ableitungen der Schubkrafte elimi- 
niert. Wird dann aus der ersten und letzten Gleichung 


{Ips hye a Nips ei (20) 
N,=— GN, — P, (21) 1] 
beriicksichtigt und zur Abkiirzung gesetzt | 
Re x” ar’ 
PE re 
erner la g 
x oe r’ x’ 
pore (Shae 
D: { D: . 9 y Px ee , 
a= Py— Fi +5 i +(5) BtAB ae (23) 
so folgt 
NS ONG GING + BYNUM eg 0. (24) 
Dies ist eine partielle Differentialgleichung zweiter Ordnung fiir die reduzierte Schnitt- 
kraft Ny. Ist sie gefunden, so folgen die Schubkrafte T aus (20) und N, aus (21), 
Die Differentialgleichung (24) laBt sich noch etwas vereinfachen. Einfithrung von 
D— rj x’ Ng (25) 
und 
q=—r jx’ q: (26) | 
ws ; xii’ Tes r’ 3 xi! raed ne 
J=F | ax | ane Ae) = (27) | 
bringt 
DO CDE KC I TD ge C8) | 


womit die das hier betrachtete Problem beschreibende Differentialgleichung gefunden ist. 
Sind am Rande der Schale die Liangskrafte N, oder Ny vorgeschrieben, wie dies meistens | 

der Fall ist, so sind bei Beachtung von (25), sowie der letzten Gleichung (17) auch sofort die Rand- 

werte von ® bekannt. Damit lat sich dann (28) mit Hilfe von Naherungsrechnungen wie z. B. 


das Differenzenverfahren ohne Schwierigkeiten behandeln. 


(Eingegangen am 30. November 1950.) 
Anschrift des Verfassers: Dr.-Ing. W. Zerna, Hannover, Technische Hochschule. 


1 Vgl. W. Zerna, Ing.-Arch. 17 (1949), S. 223. 
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___ lischen Eigenschaften der | 

< zeichnis. | Fe piin wenu 
No st (i pan te : ‘ F 

_ Aus den Besprechungen: EE tay Ace ah & 

pier ese Wer ass, Buch durchgearb rd é ¢ des allgemeinen ‘Wissens das Einzel- 

__ ergebnis vielfach in ein Ganzes einordnet. So ist das Buch ein Wegweiser zu einer allgemeinen wissenschaftlichen | 


Y 


_ immer von der Gefahr zunehmenden Spezialistentums infolge der Entwicklung der Naturwissenschaft gesprochen 


% wird, so sollte man diese auf Grund der Entwicklung der Grundwissenschaften sich herausbildenden ‘Zusammen! 


a 
WwW 
<= 


5 ler Physiker, ‘seine Theorien am Verhalten der Metalle zu iiberpriifen und sie fir das ‘ ferstiindnis. der Metall-— 
4 SPAN etre UNCP Oe ee 


_ Hochwertiges GuBeisen, 
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__ Metallurgie seiner Herstellung. 


a. 


Von Dr. 


-—__verbesserte Auflagé Mit 1063 Abbildungen. XIT, 1070 Seiten. 1951, Ganzleinen DM 135.— 
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¢- Inhaltsibersicht: I. Einleitung. — II. Die konstitutionellen Grundlagen der Eisen-Kohlenstoff-Legierungen. 
-— III. Uber molekulare Eigenheiten kohlenstoffhaltiger Lésungen. — 1V. Der Einflu8 des Siliziums auf die Gleich- 
gewichts- und Graphitisierungsvorginge. — V. Der Mechanismus der Graphitisierung siliziumhaltigen GuBeisens. — 
- VI. Die strukturelle Beherrschung der metallischen Grundmasse. — VII. Die Primirkristallisation des GuBeisens. — 
' VIII. Der EinfluB der stindigen Eisenbegleiter. — IX. Der Einflu8 der Gase. — X. Technologische Higenschaften 
_  __. fliissiger und erstarrender Eisen-Kohlénstoff-Legierungen. — XI, Die mechanischen und elastisehen Eigenschaften 
des Gufseisens. — XII. Die physikalischen Eigenschaften des Gufeisens. — XIII. Verhalten des GuBeisens bei hohen 
-_ und tiefen Temperaturen. — XIV. Technologische Eigenschaften des festen GuBeisens. — XV. Die chemischen Eigen- 
schaften des GuBeisens. — XVI. Der zusitzliche Oberflachenschutz von Gufecisen. — XVII. Festigkeitseigenschaften 

_ von GrauguS und Temperguf nach Vorkorrosion, — XVIII. Der HinfluB thermischer Nachbehandlungen auf die Ge- 
fiigeinderungen und die Higenschaften perlitischer Grundmassen (Hartungs- und Gliihprozesse), — XIX. Legiertes 


. GuBeisen. —X X. Die spanlose Verformung des GuBeisens. — XXI. Das Schweifgen von GuSeisen.— XXII. Fehler- 
-quellen bei der Gefiigeuntersuchung von GuBeisen. — XXIII. Das Schmelzen von Gufeisen im Kupolofen. — XXIV. 
_ Andere Schmelzéfen.— XXV. Einige besonders wichtige Anwendungsgebiete fir GuBeisen. — XXVI, GieBen und 
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SchweiBen. — XXVII. Anhang. — Sachverzeichnis. 
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Das einschlagige in- und auslindische Schrifttum, dessen Studium dem Fachmann ein abgerundetes Bild vom 
Stand der Technik hinsichtlich Herstellung, Higenschaften und Hinsatzméglichkeiten des heutigen hochwertigen 
GufBeisens zu vermitteln yermag, ist so tiberaus umfangreich geworden, daf nur wenige Fachleute Zeit und Mu®8e 

finden darften, davon in vollem Umfange Gebrauch zu machen. Hier setzt der Zweck des vorliegenden Werkes ein. Es 


Beziehungen zwischen dem Gefiigebau und den Higenschaften deselben sowie tiber den zweckmaBigen Kinsatz dieses 
3  Werkstoffes und seine kiinftigen Entwicklungsméglichkeiten. Die Darstellung und kritische Betrachtung des Fach- 
 . gebietes erfolgte in einer Form, die den Leser immer wieder auf offene Probleme hinweist und ihn immer wieder zu 
eigenem produktiven Denken anregt, um auf diese Weise der weiteren Entwicklung dieses Werkstoffes zu dienen. 


‘Uber die erste, 1942 erschienene Auflage schrieb die Zeitschrift des Vereins deutscher Ingenicure: 
,,Hs ist ein Verdienst Piwowarskys, daB ‘er mit dem vorliegenden Werk cine sehr umfassende, auf den letzten 


Stand gebrachte, kritisch wiirdigende und zur weiteren Mitarbeit anregende Zusammenfassung des gesamten Wissens~- 
sbietes gibt, das sich bislang in zahllosen Einzelveréffentlichungen zerstreut fand.“* 
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_ Grundlage fiir die Heranbildung der Jugend, in der es keine Trennung zwischen Leichtmetall-, Schwermetall- und — 
ee Eisenfachmann gibt. _Auch hier verwischen sich die Grenzen, wie sich in der modernen Naturwissenschaft die Grenzen, — 
_ zwischen den groBen Wissenschaftsgebieten der Physik, Chemie und Biologie zu verwischen beginnen. Wenn _ 


_ fassungen als groBen Gewinn nicht tibersehen. Diesen Gewinn veranschaulicht das Buch von Masing; der Metall- — 
_ kundler wird dadurch angeregt, sich immer wieder mit den erforderlichen Grundwissenschaften zu beschiftigen, — 


Se _legierungen und zur Schaffung neuer Legierungen nutzbar zu machen... Ee SN ausieet poe Ve 
aT ae WES oe ee) cna SE aie ees “ Eduard Houdremont in Stahl und Eisen‘*. 
os ES OSS toh Rag: Sa ae ea aL ieee ar clr ae a 5h Bes Se ee 

ie Le . : f= Net Sg ae oS 


1er ste Von -Ing. habil. Eugen Piwowarsky, ord. Professor der 
__ Eisenhiittenkunde, Direktor des Instituts fiir allgemeine Metallkunde und das gesamte 
> GieBereiwesen der Rheinisch-Westfalischen Technischen Hochschule Aachen. Zweite, — 


soll dem Fachmann in geschlossener Buchform einen Uberblick bringen iiber die Metallurgie des GuBeisens, tiber die _ 


ri 


Das Elektrostahlverfahren. Ofenban, Elektrotechnik, Metallurgie und Wirt-_ nae | 


schaftliches. Nach F..T. Sisco ,,The Manufacture of Electric Steel‘. Zweite deutsche, 
erweiterte Auflage. Von Dr.-Ing. Heinz Siegel, Diisseldorf-Oberkassel. Mit 140 Ab- 


*  bildungen. XI, 432 Seiten, 1951. ~ Ganzleinen DM 31.50 | 
Inhaltstibersicht: A. Einfiihrende Betrachtungen? Die Entwicklung des Hieibroutablvertalinens: Die Licht- © _ 
bogendfen. Die Induktionséfen, — B. Der Lichtbogenofen: Die bauliche Gestaltung der Lichtbogendfen. Die elek- =] 
trische Ausriistung der Lichtbogenéfen. Die Elektroden. Die Energiewirtschaft der Elektrostahléfen- — €. Der =| 


kernlose Induktionsofen. — D. Verschiedene Ofen. — E. Gesichtspunkte flr den Bau neuer Elektrostahlwerke. — jae 
F. Die feuerfesten Baustoffe fiir den Elektrostahlofen. — G. Die Einsatzstoffe und die Schlackenbildner. — H. Die 
Metallurgie des Lichthogenofens: Die allgemeine Schmelzungsfihrung bei festem Papseys, Die Kockvorginge beim - 
basischen Verfahren. Die Feinungsvorginge beim basischen Verfahren. Die Abarten des Einschmelzens. Die Ab- ; a 
arten des Kochens. Der Betrieb beifliissigem Einsatz. Die Besonderheiten bei der Herstellung einzelner Stahlsorten. 
Saurer Elektrostahl. — J. Die Metallurgie des kernlosen Induktionsofens. — K. Gegeniiberstellung des Betriebes - 
der Lichtbogenéfen und der kerniosen Induktionséfen. Weitere Entwicklung. — L. Selbstkostenwesen im Elektro- 
stahlbetrieb. — Sachverzeichnis. os - . on 


Der ,,Sisco-Kriz‘‘*, die erste Darstellung des Elektrostahlverfahrens in deutscher Sprache, geht zuriick auf das 
‘Buch ,,The Manufacture of Electric Steel‘‘ von F. T. ‘Sisco, das den Vorzug hatte, erstmalig den vorliegenden Gegen- 
stand in auBerordentlich gliicklicher Weise zu behandeln, weshalb damals die stoffliche Gliederung und die Betrach- = 
tungsweise in die deutsche Bearbeitung iubernommen worden sind. Da auch heute noch die angewandte Art der ee 
Darstellung als vorbildlich bezeichnet werden kann, so wurde sie auch in der Neuauflage beibehalten. Neue Ergeb- ee | 
nisse und Erfahrungen muBten allerdings auch in anderer Form zur Darstellung kommen, Da Dr, Siegel in dem < 
von Dr. Kriz ausgerichteten Betrieb lange Jahre tatig gewesen ist, war es ihm méglich, die erforderliche Neubearbei- ts 
tung vorzunehmen. So gibt das Buch auch in der zweiten Auflage einen klaren und einfachen Uberblick iiber das 5 
gesamte Betriebsverfahren und ist im iibrigen so abgestimmt, da es auch dem reinen Praktiker etwas bedeuten 


wird. Nach wie vor bleibt es das Ziel des Buches, den Leser zu unterrichten, in wéelcher Weise einwandfreier Elektro- 
ra 


stahl wirtschaftlich erschmolzen werden kann. ~ Zi : nes 


Der Niederfrequenzofen wurde ausfihrlich behandelt, obwohl er heute nicht mehr iiberall anzutreffen ist. Es a. 
wurde dabei einmal von dem Gesichtspunkt ausgegangen, da gerade in letzter Zeit yesukupeapies mit niedrigen 
Frequenzen erprobt worden sind und es also nicht ausgeschlossen ist, daB der Niederfrequenzofen in der einen oder 
anderen Form wieder auflebt. Ferner sind aber auch die gesammelten Erfahrungen zu wertvoll, als da man sie 
schon heute vergessen sollte, zumal in Qualititswerken von Ruf in Europa und Amerika immer noch Niederfrequenz- A fe 
éfen zur vollsten Zufriedenheit arbeiten, Aus ahnlichen Griinden wurde auch die Bodenbeheizung der piebebpeee: 12 
éfen noch erwahnt. ; 2 


Einfiihrung in die Probleme des elektrischen Lichtbogen- cud 


Widerstandsofens. Von Dr.- Ing. Erich Kluss. Mit 163 Abbildungen. VII, 263 Seiten. 
1951. Ganzleinen DM 28.50. 


Inhaltsiibersicht: Einleitung. — Der elektrische Stromkreis. Das Kreisdiagramm, — Die Induktivitat. Der ine 
induktive Widerstand, — Das magnetische Feld der Stréme. — Die Stromverdrangung. Wirbelstréme. — Kontakt- “k 
probleme. — Die elektrischen Verhiltnisse im Ofen. — Unsymmetrieprobleme. — Transformatoren. — Nick ikanientee 
Betriebszustinde, Ausgleichvorgange. — Schrifttum, Sachverzeichnis. ; z 


Der elektrische Ofen gehért in immer steigendem MaSe zu den GroSverbrauchern an elektrischer Energie. In 


ihm wird die zugefihrte Energie in thermische bzw. chemische umgesetzt. Die Umsetzung erfolgt | bei relativ geringen 
Spannungen mit entsprechend hohen Strémen, die bei GroBdfen in der Gré8enordnung von 100 000 A liegen. Die 


zweckentsprechende Fortleitung dieser hohen Wechselstréme zum Ofen, die durch ihre Wechselfelder bedingten Bs 
Erscheinungen der Induktion, Wirbelstréme, Stromverdrangungseffekte bilden die Kernprobleme fiir den Bau und eee 
Betrieb derartiger Ofen. Sie werden hier ausfiihrlich behandelt. Dariiber hinaus werden die elektrischen Vorginge =. -} : 
im Lichtbogen bei Lichtbogenéfen, sowie die Verhiltnisse im Strémungsfeld ‘bei Widerstandséfen gebracht, Beim ae 
Betriehb von Drehstroméfen treten die Probleme der Unsymmetrie auf, wie z. B. die Erscheinungen der »toten’* 
und ,,scharfen‘ Phase. Sie werden hier rechnerisch ausfihrlich behandelt-~GréBte Bedeutung fiir den Betrieb des mS 
Ofens gewinnt der Ofentransformator, Ersatzschaltungen desselben, sowie bei GroBéfen der Regeltransformator. 3 : 
Es ist deshalb hier den Ausfithrungen tiber den Transformator und seine Verwendung ein breiter Raum gelassen, A fe. 
Zum SchluB8 werden die nichtstationiren Vorgiinge im Ofenbetrieb, Schaltvorginge und ihre auewihesee, die i im ied 
Betrieb schon éfter unliebsame Uberraschungen brachten, behandelt. ae 
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